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FOREWORD 

During the suiraner of 1967, A Cooperative College-School Science 
Program in Computing Science Education was conducted by the University of 
Colorado Department of Civil Engineering with the support of the National 
Science Foundation- (Grant GW 1755). The program consisted of two 5-week 
terms. The first term was held in cooperation with the Boulder Valley 
School District; the second term, in cooperation with the Jefferson County 
School District. Twenty mathematics and/or science se^^ondary school 
teachers and twelve eleventh grade students- participated during each term. 

The course content for the two terms was identical and consisted 
of course work, laboratory work, and study periods each day. The course- 
work was composed of two- formal lecture courses — Computer Related 
Mathematics and Computer Programming; while the laboratory work encompassed 
demonstrations, discussion periods, and problem- solving by the participants 
using BASIC language on remote consoles connected to the General Electric 
time-sharing computer system in Phoenix, Arizona. 

During each term, many of the teachers (individually or in pairs), 
assisted in some cases by students, undertook special projects. These 
special projects consisted of writing up the details of typical high school 
lectures, incorporating computer- related concepts. In some cases, the topics 
chosen were those ordinarily - presented, with modification to reflect a 
computer oriented approach. In other cases, the topics were new computer 
related ones that' could replace certain elements in an existing course . 

This document is a compilation- of these special projects. 

The main objective of the special- projects was to begin to motivate 
the high school teacher's development toward computer oriented educational 
materials. Hopefully, by printing all projects, individual teachers will 
be able to review the ideas of the participants and be further motivated 
toward- incorporating- computer related mathematics and concepts into their 
classes . 



Boulder, Colorado C. C. Feng 

August, 1967 
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CALCULATING THE REAL ZEROS OF A REAL, POLYNOMIAL FUNCTION 

By Robert L. Bassett, Carol E. Friend, and Donald C. Lohrentz 

To find the real zeros of an Nth degree function there are five parts to 
this project: 

I. Finding the integer roots of a function. 

II. Searching for an interval where a root occurs. 

III. Finding the real zeros of a function, half interval method. 

IV. Finding the real zeros of a function, false position method. 

V. Finding the real zeros of a function, Newton’s method. 

Any one of the five parts may be used, with only a minimum of 
programming experience . 

Part I Finding the Integer Zeros of an Nth Degree Function 

Psrt I could be used as introduction to the other parts, since the 
second year algebra student is familiar with this process. Parts III, IV and 
V are more difficult mathematically, therefore graphs are used to explain 
them. Part I of this project will assume that the reader has already covered 
the following topics: 

A. The remainder theorem: The remainder when any polynomial P(x) is 

divided by (x-a) is P(a). 

B. The factor theorem: The factor theorem and its converse; (x-a) is a 

factor of p(x) if and only if P(x)=0 when x=a. 

C. The fundamental theorem of algebra: Every complex polynomial equation 

of degree n has at least one complex number as a root. 

Let’s assume that an Nth degree polynomial is given by: 
f(x) = a^^x^ + a^x^ ^ ^ ’ * ‘ + a x + a 



n 



n+1 



Problem: to find the integer values of x that make f(x)=0 or to find 

{xel 



n ^ n-1 ^ n-2 . 

^1^ + 3 ^,^ + a^x + 



+ a x + a - = 0} 
n n+1 



First set f(x)=0 and factor out an x from the first n terms. We then 



, , n-1 n-3 n-4 

have: xCa^^x + + a^x + 

f n-1 ^ n-2 ^ n-4 

xCa^^x + + a^x + 



+ a ) + a ^ = 0 

n n+1 



• + a ) = •■ (a -). If the integer truth 
n n+1 



set of f(x) =0 is not empty then there exists an x or (a.x^ ^ + a_x^ ^ 

n-3 . ^ ^ 

+ a_x + • * ‘ + a ) which is an integer factor of - (a ., ). Therefore all 
on ° n+1 

of the integers that belong to the truth set would be factors of - (a ^ ) . 

3 2 » J- 

Example f(x) = x^ - 2x - llx + 12 

q o 

Find {xel | x - 2x - llx + 12 = 0} 

Therefore all members of the truth set would have to be factors of -12. 

Hence the only possible elements of the truth set are +1, +2, +3, +4, +6, +12, -1,-2, -3, 
-4, -6, -12. Using the remainder theorem, factor theorem and synthetic 



i 

















2 



division, we would find that 1,4, and -3 are the only integers that satisfy 
3 2 

X - 2x - llx + 12 = 0. In other words the truth set is {1,4, -3}. 

Checking each one of the possible roots is very time consuming, 
especially if it is done by hand computation. The next problem is how do 
you tell the comprcer to do all the time consuming computations? To keep 
the program simple the writer assumes that the reader can factor the 



constant term in f(x). Furthermore if the term equals zero then zero 

is a root and any other integer root would be a factor of a^. On the 
following page. Fig. 1 is a flow chart telling the computer how to do 
these computations. 

PROGRAM Remarks 



10 READ N 
20 LET 1=1 
30 READ A(I) 

40 LET 1=1+1 

50 IF I<=N+1 THEN 30 

60 READ X 

70 LET F=0 

80 LET 1=1 

90 LET F=F*X+A(I) 

100 LET 1=1+1 

110 IF I <+N+l THEN 90 



N is the degree of F(X), and must be entered in 
Step 170 Data [The value of N]. In steps 20 
through 50 the computer is reading A^, 

• • • TT3 1 oVirMI "\ A 

n * n+1 ' 

in step 170 Data N, A 



These values should be 

A^, • • • A 



^ 2 ’ ^ 3 ’ 
entered 



1’ 2’ “n+1* 
the computer reads your first factor of 
X’s should be entered in 180 Data X , X, 
• • • In Steps 70 through 110, F(Xjf 

F = 0 then X is 
a root. 



In step 60 



n+1' 



All 



X, 



generated. 120 if 
F 5 ^ 0 then X isn’t 



. . "3 

is being 

a root. 120 



if 



120 IF F=0 THEN 

130 PRINT X "IS N0T AN INTEGER R00T 0F F(X)" 

140 G0 T0 60 

150 PRINT X "IS AN INTEGER R00T 0F F(X)" 

160 G0 T0 60 
170 DATA 
180 DATA 

190 END , 4. ^ 

used the computer stops 



170 Data [The value of N]. [The values of A^^, A^ 



X’s]. 



A ]. 180 Data [The values of all the 

^ After all the data in 170 and 180 is 



Using this program, to solve for the integer roots of x - 2x - llx + 12 = 0, 
we get the following: 



FIND I 11:46 PX THU 08/24/67 

1 IS AN INTEGER R00T 0F F[X] 

2 IS N0T AN INTEGER R00T 0F F[X] 

3 IS N0T AN INTEGER R00T 0F F[X] 

4 IS AN INTEGER R00T 0F F[X] 

6 IS N0T AN INTEGER R00T 0F F[X] 

12 IS N0T AN INTEGER R00T 0F F[X] 

-1 IS N0T AN INTEGER R00T 0F F[X] 

-2 IS N0T AN INTEGER R00T 0F F[X] 

-3 IS AN INTEGER R00T 0F F[X] 

-4 IS N0T AN INTEGER R00T 0F F[X] 

-6 IS N0T AN INTEGER R00T 0F F[X] 

-12 IS N0T AN INTEGER R00T 0F F[X] 

0UT 0F DATA IN 60 
FIND I 11:45 PX THU 08/24/67 
10 READ N 
20 LET 1=1 
30 READ A(I) 

40 LET 1=1+1 
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50 IF K=N+1 THEN 30 
60 READ X 
70 LET F=0 
80 LET 1=1 
90 LET F=F*X+A(I) 

100 LET 1=1+1 

no IF I<=N+1 THEN 90 

120 IF F=0 THEN 150 

130 PRINT X’’IS N0T AN INTEGER R00T 0F F(X)" 

140 G0 T0 60 

150 PRINT X"IS AN INTEGER R00T 0F F(X)" 

160 G0 T0 60 

170 DATA 3, 1,-2, -11, 12 

180 DATA 1,2, 3, 4, 6, 12, -1,-2, -3, -4, -6, -12 

190 END 

NOTICE: 170 Data 3, 1,-2, -11, 12. N=3,a^=l,a2=-2,ag=-ll,aj^=12. 

180 Data 1,2, 3, 4, 6, 12, -1,-2, -3, -4, -6, -12. These are the factors of -12, 

Part II Searching 

Part II may be used as a preliminary program lo Parts III, IV, or V to 
locate the interval where a root exists. 

A discussion of the properties and behavior of polynomial functions is 
an essential topic in the study of elementary functions. Quite often a major 
portion of this discussion is concerned with calculating the zeros of a 
real polynomial function P(X), or equivalently, solving the polynomial equation 
P(X)=0 for its roots. Polynomial equations of the first or second degree may 
be solved directly by utilizing, respectively, basic properties of the real 
numbers or the quadratic formula. However, the solution of equations whose 
degree is greater than two is not as straight-forward and requires more 
ingenuity on the part of the student. We will be restricted to a presentation 
of calculating the real zeros of a real, polynomial function whose degree is 
three or greater. We note that odd -degree real equations will always have at 
least one real root since complex roots occur as conjugate pairs. 

A classroom presentation of this subject would be preceded by the 
•following related sequence of topics: remainder theorem, factor theorem, 

fundamental theorem of algebra and synthetic division. With the above know- 
ledge the student is actually capable of solving a higher degree polynomial 
equation. However the solution process is mainly one of trial and error. 

Thus a typical classroom presentation would now involve the selection and 
teaching of time-saving techniques. Although the application of these 
techniques does not necessarily result in an actual solution they are useful 
since they give the number and approximate locations of the real roots. 

A) Descarte’s Rule of Signs: i) The number of positive roots of. a 

polynomial is equal to, or less than, by an even integer, the number of sign 
changes in the coefficients, ii) The number of negative roots of a polynomial 



lERiC 








START 












O 



FIQURE 1 



is equal to, or less than, by an even integer, the number of sign changes in 
the coefficients, if x is replaced by -x* 

B) Newton’s Relations: It may be shown by utilizing Newton’s relations 

that if the roots are real, then the maximum value of the roots, in absolute 
value , is 



C) Location Principle: If any real polynomial and real numbers a and 

b are given such that P(a) and P(b) have opposite signs, then there exists at 
least one real root of the equation P(x)=0 between a and b. 

With the exception of those equations which have relatively small 
integers as their roots, the solution process is still tedious. The assigned 
problems must be few in number and relatively simple in solution. The 
question which often arises is, what course of action is taken if one has a 
large number of these equations to solve or the roots are not small integers? 
The answer is, use a computer. 

The introduction of a computer approach to the solution of polynomial 
equations is particularly relevant at this point since it emphasizes the 

main advantage of a computer it is a timesaving device. In the remainder 

of this paper we will present four programs concerned with the general 
solution of polynomial equations. The first will be a program for searching. 
That is, we will direct the computer to locate approximate values of the 
roots by searching for numbers a and b such that P(a) and P(b) have opposite 
signs. Secondly, we will present three programs which utilize these results 
to actually solve the equaxion for its roots. 

In searching, the general procedure will be to first have the computer 
calculate the upper bound, R for the absolute value of the real roots and 
then to calculate function values at equally-spaced intervals from -R to R. 
The size of the increment, D is arbitrary. The computer will compare 
successive function values for opposite signs, and when they are found will 
print out x^, P(x^), P(x2). A flow chart. Figure 2, and general program 

with explanations appear on the following pages. 

The searching program was run with the following data statement: 

280 DATA 3,20,24,-38,-65,25. 

That is, we were searching for approximate loactions of the roots of P(x) = 

3 2 

24x - 38x - 65x + 25. The following results were obtained: 



X 



max 




^ II n-1 
where P(x) = a^x + a^x + 
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FIGURE 2 



-1.40745 P[X1] = -25.7027 -1.12596 P[X] = 15.7525 

.281489 P[X1] = 4.22754 .562978 P[X] = -19.3551 

2.25191 P[X1] = -40.0038 2.5334 P[X] = 6.67269 

PROGRAM FOR SEARCHING EXPLANATION 

10 READ N,C N = degree, C = number of equal increments 

20 F0R 1=1 T0 N+1 

30 READ A(I) coefficients of a particular pol 3 momial 

40 NEXT I 

50 LET R=SQR((A2)/A(1))*(A(2)/A(1))-2*A(3)/A(1) ) upper bound for roots 

60 PRINT "MAXIMUM VALUE=" R 

70 LET D=2*R/C calculate size of increment 

80 LET X=-R set x equal to lower bound 

90 G0SUB 190 calculate P(x) 

100 LET F1=P set FI equal to P(x) 

110 F0R J=1 T0 C beginning of loop 

120 LET X=X+D increments x 

130 G0SUB 190 calculate P(x+D) 

140 LET F2=P set F2 equal to P(x+D) 

150 IF F1*F2<=0 THEN 240 — if FI and F2 have opposite signs then go to print 
160 LET F1=F2 statement 

170 NEXT J end of loop-returns to line 110 

180 G0 T0 290 have ended search for roots between -R and R 

190 LET P=0 

200 F0R K=1 T0 N+1 lines 190 to 220 evaluate the polynomial at the 

210 LET P=P*X+A(K) last given value of x. 

220 NEXT K 
230 RETURN 

240 LET X1=X-D calculates Xl for output purposes 

250 PRINT X1,"P(X1)="F1,X,"P(X)="F2, 

260 PRINT 

270 G0 T0 160 return to main loop and continue searching 

280 DATA 

290 G0 T0 10 return to beginning if searching for roots of more 

300 END than one polynomial 



We may use results of this type in succeeding programs to actually jcalculaite 
the roots. 

Several facts should be noted at this point. 

i) This program will give approximate locations for the real roots only, 

ii) We can only assert that there is at least one root between the 
printed x values. There may be more. 

iii) If the results of running the program are not completely successful j 
this may be remedied by running the program again and using a larger number 
of equal intervals. For instance, in the above example suppose that we had 
located only two roots. If we suspect that the third root is real we could 
substitute for line 280 the following: 

280 DATA 3,50,24,-38,-65,25 

and run the program again. With this change the computer will calculate 
function values at 50, instead of 20, equally-spaced points in the interval 
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-2.81489 to 2.81489. The greater the number of intervals the more refined- 
the search. 

iv) Another remedy is to increase the value of R. This may be 
necessary if any of the roots are complex since Newton’s relations apply 
only if all the roots are real. 

Part III Finding the Real Roots of any Polynomial, f(x) = A.x’^ + 

A^x +• • *+Ax+AT,By the Half Interval Method 
2 . n n+1 

A. Basic Idea: Two values, x^ and x^ are selected and tried. If it 

is found that f(Xj^) and f(x 2 ) have opposite signs, then a real root exists 
for some value of x between x^ and x^. The new x will equal (x^ + 

see. Fig, 3. This new x will replace either x^^ or x^ and the process will 
continue until the root is found. 

B. Running the Program: After the pre-punched program has been 

entered into the computer, change the three data statements (line numbers 500, 
510, 520) to fit your program. 

Example: 500 DATA 3 (This is the value for the order of the polynomial) 

510 DATA 3, -2, 4, 3 (These are the values of the coefficients: 



^ 1 * ^ 2 * ^ 3 * ^ 4 ^ 

520 DATA -1,0 (These are the guesses for x^ and x^; the 

computer will try to find a root somewhere 
between these values) 

Run the program. 

Note - The computer will determine ONE of the following conditions 



for a given set of data: 
f(x) 




FIGURE 3 
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1. is a root. 

2. X 2 is a root. 

3 . There is no real- root between and x^ . 

4. The exact value of a root. 

5. The value of a root to an accuracy of + 0.000001. 

6. If one of the above conditions is not met af1;er 1000 divisions, the 
computer will print "NO CONVERGENCE". 

C. The program 'written in Basic Language (See Flow Chart, Figure 4) 
PROGRAM COMMENTS 

10 LET J=1 See steps 270,280 

20 READ N Value of the order of the polynomial. 

30 FOR K=1 TO N+1 Steps 30,40,50 get the coefficients into the computer. 

40 READ A(K) 

50 NEXT K 

60 READ XI, X2 Two guesses for x; the computer will determine if 

there is a root between them. 

70 LET F1=0 Steps 70,80,90,100 compute the value of the polynomial 

80 FOR K=1 TO N+1 at x. using the form: f (x. ) = [[[[A.x + A^]x + A_]x 

90 LET F1=F1*X1+A(K) ^ ^ 

100 NEXT K + • • • + A^]x + 



110 IF F1=0 THEN 310 A check to determine if x. is a root. 

120 LET F2=0 Steps 120,130,140,150 compute the value of the 

130 FOR K=1 TO N+1 polynomial at x_ (same form as steps 70-100). 



140 LET F2=F2^*X2+A(K) 



150 NEXT K 

160 IF F1*F2>0 THEN 330-If the product of f(x^) and f(x 2 ) is a positive number 

then no root can be found between them. 

170 IF F2=0 THEN 350 A check to determine if x is a root. 

180 LET X=(Xl+X2)/2 A new value for x halfway'^between x. and x„. 

190 LET F=0 Steps 190,200,210,220 compute the value of the 

200 FOR K=1 TO N+1 polynomial at the new x (same form as steps 70-100). 

210 LET F=F*X+A(K) 

220 NEXT K 



230 IF F*F1<0 THEN 390 — If the product of f(x) and f(x. ) is negative, a root 

lies between them. 

240 IF F*F1=0 THEN 370 — If f(x) is zero then the product will be zero and x 

is a root. 



250 LET X1=X Replace x by x because the product of f(x) and 

is positive. 

260 IF ABS( (X2-X1)/X2)< If interval is within 0.000001, then x« is printed. 

= .000001 THEN 410 

270 LET J=J+1 Steps 270,280,290. After 1000 tries without success 

280 IF J<=1000 THEN 180 in 260, "NO CONVERGENCE" is printed. 

290 PRINT "NO CONVERGENCE" 

300 GO TO 60 (Also 320,340,360,380,400,420) After solving for a 

310 PRINT XI" IS A ROOT" given x. and x the computer goes to step 60 and 
320 GO TO 60 begins again with a new x. and x_ read from line 520. 

330 PRINT "NO ROOT BETWEEN" ^ ^ 

XI "AND" X2 
340 GO TO 60 

350 PRINT X2 "IS A ROOT" 

360 GO TO 60 
370 PRINT X 
380 GO TO 60 
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390 LET X2=X x is replaced by x« 

400 GO TO 260 ^ 

410 PRINT X2 
420 GO TO 60 

500 DATA 3 Data for- step 20. 

510 DATA 3, -2, 4, 3 Data for steps 30,40,50. 

520 DATA “1,0 Data for step 60. 

599 END 




Part IV Method of False Position 

The method of false position utilizes the results of the searching 
program to calculate the roots of a polynomial equation. Figure 5 gives a 
geometrical interpretation of this method, (x^, f^) and (x^, are points 
on the graph of the polynomial function such that f^ and have opposite 
signs. The secant determined by these two points is introduced. By a 
similar triangle argument we may solve for x^, the point of intersection of 
the secant and the x-axis. 



^3 " ^1 
^2 “ ^1 



^ 1-^2 



==3 = 



\ * ^2-^2* 
^ 2-^2 



Note that this equation will also work if f^ < f^. The only stipulation 
is that x^ < x^o In this particular instance since f^ » f 2 < 0 we introduce 
the secant determined by (x^, f^^) and (x^, f^), and precede to solve for x^ 
in the same fashion as before. If > 0 we would have introduced the 

secant determined by (X 2 , f 2 ) and (x^, f^). We see that the points x^,x^,x^, 

“ ° , etc. are converging on the value x which is a root of the polynomial. 





A flow chart. Figure 6 above, and program with explanations follow on the 
next pages. 

This program was run with the following data statements. We are using 
the results of the searching program. 

250 DATA 3,1000,24,-38,-65,25,-1.40745,-25.7027 
260 DATA -1 .12596, 15. 7525,. 281489, 4. 22754,. 562978 
270 DATA -19.3551,2.25191,-40.0038,2.5334,6.67268 
The following results were obtained; -1.25 



.333333 
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3 2 

Substitution into P(x) = 24x - 38x - 65X + 25 reveals that these are 

the exact values of the zeros. 

The following items about the program should be noted: 

i) The maximum number of iterations, Q, is introduced to prevent an 
infinite loop. 

ii) Lines 150 and 190 calculate the difference between two successive 
approximations of the root. 

PROGRAM FOR METHOD OF FALSE POSITION EXPLANATION 

10 READ N,Q N=degree Q=counter 

20 F0R 1=1 TO N+1 

30 READ A(I) Coefficients of the polynomial 

40 NEXT I 

50 READ X1-F1,X2,F2 Values obtained from searching 

60 F0R K=1 T0 ,Q Beginning of loop 

70 LET X3=(X1*F2-F2*X1)/(F2-F1) Calculate X3 

80 LET P=0 Lines 80 to 110 calculate P(x3) 

90 F0R J=1 T0 N+1 
100 LET P=P*X3 + A(J) 

110 NEXT J 

120 LET F3=P Set F3=P(x3) 

130 IF F3=0 THEN 230 IF P(x3)=0 Then print 

140 IF F3*F1<0 THEN 190 If F3 and FI have opposite signs then go to 190 

150 IF ABS(X3-Xl)<=10E-6 THEN 230-Calculate difference between successive 
160 LET X1=X3 approximations 

170 LET F1=F3 Lines 160 and 170 reassign values to XI and FI 

180 G0 T0 220 since FI and F3 have the same sign. 

190 IF ABS(X3-X2)<=10E-6 THEN 230-Same as 150 

200 LET X2=X3 Lines 200 and 210 reassign values to X2 and F2 

210 LET F2=F3 

220 NEXT K End of loop 

230 PRINT X3 

240 G0 T0 50 Go to line 50 for calculation of next root 

250 DATA 
280 END 

If the difference is sufficiently small the process is stopped. 

iii) This method will not work if a given root is of even multiplicity. 

In this case the graph of the function is tangent to the x-axis. 

Part V Finding the Real Roots of Any Polynomial, f(x)=A^x + A^x + 

* * * + A X + A , by Newton’s Method, 
n n+1’ 

A. Basic Idea: A value for x^ is tried near the value of the desired 
root, and f(x^) is computed. The intersection of the tangent to the curve at 
(x^,f(x^)) and the x axis gives the new value of x to be computed (see Fig. 

7 ) . This process is repeated until the root is found or until the desired 
closeness is attained. 

B. Running the Program: After the pre-punched tape has been entered 

into the computer, change the three data statements (line numbers 300,310, 

320) to fit your program. 
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Example: 300 DATA 3 (This is the value for the highest power of x). 

310 DATA 3, “2, 4, 3 (These are the values of the coefficients: 

^ 1 ’ ^ 2 * ^ 3 * ^ 4 ^ 

320 DATA 0 (This is the first trial for x^; it should be 
chosen near the root). 

Run the program. 

Note - The computer will determine ONE of the following conditions 
for a given value of x^: 

1. The value of a root to an accuracy of 0.000001. 

2. If the above accuracy is not attained after 1000 cycles, 
the computer will print NO CONVERGENCE . 



f(x) 
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C. Derivation of the Function: Definition of the slope m of a 

tangent line is: 

m = f(x^)/(x^-X 2 ) 
from which we get: 

^2 = Xi - f(x^)/m 

from calculus, the slope m is the derivative of the function at 
point x^: 

m = f’(x.^) 

so, Newton’s equation for the second trial is: 

X2=Xi-f(Xi)/f »(xi) 

D. The Program Written in Basic Language (See Flow Chart, Fig. 8) 

PROGRAM COMMENTS 



10 LET J=1 See 180,190. 

20 READ N The largest exponent of x. 

30 FOR K-1 TO N+1 Coefficients of the polynomial are entered in 

40 READ A(K) steps 30,40,50. 

50 NEXT K 

60 READ XI A number near the root should be entered. 

70 LET P=0 Steps 70,80,90,100. f(x. ) is computed using the 

80 FOR K=1 TO N+1 form f(x ) = [[[[A x+A_]x + A.]x + • • • + 

90 LET P=P*X1+A(K) ^ ^ 2 3 

100 NEXT K n-'^ \+l^ 

110 LET Q=A(1)“N Steps 110,120,130,140. f’(x^) is computed using 

120 FOR K=2 TO N the form f’(x ) = [[[[A. nx+A-(n-l)]x + A_(n-2)] 

130 LET Q=Q*X1+(N-K+1)*A(K) + . . . + A h2)]x + A^]; w?ere D xn = . 

140 NEXT K n-1 n ’ X 

150 LET X2=X1-P/Q x = x - f(x )/f’(x ); see Part C, Derivation of 

the Function. 



160 IF ABS( (X2-X1)/X2)< If the interval is within + 0.000001, then x_ is 

=.000001 THEN 220 printed. 2 

170 LET X1=X2 Xj^ is replaced by x_. 

180 LET J=J+1 Steps 180,190,200. ^After 1000 tries without 

190 IF J<=1000 THEN 70 success in step 160, "NO CONVERGENCE" is printed. 

200 PRINT "NO CONVERGENCE" 

210 GO TO 60 Steps 210,230. After solving for a given x , 

the computer goes to step 60 and begins again 
220 PRINT X2 with a new x read from line 320. 

230 GO TO 60 

300 DATA 3 The order of the polynomial. 

310 DATA 3, -2, 4, 3 The values of the coefficients. 

320 DATA 0 Your first trial for x . 

399 END 1 




START 
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USE OF COMPUTERS 

IN INSTRUCTION OF BASIC STATISTICS 
by Mildred J. Beavers 
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One mathematics class available to twelfth graders at Boulder 
High School is Senior Math. This class was organized to fill a need for 
additional math beyond one year of algebra and one year of geometry for 
the student who is not primarily interested in further study of mathematics 
or science. 

The classes are composed of several segments. Most who take the 
course are college bound (95%) but will continue in other subject matter 
areas. Some take it to improve college board scores j a few, only because 
they need a credit. Included in the group are many who have had unpleasant 
experiences with mathematics earlier but realize that they should have 
greater knowledge for cultural as well as economic reasons 

Senior Math serves as a review of axioms and basic theorems of 
both algebra and geometry. Methods of solutions for linear and quad- 
ratic equations are rediscovered. New topics include fractional exponents, 
arithmetic and geometric sequences and series , permutations , combinations , 
trigonometric functions, solution of triangles, simple identities, logic, 
and others. All are taught without the rigor that is demanded of math 
majors. Many concepts are treated intuitively; others with proof. 

The principle text book used in the course is Advanced High School 
Mathematics ; Vannatta,CCacnahanI,’a-Fawtfett ; Merrill, 1961. 

This topic chosen for computer adaption is the one on descriptive 
statistics. Formulas are given based on simplified definitions. The 
assumption that students have acquired a beginning knowledge of "Basic 
Language" has been made. Therefore, the first task is to define terms: 
i.e., data, datum, statistics, descriptive statistics as different from 
Inferential statistics. 

Lesson 1 — Arithmetic Mean 



Definition — Arithmetic mean — a measure of central tendency — 
the sum of all data divided by the number of elements in the data, 
Examples 

1) The number which represents the ultimate score which is their 
grade for the quarter or semester. 

2) Average earnings per week based on a sum of 52 weeks* salary. 

Definition — Summation Symbol (This is the first time most of 
these students have seen the summation symbol. ) n 

Ih 

i = 



o 

ERIC 
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Definition — X is used to represent the arithmetic mean . 

n 



X = i J X. 
n 

i = 1 

The exercises on page 254, 1, 2, 3, 4, 5, would be assigned. 

These are artificial situations, overly simplified in nature. 

Example 1. Find the mean of 87.2, 68.5, 74.8, 94.0, 82.2, 96.1. 

The other 4 examples are about the same. 

This is a matter of 10 minutes at the most in the classroom. 

Now, in order to demonstrate the speed and efficiency of the computer, 
I would take the students through the flow chart and the program for the 
computer. These could be written on the board, or the overhead projector 
might be used. Under either set of circumstances, the two would be used 
side-by-side, explaining the flow chart first, and leading into the program. 
Problem 4 of the set is good to demonstrate the program: 

John’s per cent grades on five tests were 84, 72, 91, 64, 83. 

Find the mean. 

After the explanation using Problem 4 and assuming that a computer 
were available, the class would then have a demonstration of their prob- 
lems as solved by a computer. Since they have worked unrealistic examples, 
these examples should now be contrasted with an example of the magnitude 
they might find in the real world. In order to make the contrast, they 
would put through the computer the arithmetic mean of the automobiles per 
state in the United States. The numbers of automobiles are expressed as 1000 

111 DATA 51,1349,72,635,643,8371,878,1250,203,209,2637,1612 ; 

112 DATA 1873,3485,1541,674,1646,289,642,198,275,2644,388,5282 

113 DATA 1699,260,4392,1043,907,4326,361,896,279,1339,4393 i 

114 DATA 405,136,1508,1310,554,1530,148 j 

The arithmetic mean is 1468.69. This indicates 1,468,690 as the average 
number of cars per state in 1965. This took the computer less than 1 
second to do. 



The comparison of the two should begin to show a little of the power j 
of the computer. ; 

Lesson ^ — Other Measures of Central Tendency i 

Definitions i 



1) Array — the raw data is in no special arrangement. If these 
numbers are arranged into an ordered sequence, usually in 
order of increasing size, then this sequence is called an array. 

2) Median — the median is the mid-value of a set of data. If the 
number of data is odd, then the median is the mid-odd value. 
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10 READ N 

20 FOR I = 1 TO N 

30 READ X(I) 

40 NEXT I 



50 LET S = 0 

60 FOR I = 1 TO N 
70 LET S = S + X(I) 



80 NEXT I 



90 LET M = S/N 

100 PRINT "THE 
MEAN=" M 

no DATA 5,84,72, 
91,64,83 

111 DATA 

112 DATA 
115 G0 T0 10 



Program 1 



N is the number of elements in the data. In this 
problem N = 5, because there are 5 test scores that 
we are using. 

This statement directs the X(I) in 30 to take on all 

successive values of X. from i = 1 to 5. These X«*s 

1 1 

are those values listed in the data. 

The computer will read 84 the first time, since X is 
84; the second time through X^ ” ^2’ "^2. 

This statement says that the computer should use the 
X the first tim^ it reads X(I); X2 j the second; X^, 

the third, and so on until the list of data is used. 
In this problem the computer would go through this 
loop 5 times. 

Since the sum of the data is the first calculation 
we are after, we start with a sum of zero and add 
each datum to the previously calculated sum. 

This statement tells the computer to go use the X’s 
in proper sequence for the summing operation. 

This is key statement for the addition process. We 
started with a sum of zero in 50. 60 said read first 

X^. 70 says add 84 (XI) to 0. This is the statement 

that does the summation operation, giving us partial 
sums as the computer uses successive X^’s. 

This statement tells the computer to go back and 
read 72 and then in 70 add 72 to 84. The next time 
add 91 to the sum of 72 + 84, namely 156. New S in 
70 is now 247. This is a partial sum. The computer 
keeps doing this until the last X^ , in this case X^ 

(83), is added to the previous partial sum. 

This command tells the computer to divide the sum 
by N(5). 

Here the statement "The mean is 78.8" 



Other problems 
Other problems 

If you wish to do several problems the computer will 
go back to 10, read new N, and new data. If you 
have several data lines for one stop program change 
115 to "G0 T0 20." 



120 END 
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Thus, if there are 9 numbers in the data, the 5th number 
would be the median. If there are an even number of data, 
then the median is the arithmetic mean of the middle two 
numbers. 12 numbers in an array would yield a median which 
would be the arithmetic mean of the 6th and 7th numbers. 

3) Range — the range is the difference between the largest and 
smallest numbers of the array. 

The examples in the text are again trivial in nature. A quick 
calculation of Problems 1 and 2 will give enough practice to understand 
the principles involved. Problem 4 — "Will the median be affected if 
the numbers at the extremes of an array are changed? Will the mean be 
affected by such changes? Can the mean and median ever coincide?” This 
problem lends itself to experimentation with the data. First the sugges- 
tion could be made to the class that they try adding data to the array in 
patterns; i.e., adding numbers at either end which were greater and lesser 
by the same value, adding a value to the upper end that is twice that of 
the one added to the lower end of the array. All sorts of arrangements 
of this nature can be explored and intelligent guesses recorded in a 
temporary place. 

Then on the overhead projector or chalkboard, the flow chart and 
program could be written. In this lesson, instead of giving a complete 
explanation immediately, give the class a chance to recognize what they 
can in both. The recognition will be limited, but commands which are 
mathematical in nature as well as previously used statements can be noted. 
Then go through both flow chart and program, giving a detailed account 
of what is happening. Example 1 in the set of exercises will serve to 
illustrate. 

Example 1 — Find the median of $64, $82, $51, $90, $67, $71, $58, 

$94, $63. What is the range? (See Figure 2.) 

/< 

Again, the problem involving the number of automobiles would 
serve to illustrate the speed and competency of the computer. The time 
that the problem took the computer to form an array and calculate range 
and median should be pointed out so that students are reminded of one 
of the chief assets of the computer. To emphasize this, the instructor 
might even time the students as they tried to do the same problem. 

Let them pick any part of the problem to do but use the same set of 
data. Start with a stop watch, and call time. Using the same data 
as in the previous lesson, the computer forms an array, calculates 
the range (8299) and the median (1043) in 6 seconds. Compare the 
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Program for Forming an Array > Range , 



and Median 



10 READ N 
20 FOR 1=1 T0 N 
30 READ A(I) 

40 NEXT I 

50 LET K=1 
60 LET J=K+1 

70 IF A(K)<=A(J)THEN 110 



80 LET T=A(K) 

90 LET A(K)=A(J) 
100 LET A(J)=T 



110 LET J=J+1 
120 IF J<=N THEN 70 

130 LET K=K+1 

140 IF K<N THEN 60 
150 LET D=A(N)-A(1) 
160 FOR 1=1 TO N 



170 PRINT A(I) 

180 NEXT I 
190 PRINT 

200 PRINT "RANGE="D 

210 IF INT(N/2)=N/2 THEN 
240 

220 LET M=A((N+l)/2) 

230 GO TO 250 



N is number of elements in data. 

20-40 — this subscripted variable technique 
was explained in previous lesson. It simply 
instructs the computer to read the data in 
sequence . 

50-60 — these two steps set up a routine, 
making J one unit more than K where J and K 
are both values of I. 

Compares A(K) to A(J). If A(K) is less than 
or equal to A(J), then the numbers are in an 
order with the smaller one first, which is 
what the program is telling the computer to do — 
arranged in ascending order. If A(K)^A(J) then 
the computer goes to 110. 

80-100 — if A(K)>A(J) this means they are not 
in ascending order, and the order must be 
changed. To do this, the commands as they 
are written assign the value A(K) to a new T. 

Then A(K) takes on value A(J), whatever this 
value was. A(J) then gets the value located 
in T which was the first A(K) . These commands 
reverse the order of two pieces of data. 

Says for the computer to increase the value 
of J by one preparing the computer to continue 
the comparison of data. 

If J is^N then the computer returns to 70 to 
compare the new A(J) with A(K). If A(K)^ 

A(J), then the program continues through steps 
110,120 until either J>N or A(J)<A(K). If A(J)< 
A(K) steps 80-100 are repeated. 

Sends the computer back to start the process 
all over again, starting with the second datum 
and comparing it to the rest of the list of 
data in the same way that it did the first 
item of data. 

The computer continues as long as K is less 
than N. When K=N then the computer goes 
automatically to 150. 

At this point the data is an array. The computer 
finds the smallest datum and subtracts it from 
the largest. 

Statements 160-180 command the computer to punch 
out the array. These statements could precede 
statement 150. 



Spaces the print-out sheet so that the information 
is orderly and easy to read. 

Prints the range identiJ^ing it as such with. §< 
quotation saying this. 

Here, the notation INT(N/2) must be explained. 

If these are equal then we are ready to find 
Median 1. If not, continue to statement 220. 

Gives the computer the number of the last term 
of the first half of the terms of the data. 

Sends the computer to statement 250 which will 
print out the value of this particular A(I). 



I 

I 



1 

•J 

i 



i 

] 

s 

I 

i 

















23 




Flow Chart S-M-R 

Sorting into an Array, Range, Median 
Figure 2 











students’ accomplishments with that cf the computerc 
Lesson 3 — The Mode 

The mode is such a simple concept that five minutes of discussion 
is usually enough » 

Definitions 

1) Mode — The mode of a set of numbers is that number which appears 
more frequently than any other in the sere 

2) Bi~Modal — If there are two modes in a set of data, the data is 
said to be bi-modalo 



The exercises in the text are adequate to illustrate the ideao 
Exercise 3 is one that will serve as a basis for discussion: On a very 

easy arithmetic test given to a large number of students, would the mode 
of the scores likely be higher or lower than the mean? On a very difficult 
test, would the mode of the scores likely be higher or lower than the mean? 

Included is a program which will calculate the mode as well as 



form an array and find the range and mean^ These means are not included in 



these lessons, because of the nature of the class; but it may be of interest 



to others using a programming approach to the teaching of statistics « 



10 READ N 
20 F0R 1=1 TO N 
30 READ A(I) 

40 NEXT I 
50 F0R I=1T0 N 
60 LET Q(I)=A(I) 

70 NEXT I 
80 G0 SUB 450 
90 F0R 1=1 T0 N 
100 LET A(I)=Q(I) 

110 NEXT I 

120 LET D=A(N)-A(1) 

130 F0R 1=1 T0 N 
140 PRINT A(I): 

142 NEXT I 
144 PRINT 

146 PRINT "RANGE="D 

150 IF INT(N/2)=N/2 THEN 180 

160 LET M=A(N/2) 

170 T0 T0 190 

180 LET M=(A((Ntl)/2)+A((N-l)/2))/2 

190 PRINT "MEDIAN="M 

200 FOR 1=1 T0 N 

210 LET S(I)=0 

220 NEXT I 

230 LET 1=1 

240 F0R K=2 T0 N 

250 IF A(I)=A(K) THEN 290 

260 LET I=K 

270 NEXT K 



280 G0 T0 310 
290 LET S(I)=S(I)+1 
300 G0 T0 270 
310 F0R 1=1 T0 N 
320 LET Q(I)=S(I) 

330 NEXT I 

340 G0 SUB 450 

350 IF Q(N)=1 THEN 420 

360 F0R 1=1 T0 N 

370 IF Q(N)=S(I) THEN 400 

380 NEXT I 

390 G0 T0 iO 

400 PRINT A(I) "IS A M0DE" 

410 G0 T0 380 

420 G0 T0 10 

430 PRINT "N0 M0DE" 

440 G0 T0 10 

450 LET K=1 

460 LET J=K+1 

470 IF Q(K)<=Q(J) THEN 510 

480 LET T=Q(K) 

490 LET Q(K)=Q(J) 

500 LET Q(J)=T 
510 LET J=J+1 
520 LET K=K+1 
530 RETURN 

560 END 












Lesson ^ — Other Means — Geometric , Quadratic , Harmonic 

This section is treated as an optional section in teaching, depending 

upon the background and interest of the students in a particular class. 

Included are programs which will calculate these three means if they are 

of importance to the reader. 

Geometric Mean 

10 READ N 
20 F0R 1=1 T0 N 
30 READ X(I) 

40 NEXT I 
50 LET S=1 
60 F0R 1=1 T0 N 
70 LET S=S*X(I) 

80 NEXT I 
90 LET G+S (1/N) 

100 PRINT "GE0METRIC MEAN="G 
no DATA 5,243,81,27,9,3 
120 DATA 5, 3, 2, 6, 8, 5 
130 G0 T0 10 
140 END 

Quadrati c Mean 

10 READ N 
20 F0R 1=1 T0 N 
30 READ X(I) 

40 NEXT I 
50 LET S=0 
60 F0R 1=1 T0 N 
70 LET S=S+X(I)*X(I) 

80 NEXT I 

90 LET Q=SQR(S/N) 

100 PRINT "QUAD. MEAN="Q 

no DATA 9,6,6,7,8,8,8,11,11,12 

120 DATA 8, 2. 5, 2. 5, 2. 5, 3. 5, 3. 5, 3. 8, 4, 4. 2 

130 DATA 5, 3, 2, 6, 8, 5 

140 G0 TO 10 

150 END 

Harmonic Mean 

10 READ N 
20 F0R 1=1 T0 N 
30 READ X(I) 

40 NEXT I 
50 LET S=0 
60 F0R 1=1 T0 N 
70 LET S=S+1/X(I) 

80 NEXT I 
90 LET H=N/S 

100 PRINT "HARM0NIC MEAN ="H 
120 DATA 5, 3, 2, 6, 8, 5 
130 G0 T0 10 



Lesson ^ — Mean Deviation — Semi-Interquartile Range 
Definitions 



1) Mean Deviation — a measure of variability. The mean deviation 
is the arithmetic mean of the absolute values of the deviations 
from the mean of a set of data. 

The formula for finding the mean deviation is 

n 



M.D. - - i 
n ^ 



X. - X 
1 



1 = 1 

2) Semi-interquartile Range — equal to one-half the range between 
the number, which is the last value in the first quartile and 
the last number in the third quartile. A quartile is one part 
of an array when the array is divided into four equal parts. 
This measure of variability is the average of the quartile 
deviations. The formula for finding the semi-interquartile 
range is 



Q = 



<53- 



The semi- interquartile range can be calculated by hand after 
using the computer to form an array, or it could be calculated 
by a program designed to follow the earlier lesson involving 
the array. 

Since the problems in the textbook are simplified situations, they 
would be used to illustrate the principle involved. Hand calculating 
these sets of data will lead into an explanation of how the computer pro- 
gram works. 



Computer Program for Mean Deviation 

10 DIM X(100),D(100) 

20 READ N 
30 F0R 1=1 T0 N 
40 READ X(I) 

50 NEXT I 
60 LET S=0 
70 F0R 1=1 T0 N 
80 LET S=S + X(I) 

90 NEXT I 
100 LET M=S/N 
110 LET T=0 
120 F0R 1=1 T0 N 
130 LET D(I)=X(I)-M 
140 LET T=T + D(I) 



150 NEXT I 

160 LET T1 = T/N 

170 PRINT"NUMBER OF VALUES"N 

180 PRINT 

190 PRINT 

200 PRINT "S X VALUES", "DEVIATIONS" 

210 F0R 1=1 T0 N 
220 PRINT X(I),D(I) 

230 NEXT I 
240 PRINT 
250 PRINT 

250 PRINT "MEAN = "M,"MEAN DEVIATI0N="T1 
270 DATA 
300 END 



5 



-1.12197 XIO 



If the problem involving automobiles is used again, the mean deviation 
-7 

, Computer time to calculate is one second. 






Lesson £ — Measure of Variability — Standard Deviation 
Definition 

Standard Deviation — The standard deviation is a measure of the 
average amounts by which individual items of data vary or deviate 
from the arithmetic mean of all the numbers that compose the data. 
Standard deviation is the quadratic mean of the individual devia- 
tions from the arithmetic mean. (If the class has done the section 
on quadratic means, then this definition is simple enough. If the 
class has not, then stop to define a quadratic mean at this time.) 
The formula for the standard deviation is; 



0 





The problems in the text are simple, as indeed they must be. The 
hand calculations are still tedious. If the computer is available, then 
one problem from the text would be enough to understand the principle 
involved. Exercise 1: Compute the standard deviation of .6,. 6,. 7,. 8, 1.0, 

2,1.4. 



After this simple problem has been finished by the class, ask for 
suggestions as to how the previous computer program for mean deviation 
could be adapted to find standard deviation. The program for mean deviation 
should be shown on the overhead or written on the board with enough spaces 
between lines so that student suggestions could be inserted and tried to 
see if they will work. Hopefully, the class would discover a step or 
steps that would have to be inserted to handle the squaring process and 
finding the square root. 



Program for finding the standard deviation that has been developed is 



as follows (this is the one we did in 

10 DIM X(I),D(I) 

20 READ N 
30 F0R I = 1 T0 N 
40 READ X(I) 

50 NEXT I 
60 LET S = 0 
70 F0R I = 1 T0 N 
80 LET S = S + X(I) 

90 NEXT I 
lOOLET M = S/N 
110 LET T = 0 
120 F0R I = 1 T0 N 

Again, the example using the 



class ) : 

130 LET D(I) = X(I)-M 
140 LET T=T + D(I)*D(I) 

150 NEXT I 
160 LET T1 = T/N 
170 LET T2 = SQR(Tl) 

180 PRINT N 
190 FOR I = 1 TO N 
200 PRINT X(I),D(I) 

210 NEXT I 
220 PRINT M,T2 
230 DATA 
250 END 

number of cars in the United States 



or some similar problem will illustrate the capabilities of the computer. 
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The use of standard deviation in the interpretation of data 5 the 
normal probability curve , percentages of the population falling within a 
certain number of standard deviations > etc., is discussed later in the 
chapter. This particular unit of work takes up frequency distributions, 
class intervals, ccmiputation of these previously discussed measures from 
the frequency distribution, binomial expansion and normal distribution, 
sampling techniques, standard error. All of this is handled in a most 
general way. 

In summary, these lessons have been designed in a very limited 
manner to show the capacities of a computer. At all times during these 
lessons one should stress that there is nothing magic about what happens 
in the machine; it can only do what a human programmer tells it to do. 
Since these students will probably never be programmers, details of how 
the program works have been explained; but no requirements as to knowledge 
of programming should be required. 

However, the occasional student could become interested in learning 
programming, and the instructor should be prepared to help this student 
progress in understanding and technique as far as possible. 
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USING THE HARDY-WEINBERG LAW ON A COMPUTER 
by Robert L. Brown 
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The purpose of this paper is to show how the use of a computer 
could make the Hardy-Weinberg Principle of Biology more meaningful to a 
high school biology student® This could be accomplished in two ways: 

1) by making what can be very difficult and tedious numerical computations 
easier (especially with a large batch of data) and 2) by showing some 
insight into the mathematics of the Hardy-Weinberg law itself. To under- 
stand the terms used, this material should not be introduced to students 
until the subject of genetics has been studied. 

The Hardy-Weinberg law was first proposed independently in 1908 
by the English mathematician, G® H® Hardy, and German physician, W. 

Weinberg. The law states that, given two alleles (A and a) which occur 
in the frequencies of p and q, after one generation of random mating the 
proportions of the three possible genotypes (AA, Aa and aa) will be in 
equilibrium and will stay in equilibrium generation to generation thereafter® 
The Hardy-Weinberg law explains one of the most striking and fundamental 
facts about nature: that, while bi-parental populations always contain 

a variety of genotypes, the population as a whole, including its variaxions, 
may continue for generation after generation without significant change. 
Variation, which makes evolutionary change possible, is maintained even 
when evolutionary change is not occurring. 

For the population to stay in equilibrium, one must make some 
assumptions. They are: 1) the population is large, 2) that mating is 
random, 3) that the alleles involved have no selection value and 4) that 
mutations do not occur, and if they do the mutation rate is in equilibrium. 
Not only does the Hardy-Weinberg law enable one to study genetic traits 
that have no selection advantage in a population, but it enables the 
biologist to study alleles for traits that do have selection advantages 
and to compute the rate at which the frequencies of the alleles change. 

To get the law in numerical form we let p = the frequency of the 
*’A" allele (the dominant allele if there is one) and q = the frequency 
of the "a" allele (the recessive allele). The choice of gametes pro- 
duced by a sexually reproducing organism can be represented as p + q = 1. 

The probability of the genotypes produced by two parents can be represented 

2 2 

by the product of their gametes or (p + q) = 1 , which yields the Hardy- 

2 2 2 

Weinberg law: p + 2pq + q =1 where p represents the AA genotypes, 

2 

2pq represents the ^ genotypes and q represents the ^ genotypes in 
the next generation. When there is dominance involved the AA and Aa 
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genotypes will have the same phenotype (will look the same)# but the 
aa homozygous recessives can be picked out of the population. This. 

means that by finding the aa’s in the population you have found the 

2 o 

q *s# and by taking the square root of q you have found q. Since 

p + q = 1, then p = 1 - q; and you now have the frequencies of p and 

q(A and a alleles) with which you can now compute the ratios of p^:2pq:q^ 

(AA:Aa:aa) in the population. 

As an example, take the trait found in humans called mid-digital 

hair (hair on the middle digit of the fingers). Genetic studies have 

indicated that the allele for mid-digii:al hair (M) is dominant to 

the allele for no mid-digital hair (m). A person with the trait would 

have the genotype MM or Mm, and a person without the trait would be mm. 

We will assume that mid-digital hair has no selection value and that 

matings involving this trait are random. 

Suppose in a sample of a given population we find that 36% of the 

people do not have mid-digital hair. We then know that mm = q^ = .36 

and that q = .6 and then p = .4. It follows that MM = p^ = (.4)^ = 

.16 = 16% of the population, and Mm = 2pq = 2x.4x.6 = .48 = 48%; and 

2 

we already know that mm q = 36% (.16 + .48 + .36 = 1 or 16% + 48% + 

36% = 100%). 

F 3 ?om the above we can compute the odds for the offspring from dif- 
ferent types of marriages. For example, what would be the chances of 
two parents, both with mid-digital hair, producing a child without 
mid-digital hair? For this to happen, both parents would have to be of 
the genotype Mm or heterozygous for this trait; and the chances of their 
both being this way would be the answer. For the father, the chances 
are Mm/(MM + Mm) = 48/(16 + 48) = 3/4. The chances for the mother are 
the same. The chances for their both being Mm is 3/4 x 3/4 or 9/16. 

The odds can also be computed for other types of marriages. See Chart 1 

at the end of this paper for a canputer program which, when given q^ will 

2 2 

then compute p, q, p“, 2pq and q . The program is set up such that 
2 2 

one q or many q values can be fed in. 

This can further be expanded to predict the frequency of marriages 
or matings in the population involving a given pair of alleles. Using 
mid-digital hair as our example, this would involve all the possible 
marriages of people with the genotypes MM, Mm and mm. Using p = .4 
and q = .6, the frequencies of all the possible marriages would be as 
follows (in using this a geneticist would probably be interested in 
a certain type of mating rather than all of the possible matings): 



MM X MM marriages = x ~ (c4)^ 

(since genotypes are the same, no need to specify sex) 

MM X Mm 



= .0256 



2 3 3 

marriages - (p x 2pq)2 = 4p q = 4(.4) x .6 = .1536 



(since genotypes are different, we must allow 
Mm X MMJ for parents being of different sex; hence, 

(p^ X 2pq) X £ above) 

2 2 2 2 

Mm X Mm marriages = 2pq x 2pq 4p q = 4(o4) (o6) 



MM X mm 
mm X MM 



'} 

:} 



marriages = (p^ x q^)2 = 2p^q^ •= 2(.4)^(r6)^ 



Mm X mm\ . 2._ 3 

mm X MmJ == (2pq x q )2 = 4pq 

4 4 

mm X mm marriages = q = ( . 6 ) 



4(.4)(c6)' 



total 



.2304 

.1152 

.3456 

.1296 

1.0000 



When all the decimal frequencies total one, you know that you have 
not made a mistake nor left out one of the possible matings. From 
the above we can see that the frequency of marriages between mates 
with mid-digital hair and those without mid-digital hair is oll52 
plus .3456 or 46% of the marriages (the frequencies for p and q used 
in this example were picked for convenience and do not represent the 
national frequencies for this pair of alleles). 

2 

In reality the frequencies are not perfect squares (q ), and 

computations involving the Hardy-Weinberg law are easily adapted to 

computers . Charts two and three at the end of this paper show two 

computer programs that will determine the frequency of various possible 

matings when given the percentage of aa or homozygous recessives in 

2 

the population (q ). 

For traits that are sex linked, like red-green color blindness 

and hemophilia, we must modify our approach in finding the frequencies 

of p and q. A color blind woman (X^^) would be represented by q^, but 

a color blind man (X^Y) has only one allele for this trait and would be 

represented by q. Since 0.6% of the females in the United States are 
2 

color blind (q = 0.006), then we would expect 8% of the males (square 

root of 0.006 = 0.08) to be color blind. Population samples have shown 

this to be true, and understanding the Hardy-Weinberg law explains why 

the percentage of color blind or hemophiliac males is so much greater 

than the same traits in females — hence, the term sex linked traits. 

From the above data we can compute the odds that any woman picked at 

A di 

random is a carrier for the color blind allele (XX). Since q = 0.08, 
then p = 0.92; and the carrier women would be the 2pq’s or 2 x .92 x .08 
= .1472 or about 15%. The odds are 15/100 that any woman picked ax 
random is a carrier for color blindness. 
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The Hardy-Weinberg law is also very useful in working with popu- 
lations with alleles where one allele does have selection value over 
the other. In this case, the population would not be in equilibrium for 
this trait; in fact, it would be evidence that natural selection is 
favoring certain genotypes. 

Next we will look at selection against dominant and recessive 
genotypes. To simplify matters, we will assume that the alleles are 
fully penetrant and that the genotypes are lethal or such that the 
organism never reproduces (dies before reproductive maturity, or if it 
reaches maturity is unable to reproduce and thus unable to pass the 
"bad" allele on to the next generation). 

In the case of the bad dominant allele that appears through 
mutation (remember if either of the parents would have shown the trait, 
they would not have reproduced), the allele would be eliminated in one 
generation; since the organism would never reproduce and thus never pass 
the bad allele on to the next generation. The genotypes are almost 
always heterozygous or Aa for the trait, since the A allele would arise 
by mutation. Thus, in one generation the bad trait is eliminated except 
for new mutations from the recessive to the dominant allele. In reality 
for this case the frequency of the bad allele will be nearly equal to 
the mutation rate of a ^ A. In any case, the allele is quickly eliminated 
from the population gene pool. 

Selection against a bad recessive allele is much less effective. 

The lower the frequency of the bad allele, the less the change in 
frequency of the genotype from generation to generation. Remember, we 
are assuming that these aa genotypes do not reproduce. Most of the aa 
bad genotypes are produced by Aa heterozygous parents who are normal 
for the trait. The exceptions would have to involve mutations. When 
we assume that certain types of mental deficiencies are inherited, we 
find that most of the mentally defective children come from normal 
parents — not from mentally defective parents. The sterilization of all 
mentally defective people at birth would not lead to a significant 
reduction in mentally defective offspring the next generation. 

If with full penetrance the frequency of aa is q^, and if no 
affected person reproduces , the number of aa individuals in the next 
generation can be calculated from the frequency q^^ of the £ alleles 
in the reproducing population. The total frequency of the reproducing 
population is the sums of AA and Aa = p + 2p q , while the a alleles 
in this population is 1/2 of the frequency of the Aa genotypes, or PqQq* 

O 

ERIC 
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The ratio of a among all alleles after disregarding the non-reproducing 
genotypes is: 






P Q 
^ 0^0 



p +2p q p +2q 
^ ^0 0 o ^o 



where q^ = initial frequency. 



q^^ = 1st generation, = 2nd generation, etc. 

Since p = 1 - q this becomes : 

^o o 






3 

1 

1 



I 



^1 1 + q 



(substitute p^ above) 



and the frequency of aa genotypes in the next generation is: 

, ,> 1 -^ " 

‘(^1 + 

If we start with an initial frequency of 1% (.01) for the aa (a^)* 

_ 2 / /7oT ' 

then q, =| ■ ■ - 

\1 + /Toi/ ‘ \1 + /.0083 

Thus, in our first generation we only have a reduction of .17%; and 
moving to the second generation we get a reduction of only .14%. The 
reduction becomes less and less as the frequency of q becomes less. 
See Graph 1 at the end of the paper. 

In general terms, after two generations: 



(Dl) 



^1 " 1 + q 



(D2) 



"^2 " 1 + q , 



‘o "1 

by substituting the right side of (Dl) for q^ in (D2) above, we get: 





1 + q. 



^2 = 



1 + 2q 



1 + 



1 + q. 



or a general formula to find the reduction for several generations: 

.2 






2 . 






where (q^) is equal to the frequency 



of aa genotypes expected from the nth generation. See Chart 4 for a 
computer program that carries this out for 10 generations. The graph 
in Graph 1 was plotted from this program starting with an initial 
frequency for q^ of 1%. After 10 generations, the frequency of aa 
genotypes has been reduced from 1% to only .25%. Graphs 2 and 3 show 
the result of 10 generations starting with initial frequencies of .1% 



o 
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and .01% respectively for aa genotypes. 

This result shows that the reduction in frequency of aa genotypes 
becomes less and less in each successive generation and that after 10 
consecutive generations of complete selection, it is still 1 /m- of the 
initial frequency of 1%. To reduce it to 1/10 (.1%) would require 
22 generations, from .1% to c01% would require 68 more generations, 
from .01% to e001% would require 216 more generations, and from .001% 
to .0001% would require 684 more generations = When one attempts to com- 
pute the above percentages one soon realizes the advantages of computers, 
which will carry out the numerical operations for whatever generation 
desired in less than a second. 

This has been an attempt to show how computers could aid in the 
teaching of high school biology. There are many ideas in biology which 
could be aided by computer techniques. Some areas involved today are 
information storage, growth forms and patterns, growth rates, taxonomic 
problems and many types of biochemical and physiological problems, to 
mention only a few. 

A good reference for applications of the Hardy-Weinberg law to 
human populations is: STERN, Curt, 1960. Principles of Human Genetics , 

733 pp. , Wo H. Freeman and Co. , San Francisco. 


















Chart 1 



Program 



10 


PRINT "Q 


M npn 

9^9 


20 


READ Q2 




30 


LET Q = 


SQR(Q2) 


40 


LET P = 


1 - Q 


50 


LET P2 = 


P*P 


60 


LET P3 = 


2*P*Q 


70 


PRINT Q, 


P, P2, 


80 


G0 T0 20 




90 


DATA .00005, .04 


100 END 





’2PQ", ”Q2’ 



Q2 



.09, .16, .25, .36, .49, .64 



81 



Answers 

HARDY 13:41 PX FRI 06/30/67 







AA 


Aa 


aa 


1 


P 


P2 


2PQ 


Q2 








.0140421 




7.07107 E-3 


.992929 


.985908 


1.40421 E-2 


.00005 


.2 


.8 


.64 


.32 


.04 


.3 


.7 


.49 


.42 


.09 


.4 


. 6 


.36 


.48 


.16 


.5 


.5 


.25 


.5 


.25 


.5 


.4 


.16 


.48 


CD 

CO 

• 


.7 


.3 


.09 


.42 


.49 


.8 


.2 


.04 


.32 


.64 


.9 


.1 


.01 


.18 


.81 



0UT 0F DATA IN 20 

A computer program written in "BASIC" computer language to find 

2 2 

the values of q, p, p , and 2pq when given q using the Hardy-Weinberg 

2 

law. To enter new values for q you would change the data in line 90, 
the DATA statement. Note that the value of q has to reach 0.7 for 50% 
of the population to show the recessive trait. For more information on 
setting up a program in "BASIC" language, see the 1967 edition of the 
General Electric "BASIC" Language reference manual put out for their 
computer time-sharing service. 
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Chart 2 



Program 

10 READ Q2 

20 LET Q = SQR(Q2) 

30 LET P 1 - Q 
40 LET PI = Ps’fpApAp 
50 LET P2 = 4*P*P*P*Q 
60 LET P3 = 4*P*P*Q*q 
70 LET P4 = 2*P*P*Q*Q 
80 LET P5 = 4*P*Q*QAQ 
90 LET P6= QAQ*Q*Q 



100 PRINT Q2,Q, P, 
110 DATA .3, .04 
120 G0 T0 10 
130 END 


PI, P2, 


P3, P4, P5, P6 




Answers 




Population 30% aa 


Population 4% aa 


(Q2) 




0.3 (30%) 


0.04 (4%) 


(Q) q 




0.547723 


0.2 


(P) p 




0.452277 


0.8 


(PI) AAxAA marriages 




0o0418427 


0.4096 


(P2) AAxAa marriages 




0.202691 


0.4096 


(P3) AaxAa marriages 




0.245466 


0.1024 


(P4) AAxaa marriages 




0.122733 


0.0512 


(P5) Aaxaa marriages 




0.297267 


0.0256 


(P6) aaxaa marriages 




0.0900 


0.0016 



A computer program to find the frequencies of all the possible 
matings involving alleles which have no selection value. The ( ) 
represent the symbols for the above statements in the computer program. 
For other values of q enter new data for line 110, the DATA statement 
(convert all percentages into their decimal equivalents). The program 
is set up in such a way that one or a hundred bits of data could be en- 
tered at one time. 



37 




Chart 3 



Program 



o 
1 — 1 


READ Q2 




20 


LET QP SQR(Q2) 




30 


I.ET P = 1 - Q 




40 


LET P1 = (P*P*P*P) + 


(4*P&P*P*Q) + (4*p*p*Q*Q) 


50 


LET P2 = (2*P*P*Q*Q) 


+ (4*P*Q*Q*Q) 


60 


LET P3 - Q*q*q*q 




70 


PRINT Q2, Q, P, PI, 


P2, P3 


80 


DATA .23, .04, .09, 


.16, 125, .36, .49, .64, .81 


90 


G0 T0 10 




IOC 


1 END 





Answers 
2 

q =.23 
P3=.0529 

.04 
.0016 

.09 
.0081 

.16 
.0256 

.25 
.0625 

.36 
.1296 

.49 
.2401 

.64 
.4096 

.81 
.6561 

0UT 0F DATA IN 10 



A computer program to find the frequencies of all the possible 
matings of three types: 1) both parents show the dominant trait, 2) one 

parent shows the dominant trait, and one parent the recessive trait; and 
3) both parents show the recessive trait. PI in the program is that 
part of the expected matings where both parents show the dominant trait; 
P2 is one parent dominant and one parent recessive; P3 is both parents 
recessive. The answers repeat for a given q in the order indicated 
(q , q, p, PI, P2 , and P3). 



q=. 479583 


p=. 520417 


Pl=.5929 


P2=.3542 


.2 


.8 


.9216 


.0768 


.3 


.7 


.8281 


.1638 


.4 


.5 


.7056 


.2688 


.5 


.5 


.5625 


.3750 


.6 


.4 


.4096 


.4608 


.7 


.3 


.2601 


.4998 


.8 


.2 


.1296 


.4608 


.9 


.1 


.0361 


.3078 



4 

I 

I 



i 



] 

I 



I 



i 

I 



i 
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Chart 4 



Program 






5 PRINT "GENERATI0N", "FREQUENCY 0F GEN0TYPE IN PERCENT" 
10 READ Q2 

20 LET Q(0) = SQR(Q2) 

30 FOR N = 1 T0 10 

40 LET Q(N) = (Q(0)/(1 + N*Q(0) ) ) t2*100 
50 PRINT N, Q(N) 

60 NEXT N 

70 DATA .01, .001, .0001 
80 G0 T0 10 
90 END 



Answers 



GENERATI0N 


FREQUENCY 0F 


1 


.826446 


2 


.694444 


3 


.591716 


4 


.510204 


5 


.444444 


6 


.390625 


7 


.346021 


8 


.308642 


9 


.277008 


10 


.25 


1 


9.39633 E-2 


2 


8.84571 E-2 


3 


8.34212 E-2 


4 


7.88033 E-2 


5 


7.45586 E-2 


6 


7.06477 E-2 


7 


6.70368 E-2 


8 


6.36957 E-2 


9 


6.05983 E-2 


10 


5.77215 E-2 


1 


9.80296 E-3 


2 


9.61169 E-3 


3 


9.42596 E-3 


4 


9.24556 E-3 


5 


9.07029 E-3 


6 


8.89996 E-3 


7 


8.73439 E-3 


8 


8.57339 E-3 


9 


8.41680 E-3 


10 


8.26446 E-3 


0UT 0F DAT^. IN 


10 



GEN0TYPE IN PERCENT 



A computer program to calculate 
complete selection against a recessive 



genotype. These answers were used to 

2 



plot Graphs 1, 2, and 3o Each q" 
listed in the data (% of aa genotypes 
changed to decimal equivalents) was run 
for 10 generations. To increase or de- 
crease the number of generations, you 
would change the 10 in line 30 to the 
number of generations desired. The E-2 



found in the answers means the i' numerical... 

-2 



value X 10 , as the computer puts the 

answers in this form to give the benefit 
of more significant digits. In line 40 
the equation , \2 






was multiplied by 100 to put the answers 
in % form to make graphing easier. The 
computer used one second of time to com- 
pute the 30 answers given. These were 
stored in a memory unit and then read out 
by the teletype after all values had been 
computed. To get these answers by hand 
would require several hours. 



/ 



ERIC 
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Chart 5 



90 DATA 


.0001 


,.0002, .0003,. 


0004, .0005, .0006, 


.0007, .0008, .0009 




91 DATA 


.001, 


.002, .003, .004 


,.005, .006, .007,. 


008, .009 






92 DATA 


.01,. 


02, .03, .04, .05 


,.06, .07, .08, .09 








93 DATA 


.1,.2 


,.3, .4, .5, .6,. 


7, .8, .9,1 








RUN 














HARDY 


14 


:41 PX M0N 


07/03/67 








Q 




P 


P2 


2PQ 




Q2 


.01 




.99 


.9801 


.0198 




.0001 


1.41421 


E-2 


.985'858 


.971916 


2.78843 


E-2 


.0002 


1.73205 


E-2 


.982679 


.965659 


.034041 




.0003 


.02 




.98 


.9604 


.0392 




.0004 


2.23607 


E-2 


.977639 


.955779 


4.37214 


E-2 


.0005 


2.44949 


E-2 


.975505 


.95161 


4.77898 


E-2 


.0006 


2.64575 


E-2 


.973542 


.947785 


.051515 




.0007 


2.82843 


E-2 


.971716 


.944231 


5.49685 


E-2 


.0008 


.03 




.97 


.9409 


.0582 




.0009 


3.16228 


E-2 


.968377 


.937754 


6.12456 


E-2 


.001 


4.47214 


E-2 


.955279 


.912557 


8.54427 


E-2 


.002 


5.47723 


E-2 


.945228 


.893455 


.103545 




.003 


6.32456 


E-2 


.936754 


.877509 


.118491 




.004 


7.07107 


E-2 


.929289 


.863579 


.131421 




.005 


7.74597 


E-2 


.92254 


.851081 


.142919 




.006 


.083666 




.916334 


.839668 


.153332 




.007 


8.94427 


E-2 


.910557 


.829115 


.162885 




.008 


9.48683 


E-2 


.905132 


.819263 


.171737 




.009 


.1 




.9 


.81 


.18 




.01 


.141421 




.858579 


.737157 


.242843 




.02 


.173205 




.826795 


.68359 


.28641 




.03 


.2 




.8 


.64 


.32 




.04 


.223607 




.776393 


.602786 


.347214 




.05 


.244949 




.755051 


.570102 


.369898 




.06 


.264575 




.735425 


.54085 


.38915 




.07 


.282843 




.717157 


.514315 


.405685 




.08 


.3 




.7 


.49 


.42 




.09 


.316228 




.683772 


.467544 


.432456 




.1 


.447214 




.552786 


.305573 


.494427 




.2 


.547723 




.452277 


.204555 


.495445 




.3 


.632456 




.367544 


.135089 


.464911 




.4 


.707107 




.292893 


8.57864 E-2 


.414214 




.5 


.774597 




.225403 


5.08067 E-2 


.349193 




.6 


.83666 




.16334 


2.66799 E-2 


.27332 




.7 


.894427 




.105573 


1.11456 E-2 


.188854 




.8 


.948683 




5.13167 E 


-2 2.63340 E-3 


9.73666 


E-2 


.9 


1 ■ 




0 


0 


0 




1 



0UT 0F DATA IN 20 
TIME : 3 SECS . 

This is the same program as listed in Chart 1 with more data. Only 

the new data and answers are listed. This chart contains or approximates 
2 

many of the q values used in Hardy- Weinberg problems. Remember that the 
2 

q percentages from the population must be converted into their respective 
decimal equivalents. 
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Graph JL 



Percent 




Complete selection against a recessive genotype e Initial frequency 
of genotype in the population was 1%. After 10 generations, the fre- 
quency has been reduced from 1% to .25%. After 10 generations, the 
frequency is still 1/4 the initial frequency. 






Graph £ 



II 



lo 

0 



0 






Percent 




Complete selection against a recessive genotype c Initial frequency 
of the genotype in the population was 0.1%. After 10 generations, the 
frequency has been reduced from 0.1% to 0.058%. After 10 generations, 
the frequency is still over 1/2 the initial frequency. 
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Graph 3 




Generation 

Complete selection against a recessive genotype. Initial frequency 
of the genotype in the population was 0.01%. After 10 generations, the 
frequency has been reduced from 0.01% to 0.0083%. After 10 generations, 
the frequency is still over 4/5 the initial frequency. (The lines do 
not necessarily become closer as the % drops due to rounding off. 

See Chart 4 for the actual data before rounding off. ) 



PYTHAGOREAN TRIADS 

by Dewayne N. Burgdorff and Jay E» Niebur 

Introduction 

The main purpose of this lesson is to show the relationship that 
exists between the numbers which satisfy the Pythagorean Theorem., These 
numbers are called Pythagorean triples or triads. All triads or triples 
are integers such that the sum of the squares of two will always equal 
the third. In other words, if the integers «re a, b, and c, then by 
applying the Pythagorean Theorem a^ + b^ should equal c^. 

Historical Aspects of Theorem 

The Pythagorean Theorem is one of the most important theorems 
in the realm of geometry. Not only does it apply to plane or solid 
geometry, but it has application in analytic geometry, as well as in 
the field of trigonometry. The origin of this famous theorem dates 
back before Pythagoras generalized it. 

The Egyptians used the relationship of the 3, 4, 5 right triangle 
to lay out their land after it was flooded by the Nile River each year. 
They used ropes measured off with these dimensions and from this acquired 
the names of "rope stretchers." Also, about the same time the Egyptians 
were using this relationship, the Hindus in India had developed more 
advanced ideas concerning right triangles. Besides using the 3,4,5, their 
right triangles included such measurable values as the 12,16,20; 15,20,25; 
^»i2,13; 15,36,39; 8,15,17; and the 12,35,37. Neither the Egyptians 
nor the Hindus tried to prove why this relationship held true. 

It was about 540 B.C. before Pythagoras used the Egyptian rule of 
2 2 2 

thumb, 3 +4 = 5 , to generalize this now famous theorem. It has not 

been proven that Pythagoras truly was the one who developed this theorem 
or whether some other Greek found the answer. Be that as it may, this 
mathematical rule will go down in history as the Pythagorean Theorem. 

Since the Pythagorean Theorem was formulated, many individuals 
have tried to prove this theorem. For example, in 1940 Elisha S. Loomis, 
a mathematics professor at Baldwin Wallace University, wrote a book 
titled. The Pythagorean Proposition . It contained over 370 formalized 
proofs . 

Probably the most famous of all these proofs was a geometrical 
proof discovered by Euclid. It has been stated that there is no better 
logically derived proof than that of Euclid. Following is this. famous 
proof and three other simple algebraic proofs. 









THE PYTHAGOREAN THEOREM 
(EUCLID’S PROOF) 
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Given ; Triangle ABC, Angle ACB is 
on Sides AB, BC, and AC. 
Prove ; Square ABHK = Square ACDE + 



1. 


Draw CM // AK 


1. 


2. 


Draw CK and BE 


2. 


3. 


In triangles KAC and BAE, 
AK = AB; and AC = AE. 


3. 


4. 


/_ KAC = /^ EAB 


4. 


5. 


Triangles KAC and ABE are 
congruent . 


5. 


6. 


DCB is a straight line. 


6 . 


7. 


Triangle BAE and rectangle 


7. 




ACDE have same base and alt 
BAE = rectangle ACDE 


• 


8. 


8. 


9. 


KAC and rectangle KMLA have 
same bases and altitude. 

^ rectangle KC^LA 


9. 


“1 n 

XV • 


10. 


11. 


X rectangle KMLA = x squa:oe 
ACDE. 


11. 


12. 


Rectangle KMLA = Square 
ACDE. 


12. 



i right angle, squares constructed 
Square BFGC 

A line can be drawn // to a given 
line. 

Two points determine a line. 

Adjacent sides of a square are 
equal . 

If equals are added to equals 
their sums are equal. 

S.A.S. = S.A.S. 

If two angles are supplementary, 
their exterior sides lie in a 
straight line. 

Parallel lines are everywhere 
equidistant . 

If a triangle and a parallelogram 
have same bases and altitudes, the 
area of the triangle equals one-half 
that of the parallelogram. 

Same as Number 7. 

Same as Numbers 5 and 8. 

Substitution axiom. 

If equals are multiplied by equals 
their products are equal. 
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13. 


Draw AF and CH 


13 


14. 


In triangles ABF and HBC, 
BC = BF; and AB = BH. 


14 


15. 


/ ABF = / HBC. 


15 


16. 


A ABF 2 A HBC. 


16 


17. 


ACG is a straight line. 


17 



18. 


A ABF and rect. BFGC have 
same bases and altitude. 

A ABF = i rect. BFGC. 


o 

00 
1 — 1 


19. 


19. 


20. 


A HBC and rect. MHBL have 
same bases and altitude. 

^ HBC = ~ rect. MHBL. 

~ sq. BFGC = ~ rect. MHBL. 


20. 


21. 


21o 


22. 


22. 


23. 


Sq. BFGC = rect. MHBL. 


23o 


24. 


KMLA + MHBL = ACDE + BFGC. 


24, 


25. 


KMLA + MHBL = sq. KHBA. 


25. 


26. 


Sq. KHBA = Sq. ACDE + BFGC. 


CO 

CM 



Another means by which this theorem 



Two points determine a line. 
Adjacent sides of a square are 
equal o 

If equals are added to equals, 
their sums are equal. 

So AoSo = S i A= So 

If two angles are supplementary, 
their exterior sides lie in a 
straight line. 

Parallel lines are everywhere 
equidistant . 

Same as Number 8. 

Same as Number 18. 

Same as Numbers 5 and 8 inclusive. 
Substitution axiom. 

Multiplication axiom. 

Addition axiom. 

The whole is equal to the sum of 
its parts. 

Substitution axiom. 

may be expressed is: 
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Method Given four identical right triangles placed in such a location 
to form a given square ABCD (as shown in Figure 3). 



Proof: 




The square ABCD has sides c, and its area is c^. The area of 
the square is made up of the four triangles plus the smaller 
inside square. The area of the smaller square is equal to 
(a “ b)^. The area of the four triangles is 4(^ ab). The 
area of the large square, c^, equals the sum of the smaller 
square, (a - b)^ and the four triangles, 4(^ab). Algebraically: 



which simplifies to 

c^ 

Finally 



(a “ b)^ + 4(-^ ab) 

(a^ - 2ab + b^) + 2ab 
a^ + b^ 




Proof: Use the same four triangles used in Method I. Then, on the 

basis of the area, you have the given relationship: 



(a + b)^ = 4(-^ ab) + c^ which simplifies to 
a^ + 2ab + b^ = 2ab + c^ Finally 
a^ + b^ = c^ 



E "erjc 
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Method III , Given right triangle ABC with altitude CD drawn to the 
hypotenuse of the triangle (as shown in Figure 5), The 
two smaller triangles are similar. 



Proof: 




Since triangles ACB, ADC, and CDB are similar, their corresponding 
sides are proportional. From these proportionalities, the 
following proportions are formed: 



1) x/a = a/c 
but it is also true that: 



2) y/b = b/c 



3) X + y = c. 

From the equations 1) and 2), it can be found that 
X = a^/c and y = b^/c. 

By substituting the expressions for x and y into equation 3), 
You obtain: a^/c + b^/c = c 



Then by clearing, fractions and multiplying both sides by c, 
you will have the Pythagorean Theorem, 

a^ + b^ = c^ 



Pythagorean Triads 

Pythagorean triads are three integers which represent the lengths 
of three sides of a right triangle. In this discussion the letters a 
and b are the legs of any right triangle, and the letter c will corres- 
pond to the hypotenuse of the same right triangle. 

There are two main classifications of Pythagorean triads. They 
are classified as multiple values or unique. The multiple valued niunbers 
are obtained from unique triads which have been multiplied by a constant 
k such that ka^ + kb^ = kc^. 

In discussing triads which are unique, there exist certain rela- 
tionships which hold true. They are as follows: 

1, In unique triads, no two numbers may have a common factor. 

They are relatively prime numbers. 
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2. If a, b, and c are triads, then a and b cannot both be odd. 

3. If a, b, and c are triads, c must always be odd; and of 
the other two integers, one must be even and the other odd. 

4. When the sides of a right triangle are integers, the perimeter 
of the triangle will be an even number; and its area will also 
be even. 

Since the development of the Pythagorean Theorem, mathematicians 
through the ages have tried to develop a general formula for the formation 
of these triads. Some of the men worthy of mention in addition to 
Pythagoras (582-501 B.C.), are Plato (429-348 B.C.), Euclid (300 B.C.), 
Maseres (1721-1824), Leonard E. Dickson, and Dr. Artemus Martin. 

Following is a list of the rules formulated by the preceding men. 
Although there is some similarity in all of the rules, each of these men 
are given credit for finding formulas which will give Pythagorean triads. 

The rules of finding integral values of a, b, and c are stated in 
order of their discovery.* 

!• Rule of Pythagoras : Let n be any odd number, then n, (n^ - l)/2 

2 0 0 
and (n + l)/2 are three such numbers. For n^ + ((n - l)/2)^ = 

((n^ + l)/2)^. 

2» Platons Rule : Let m be any even number divisible by 4; then 

2/2 9 

m, m /4 - 1 and m /4 + 1 are three such mimbers. For m’^ + 

((m^/4 - 1)) = ((m^/4 + 1))^. 

3. Euclid *s Rule : Let x and y be any two even or odd numbers 

such that X and y contain no common factor greater than 2 
and xy is a square. Then ✓xy , (x - y)/2 and (x + y)/2 are 
three such numbers. The rule may then be expressed as: 

(v/^)^ + ((x - y)/2)^ = ((x + y)/2)^. 

4. Rule of Maseres : Let m and n be any two even or odd numbers, 

• 2 2 9 

m IS greater than n and (m + n )/2n is an integer. Then m , 

2 2 2 2 

(m - n )/2n and (m + n )/2n are three such numbers. For 
+ (m^ - n^)/2n = ((m^ + n^)/2n)^. 

5. Dickson’s Rule : Let m and n be any two prime integers, one even 
and the other odd, m greater than n and 2mn a square. Then 

m + >^2mn , n + /2mn , and m + n + /2mn are three such numbers. 

For (m + /2mn) + (n + >/2mn ) = (m + n + /2mn ) . 



*Elisha S. Loomis, The Pythagorean Proposition , (Michigan; Edward 
Brothers, Inc., 1940), pp. 19-20. 
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2 2 

Martinis Rule; Let (x+y) =(x-y) + 4xy. 



Assume 4xy 



2 2 

is a square. If you assume x = mp and y = mq and we have 
2 2 2 

4xy = 4m p q , which is a square number for all values of 

2 2 2 

m, p, and q, then by substitution we arrive at (p - q ) + 

2 2 2 2 

(2pq) = (p + q ) . 

A flow chart is a diagram of the steps followed by a computer in the 
solution of a problem. 

In the flow chart for the problem called TRIADS (see Figure 6) the 
odd integer 3 is placed in the computer as N and is used to calculate 
A, B, and C from the formulas A = N^, B = ( ^ and C = ( "" . 

If C does not equal A + B, then the values calculated for N = 3 would be 
discarded; and the computer would take the next odd integer (N = 5) and 
perform the same calculations. If C does equal A + B, as it does when 
N = 3 (C=25, B=16,A= 9 so 25 = 16+9), then the computer will 
output the triad /K, }/b, andv^; or in this case 3, 4, and 5, which is 
a right triangle. This is a unique triad, and all multiples of the 
3, 4, 5 combination will be right triangles also. Therefore, beginning 
with a multiplier (I) of 2, we take I x I x /b and I x y/c, or 
2 X 3, 2x 4, and 2x5, which is another right triangle with sides 6, 

8, and 10. 

Letting I = (previous)! x 1, a new I is generated so that I now 
will equal 3. This multiplier 3 is taken times the unique triad, pro- 
ducing a 9, 12, 15 right triangle. I will be incremented again and more 
multiples found until A1 becomes greater than 100. In that case, the 
computer goes to the next step and decides if N < = 100. N is still 3, 
which is an odd number less than 100; so the program goes to connector A 
and selects a new N. This new N will be 5, and the computer will perform 
all the previous steps but using N = 5. After N = 99 has been calculated, 
the program will follow connector C and will receive an even number. 

N = 4 will be the first situation, and values for A, B, and C will be 
calculated using formulas A = N^, B = ( ^^ - — ^)^, and C = ( ^^ ^ 

The program will proceed as before until N > 100. 

In the flow chart called GUESS (see Figure 7), an integer A, 
then an integer B will be placed in the computer. Calculating X = A^ and 
C = /a^ + B^, the computer decides if - (C - 1)^ > X. If the answer 
is yes, a new value for A will be taken and new calculations begun. 

If the answer is no, the computer will decide if C is less than or 
greater than the integer C. If the answer is yes, meaning C is not an 



o 
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Integer, then a new value for B will be selected, and different calcula- 
tions made using the new B. If the answer is no, meaning C is an integer, 
then the values for A, B, and C v/ill be printed; and these three numbers 
will be triads. 

Through the ages, many hours have been spent in rigorous methods of 
solutions for the calculations of triads. Now, by means of a computer, 
one is able to obtain, in a relatively short amount of time, many combi- 
nations of triads in relation to the Pythagorean Theorem. The data was 
obtained by using the GE 265 computer. By constructing a flow chart, the 
following problems were programmed. 

Problem 1 — By applying the rules of Pythagoras and Plato, try to obtain 
all possible combinations of the Pythagorean triads. 

In the programming of this problem, place the shortest leg of the 
triangle so that it will range from 3 to 100. By using the computer, 
the results were obtained in only 12 seconds . It can be seen thar this 
rule does not solve for all values because of the limits set on a. It 
was most interesting to observe that rules dating from 500 B.C. and a 
modem computer can obtain the results in such concise form. 

Problem 2 — By developing a new rule, try to obtain the possible com- 
binations of triads . 

A searching technique program was constructed for the computer. 

Certain values for a and b were fed into the computer, a ranged from 3-50 
and b from 4 to 1000. In the process of running the program, the computer 
would print all triads in combination with a, b, and c. The results were 
obtained for values of a from 3 through 50. An interesting aspect of this 
program is revealed in the determination of the triads . The amount of 
computer time taken for this program was 75 seconds. 

It can be concluded, through experimentation with the preceding rules, 
that no one rule can apply in all instances. In comparing the two preceding 
problems, it is quite evident that the older, established rules give 
better, more concise results. The computer enables its user to save many 
man hours of work in arriving at solutions of new, as well as old, problems. 
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TRIADS 




Problem 1 



LISTNH 
5 PRINT 

10 F0R N=3 T0 99 STEP 2 
20 LET A=N+2 
30 LET B=((N+2-l)/2)t2 
40 LET C-((N+2+l)/2)+2 
50 IF C=A+B THEN 75 
70 G0 T0 80 

75 PRINT SQR(A),SQR(B),SQR(C) 

76 G0SUB 220 
80 NEXT N 

110 F0R N=8 T0 100 STEP 4 
120 LET A=N+2 
130 LET B=((N+2-4)/4)+2 
140 LET C=((Ni2+4)/4)+2 
150 IF C=A+B THEN 175 
170 G0 T0 180 

175 PRINT SQR(A),SQR(B),SQR(C) 

176 G0SUB 220 
180 NEXT N 
220 LET 1-2 

230 LET A1-I*SQR(A) 

240 LET B1-I*SQR(B) 

250 LET C1-I*SQR(C) 

260 PRINT "MV"A1,B1,C1 

270 LET 1=1+1 

280 IF AK=100 THEN 230 

290 RETURN 

350 END 

RUN 



FLOW CHART 
GUESS 



START 



3 



Yes 



— 

B 








X = A^ 






j 


C = /A^ + B'^j 

1 


2 - (C - 


■ 1)2 > 



- 



X 



TRIADS 



11:07 



PX WED 07/12/67 



No 



A 


B 


C 


3 


4 


5 


MV 6 


8 


10 


MV 9 


12 


15 


MV 12 


16 


20 


MV 15 


20 


25 


MV 18 


24 


30 


MV 192 


• • 

4606 


• • 

4610 


100 


2499 


2501 


MV 200 


4998 


5002 


MV 200 


4998 


5002 



RETURN BEF0RE G0SUB IN 290 



TIME: 12 SECS. 



CoINT(C) 



Yes 





No 


A, B, C 







STOP 



BYE 



Figure 7 



0FF AT 11:30 PX WED 07/12/67. 
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Problem 2 



CU-CE-T-SHAR 

EXPERIMENTAL SYSTEM L0ADED UNTIL 9:30. PLEASE LIST "EXPERT***", 

0N AT 8:27 PX FRI 07/14/67 TTY 3 

USER NUMBER— P61000 
SYSTEM —BAS 
NEW 0R 0LD— 0LD 
0LD PR0BLEM NAI'IE— GUESS 



READY. 

FINKI 

LIST 



GUESS 8:26 PX FRI 07/14/67 

5 PRINT "A"; "B"'; "C" 

10 FOR A = 3 T0 25 
020 F0R B-4 T0 1000 
25 LET X = A*A 
30 LET C = SQR(X + B*B) 

40 IF C*C -(C-’1)*(C-1)>X THEN 80 
50 IF C<>INT(C) THEN 70 
60 PRINT A;B;C 
70 NEXT B 
80 NEXT A 
90 END 



10 F0R A = 3 TO 100 
10 F0R A = 3 TO 50 
RUN 



GUESS 8:29 PX FRI 07/14/67 

ABC 

3 4 5 

5 12 13 

6 8 10 

7 24 25 

8 6 10 

8 15 17 

• • • • • 

28 195 197 

29 420 421 

•• ••• ••• 

48 575 577 

49 168 175 

50 120 130 

50 624 626 



TIME: 1 MINS. 15 SECS. 



RAN 450/6 SEC. 

ST0P. 

READY. 

UNSAVE 

READY 

BYE 

*** 0FF AX 8:36 9 MIN. 



CHARGE. 



O 




NUMERICAL ANALYSIS AND DEMONSTRATIONS FOR 
P.S.S.C. LABORATORY EXPERIMENTS 
by Wayne B. Daniels 

Numerical Analysis of P.S»S«C» Laboratory Experiment 
I-l Short Time Intervals 
References 

Physics Laboratory Guide , 2nd Edition; Physical Science Study 
Committee, D. C. Heath and Company, 1965, pp. 103. 

Physics , 2nd Edition, Physical Science Study Committee, D. C. Heath 
and Company, 1965, pp. 12-14 "Repetitive Motions; the Stroboscope." 

An objective of this laboratory experiment is to calibrate the vibrat- 
ing clapper of an electric bell by determining its period of vibration. The 
vibrating clapper and a tape moving through the clapper will be used to study 
the motions of objects used in later laboratory experiments. A second 
objective of this experiment is to gain some appreciation for the disc- 
stroboscope as a simple device useful in measuring rapid periodic motions. 

The numerical analysis part of the experiment requires that the student ac- 
tually analyze the relationship between the vibrations of the clapper and 
the rotations of the disc-stroboscope and express this relation in an 
algorithm. The learning experience gathered from the logical derivation of 
this algorithm and from the building of the accompanying flow chart is 
the primary objective of this revised experiment. 

The first periodic motion measured is that of a weighted vibrating 
steel blade clamped to the edge of the table. The motion is slow enough 
to allow the observer to "stop" the motion of the end of the blade by 
rotating the disc-stroboscope with only one of its twelve slits open. 

The rate of vibration of the weighted blade is increased by adjusting the 
weight closer to the secured end of the blade. 

The bell clapper is powered by a 1.5 volt battery and is attached 
to a desk top. To reduce the rate of vibration, the experimenter might 
weight the end of the clapper with a clothespin for the first stroboscope 
reading. By a trial and error process, the experimenter must rotate the 
stroboscope at different speeds and open the slits one by one until the motion 
of the clapper has apparently stopped. The student should detect a rela- 
tion between the frequency of the clapper, the frequency cf the stroboscope, 
the number of slits opened in the stroboscope, and the apparent stopping 
of the motion of the clapper by the rotation of the stroboscope. 

The algorithm and flow chart shown in Figure 1 illustrate how the 
student’s trial and error technique could be replaced by a computer’s 
search for the number of slits and the frequency of the stroboscope 



needed to "stop" -^he motion of the clapper. 

Let = the frequency of the clapper 

f = the frequency of the stroboscope needed to stop the 
motion of the clapper 
n = the number of slits in the stroboscope 
k = any whole integer greater than zero 

Then, (k)(F )=(n)(F^) 

O o 

KF 

F = — ^ 
s n 

Values F , the fr .quency of the clapper , and values F , the 
^ s 

frequency of the stroboscope, must enter the computer as input. Since 
the frequency of the stroboscope depends upon hand turning, its range of 
frequency should extend from thirty to ninety rotations per minute. The 
maximum number of slits available on the stroboscope is twelve. This 
flow chart lets the number of slits begin at six and allows the mamber 
to range down to one and up to twelve depending upon the value of k which 
is established in the algorithm. 

Consideration must be taken that if the frequency of the clapper is 
any whole integer multiple (k) greater than the frequency of the one-slit, 
disc-stroboscope, the stroboscope will appear to stop the motion of the 
vibrating clapper. In this case, however, a unique solution for the 
frequency of the clapper cannot be determined. If the frequency of the 
clapper is any unit fraction (i.e., ~, ~, -I*, etc.) of the frequency 

of the stroboscope, then the stroboscope will appear to stop the motion of 
the clapper at several locations of the rotating slit. Again, no unique 
solution for the frequency of the clapper is possible. To eliminate the 
possibility of such multiple frequencies, the value of (k) is tested by 
the computer. Should the value of (k) be less than one or greater than 
two, fewer or additional slits are selected ranging from one to twelve, 
forcing (k) within the required limits. 

The cOTiputer will supply the experimenter with the frequency of 
the clapper being tested, the lowest frequency of the stroboscope needed 
to "stop" the motion of the clapper, and the number of slits needed. 

Summary Questions 

1. What adjustments may be made to the disc-stroboscope to allow 

the experimenter to measure frequencies greater than that of the stroboscope? 

2 . The flow chart shows an input of clapper frequencies and 
stroboscope frequencies. How does the output differ from the input? 
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3. if “he rcnr.purer should ssle:*: more slits than are available on 
the strcOcscope , what warning v/cuid the computer gi /e the experimenter? 

4, Under what ocnditicns would the computer select more slits 
than a'e available on the strobe scope? 

5.. V7hat provision has been made in the flew chart. Figure 1, to 
make sure tnat the cemputer does not match stroboscope frequencies with 
multiples and unit fr-actiins of clapper frequencies? 

Answers to Summary Questions 

1. The number of slits opened in the stroboscope must be increased ^ 

2. The outpur from the ccm.puter will list the low^est frequency 

and the least number of sli“s needed to match the frequency of the clapper ^ 

3. The cemputer will print cut "upper limit," and then proceed 
to the next clapper frequency. 

4^ The computer- would select mere slits than are available on the 
stroboscope if the clapper v/ere vibrating with a frequency greater than 
twelve times the maximum frequency of the strcbosccpec 

5c Two logic boxes located in the flow chart test (k) for a value 
greater than two or less than onec If (k) is greater than two or equal 
to two, the number of slits is increased by one. If (k) is less than one, 
the number of slits is reduced by one,. When 1 £ k > 2, F is multiplied 
by (k) to find the lowest stroboscope frequency and the least number of 
slits needed to "stop" the motion of the clapperc 
Numerical Demonstrations for P . S o S o C . Laboratory Experiments 

II-ll Waves from Two Point Sources 



11-12 interference and Phase 



11-13 Young ' s Experiment 
References 



Physi.es Laborarory Guide , 2nd Edition; Physical Science Study 
Committee; D. C, Heath and Company, 1965, pp. 28-31 

Physics , 2nd Edition; Physic.al Science Sxudy Committee; Dc C« 

Heath and Company, 1965, Chapter 17, "Interference," pp« 282-292 c 

The main purposes of the three experiments are to provide a situation 

where the experimenter can observe the behavior of waves generated at two 

point sources in a ripple tank and apply the resulting conclusions to 

light waves passing through two slits in a coated glass slide » Waves 

generated at two poinr sources will interfere with one another, producing 

a pattern of nodal lines whose changes are the primary concern of this 

experiments The angle tha^; a given nodal line makes with the line 

connecting the two point sources is a function of the wavelength of the 

wave, the distance between the point sources, and the phase delay of the 
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waves at the second source. The purpose of the numerical demonstration 
of this experiment is to show how the angle of a nodal line changes as 
values of wavelength, distance between point sources, and phase delay change. 

Figures 2 and 3 illustrate the mathematical derivation of the 
algorithm, showing the dependence of the angle of the ndfal line. 

~ ^ point located on the nodal line (n). 

X = the distance of the point fran the bisector. 

6 = the angle formed by the bisector and a line connecting 

the midpoint between the point sources and point P . 

n 

d = the separation between the point sources. 

A = the wave length of the waves at both point sources . 

p = the phase delay or that fractional pare of the wave’s 
period by which it is delayed. 

pA = the distance by which one wave is delayed compared to 
the other. 

n = the number of the nodal line right or left of the bisector. 

Fran Figure 3 

PSiPn - PS^P^ = (sin 6)(d) = i 

d 

(n - ~ + p)A 
sin 6 = 

Q 

The flow chart. Figure 4, shows the manner by which the angle 0 shall 
change. Nodal line (two) is selected at random for demonstration purposes. 

As the distance (d) between the point sources is decreased from four centi- 
meters to one centimeter by steps of one centimeter, the phase delay is 
lengthened from zero to one by steps of 0.2, and the wave length is increased 
from one to three centimeters by steps of one centimeter. Sin 6 is printed 
corresponding to each combination of d, p, and A. 

Summary Questions 

1. Which variable in the algorithm is instructed by the flow 
chart to change first? 

2. How many values of sin 0 would be printed by the instructions 
of this flow chart? 

3. Rewrite the flow chart in a way in which it would demonstrate 
whether or not there is equal angular spacing between nodal lines when 
the variables are held constant. 

Answers to Summary Questions 

1. That variable read- in by the inner loop of the flow chart 
changes first. It is the variable (d). 

2 . Seventy-two . 








1 





Figure 1 
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(PSj) < 4 > (PSj^) 



Figure 3 
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Figure 4 
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SOME USES OF THE COMPUTER IN PHYSICS TEACHING 
by William E. Einert and Hepburn Ingham 

Since computers can do prodigious amounts of arithmetic with great 
speed and accuracy, they have important application in solving some physics 
problems . In physics the attempt is made to understand nature in a quanti- 
tative way. A physical variable is related to another variable or a constant 
in a precise way . For a value of x one tries to determine by suitable arith- 
metic the value of f(x). One can work with a known functional relationship 
to determine f(x), or can work with data to determine what the functional 
relationship is. 

Many topics taught in the customary high school physics course could 
be augmented, enriched, or more easily understood by computer techniques. 

The purpose of this paper is to identify a few of these topics and to show 
how they might be introduced and taught through the use of the computer and 
allied techniques. 



Of equal importance is the gain which can be made by the student in 
learning about computer methods and in establishing foundational understanding 
on the use of computers. In other words, here the computer is actually used 
as an adjunct to problem solving and an aid to understanding. It is not 
taken in isolation (as it would be in a separate course) or as an interesting, 
novel "gimmick” (as it might be in a mathematics class). 

The following topics are presented mostly as hints to the knowledgeable 
teacher. They are not meant as finished teaching units, but could serve as 
curriculum guide lines for the teacher in presenting and utilizing each con- 
cept . Most of the material could be given to the student in printed form 
as it is written here if it is accompanied by a well planned lecture , class 
discussion, laboratory experience, and appropriate exercises. The material 
is written with the assumption that the student has had an introduction to 
BASIC language and some experience with flow charting. 

An interesting set of physical relationships is met early in 
the P.S.S.C. physics. These exist on the subject of scaling. Intui- 
tively most students believe that, for any physical structure, the full- 
sized object will have all of the properties of a similar miniature 
model. That is, they are of the same material; and the large structure 
is simply increased by some constant in every dimension over the small 
model. This belief turns out to be surprisingly false. If the dimen- 
sions of the model increase over a given range, the relationship of 
support strength to weight varies widely for an increase in dimension (R). 

This is shown by the graph in Figure 1. The support strength is directly 



proportional to the cross-sectional area of the supports for the structure, 
which is related to (R^). On the other hand the weight is directly 
related to the volume, which is related to (R^)o One can see that nhe 
weight increases much more rapidly than suppom strength for an increase 
in dimension (R). 

An illimiinating exercise would be to increment (R) in unit 
steps on the computer and to get the strength-to-weight ratio. This 
ratio would decrease to some critical value (Q) , an allowed safety 
factor, and would give some idea of maximum size of the structure. 

The problem could be extended and (R) allowed to increase to where the 
strength-to-weight ratio equals unity, at which the structure would 
fail. The flov? chart for the program could be as shovm in Figure 2. 

The parameters in the program are defined as follows: 

R = dimension of the structure 

“ 3 constant relating dimension (R) to crossectional area 
or support strength 

1^2 “ ^ constant, relating dimension (R) to volume and weight 
of the structure 

Q = a safety factor, a minimum safe ratio of strength to weight 

A problem that the student might like to consider is that of a 
cylindrical support column whose material has x Ibs/in^ compress ional 
strength. The value of (x) for various materials can be found in a hand- 
book along with the density (p) of the material. The relationships 
would be: 

Strength « Area = (R in)^ (x Ibs/in^) 

and = (tt)(x Ibs/in^) 

Weight a Volume = (TfR^h in^Kp Ibs/in^) 

and = Trp Ibs/in^ 

A general area in which computers do valuable work is in the 
handling of experimental data. If one is looking for the function that 
exists between two variables it is desirable to have some aid to do arith- 
metic. This is especially true where a large number of data values are 
gathered in a short time. A specific example occurs in the refraction 
of light. Here the student can rapidly measure many angles of incidence 
and refraction (0^) and (6^) respectively. In searching for the functional 
relation that exists between these variables, the student needs to do 
such operations as (0^/0^), (0^ x 0^), (sin 0./sin 0^), etc. If a 
computer could be utilized the student could be encouraged to take 
many more data measurements in limited labor’atory time. 
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Another problem in physics is getting the slopes of curves at 
several points. For example, taking the slope of a velocity-time 
graph to get the rate of change of velocity. By definition this is the 
acceleration. One could reproduce an acceleration graph from a velocity- 
time graph as illustrated in Figure 3. 

The computer could be programmed to compute the slopes at a large 
number of points. By so-doing, the student could get a good feeling for 
the significance of slope and could make a much more accurate represen- 
tation of the acceleration-time graph. 

An important mathematical function in physics is the inverse square 
relationship. This is encountered when studying the effects of radia- 
tion, such as light, on a surface at some distance, (x), from a point 
source. The function also applies to electrical and gravitational force 
fields, which act from a static charge center or a gravitational center. 

In general, the relationship can be expressed as y « l/x^. In the case 
of radiation, the intensity Cl) on a surface at a distance (R) from a 
point source is given by: I = constant /R^. The functions for the 

electric and gravitational force fields are defined by F = KQ^Q.^/R^ 
and F = GM M /R^. in a given problem the terms KQ^Q^ and GM.M^ can 
be constant, and a computer used to find the force (F) for anj^ separation 
of the centers (R). A flow chart for a gravitational problem could be 
as in Figure 4. 

In a satellite problem (see Figure 5) the masses of the satellite, 

earth, and the universal gravitational constant combine in the equation 

to be For a measured or given set of radii, R^^, R^, . . . R^, 

the corresponding forces F , F^, . . . F can be calculated. In this 

1 z n 

case the gravitational force provides the necessary centripetal deflection. 
Knowing the centripetal force, satellite mass, and radius one can also 
calculate the velocity and time of revolution of the satellite. A consis- 
tent set of units must be used. 

Calculation of the total resistance of parallel electric circuits 

would very quickly bring the student to see the effectiveness of the 

1 ^1 

computer. The reciprocal relationships ~ = Z — become very burden- 

R^ K. 

some handwork when the r’s are large or non-iiiteger, or when n is large. 

The summation of these unit fractions is customarily tackled by the 
student as if they were ordinary fractions. The student often laboriously 
finds a common denominator and proceeds in that fashion. 

A series of such problems, starting with small integer values 
for r and proceeding to become more computationally lengthy, should 
quickly bring the student to his knees begging for computer help. 
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When the student wants the computer he is more likely to appreciate it, 
and the motivational level for learning is greater. 

In addition, the availability of a computer would make work with 
more complex circuitry possible. These are generally enjoyable to the 
student and can serve as springboards into other concepts and topics. 

Even if hardware computer is not available, the software of the . 
flow chart and program for this would be illuminating. See Figure 6 
for the Flow Chart and Program Notes. 

Finding the area under curves (i.e., the measure of the interior 
of the figure bounded by y = f(x), x = a, x = b, and y = 0) comes up 
in the study of distance, velocity, acceleration, and. time inter-rela- 
tionships, Customarily this problem is solved by the rectangular step 
method. If many steps are used (a small increment selected), or if 
the function under discussion is non-linear, the arithmetic work quickly 
becomes quite lengthy. For this reason, students (and teachers) never 
extend their work to the place where the results would be considered 
satisfactory. Furthermore, the student is apt to miss the whole basic 
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Program 



Flow Chart 






^ R1,R2, RN 



D = 


1 

Rx 










No 








Rx.S 



Yes 



10 READ N 
20 For X = 1 to N 
30 READ R(x) 

40 LET D - D + 1/R(x) 

50 NEXT X 
60 LET S - 1/D 
70 PRINT S 

80 DATA 5, 2, 3, 4, 5, 6 
90 DATA 8, 2, 4, 6, 8, 10, 12, 14, 16 
100 TO TO 10 
no END 



Notes 



Line 10 reads 5 into the machine <> 
Lines 20 and 30 pick up the resis- 
tances 2, 3, 4, 5, and 6 from the 
data in line 80 . Lines 40 through 
60 perform the actual calculations <> 
Line 100 tells the computer to go 
back and pick up a new set of data. 

Figure 6 



concept of setting an increment which approaches a limit — and then is 
denied a beautiful insight into the concept and workings of integral calculus. 

In the absence of hardware, the flow chart alone would be of some 
benefit in bringing about some of the objectives. Students should be 
allowed to construct their own flow charts and should be encouraged to 
build in a loop which runs the program through with successively smaller 
increments (greater number of rectangle steps). This would tend to 
promote a feeling for convergence to the current value. See Figure 7, 

A computer program for projectile motion could be a useful tool 
in studying some of the important relationships in kinematics. Given 
various muzzle velocities and angles we can compute the missile’s time 
in air, its maximum height, and its range. The calculation for such 
a problem can be very time-consuming, which detracts from gaining over- 
all knowledge of the kinematics involved. This is particularly true if 
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4x( .l)x/n^-x^ 



> 




Sectof> 1 is divided into rectangles 
of width .1 arid length equal to the functi|mal 
value at incremented x values. 4 times the 
summation of these rectangles would give the 
total area of the circle. 



Figure 7 



we want to make numerical comparisons for a nitnber of muzzle velocities 
and angles. Included in this particular program are: (1) the problem 

of converting from degrees to radian measure . (the computer uses only- 
radian measure); (2) the vector relationship of the muzzle velocity to 
its vertical and horizontal components — incorporated in the program by 
use of the sine and cosine functions. The computer program is straight 
forward and probably of more value than a flow chart. Included are notes 
with the program (see Figure 8) explaining the various variables and ideas 
involved. 
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Projectile Program 



100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 



(V IS MUZZLE VELOCITY, M, N FRACTIONAL 
PARTS OF n TO CHANGE FROM DEGREES TO 
RADIANS IN LINE 110) 

120, 130 GET HORIZONTAL AND VERTICAL 
VECTOR COMPONENTS OF V 
140-170 ARE THE MOTION FORMULAE 

PRINT "V2="V2, "V1="V1, "T="T, "D="D, "H="H, "A="A, 

PRINT "V="V 
GO TO 100 

DATA 1000, 1, 4, 500, 1, 4, 1000, 1, 2, 1000, 2, 3, 1000, 0, 2, 1000, 
DATA 2, 0, 

END 



READ V, N, M 
LET A=N*3.14/M 
LET V2=V*SIN(A) 
LET Vl=V*COS(A) 
LET T1-V2/32 
let T=2*T1 
LET D=V1*T 
LET H-16*T1*T1 



RUN 

WAIT. 



HEP12* 14:23 PX TUE 07/18/67 



V2= 706.825 


Vl= 707.388 


T- 44.1766 


D= 31250 


H- 7806.28 


A= .785 


V= 1000 








V2= 353.413 


Vl= 353.694 


T= 22,0883 


D= 78125 


H= 1951.57 


A= .785 


V= 500 








V2= 1000. 


Vl-0 


T= 62.5 


D= 0 


H= 15625. 


A= 1.57 


V= 1000 








V2= 866.556 


Vl=-499.08 


T= 54.1597 


D=-27063.3 


H= 11733.1 


A= 2.09333 


V= 1000 








o 

11 

(N 

> 


Vl= 1000 


II 

o 


D= 0 


o 

11 


> 

II 

O 


V= 1000 








V2=-4. 65661 E-7 




Vl=-4. 65661 


E-7 


T=-2. 91038 



D=2 0 
V= 1000 



H= 3.38813 E-15 



A= 5.78960 E 76 



Figure 8 
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Introduction to Teacher 

a) TOPIC: G.C.F. approached through prime factorization, factors 

and Euclid’s Algorithm. 

b) LEVEL OF EXERCISE: Seventh grade through Algebra I. 

c) BASIC KNOWLEDGE NEEDED: 

^ability to write a simple program with flow chart. 

*commands such as LET, PRINT, GO TO, READ, END SQR(A), INT(A), 
FOR-NEXTo 

d) SUGGESTED INSERTION INTO CLASS WORK: after study of primes, 

factors, and G.C.F. has been covered in class. 

e) APPROACH USED: The discovery approach is suggested for this unit. 

By asking the correct questions and showing examples, the teacher should 
try to guide the students into a knowledge of what they need to know and 
need to change. To find the G.C.F. we have approached the program through 
a knowledge of factors and have finished by using Euclid’s Algorithm. 

Given as supplementary topics are the programs for finding prime 
factorization of any number and finding the L.C.M. for any two positive 
integers . 



I. PROBLEM: To find the greatest common factor of any two positive 

integers A and B using both the factor approach and Euclid’s Algorithm 
approach. This is designed to be a class project to write and rewrite the 



programs together. Also included at the end are the flow charts, example 
boxes, suggested homework assignments and quizzes. 

The class might be asked how they would program a simple division 
problem, A/B. They should be encouraged to experiment with a trial and 
err or ■ approach. They might start thus: 



PROGRAM I: 



10 READ A,B 
20 LET Q=A/B 
30 PRINT Q 
40 G0 T0 10 
50 DATA 5, 2, 7, 1,7, 3 
60 END 



SANl 11:21 PX TUE 08/22/67 



2.5 

7 

2.33333 

0UT 0F DATA IN 10 
TIME: 0 SECS. 



At this point the teacher might put example boxes on the board for A,B, and 
Q, and show what the computer will find. (Fig. 1) Someone should discover 
that the computer automatically calculates out to 6 digits. We would like 
to have the whole number and the remainder. Here the teacher should 
explain the integer function of the machine. Then the students should 
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start again, inodif 3 ''ing program 1. Discussion on how to find the remainder 
of a division problem should result in step 30. 

PROGRAM II: 

11:29 



10 READ A,B 
20 LET Q=INT(A/B) 

30 LET R=A-Q*B 
40 PRINT R,Q 
50 G0 T0 10 
60 DATA 10,3,7,2,8,5 
70 END 



SAN 2 



PX TUE 08/22/67 



1 

1 

3 



0UT 0F DATA IN 
TIME: 0 SECS. 



3 

3 

1 

10 



Again the teacher shoula put boxes on the board (Figure 2) and show what 
the computer vjill do with this program. At this point the teacher could 



ask how to find factors of a number. If B is a factor of A , the 



remainder of the division A/B is zero. If you wanted to find all factors 
of A , you could start with B= 1 and check all integers up to A , Any 
whose remainder is zero should be printed out. The students should either 
be familiar with or be encouraged to think of how to get the B’s to 
increase by one each time. 

NOTE: in step 30, B+1 is now the new value of B , 
This is not the same as B+1 = B 



10 LET B=1 
20 PRINT B 
30 LET B=B+1 
40 GO TO 20 
50 END 



Combine the above with program 2 and add a check for R=0. (Step 50) You 
should have a program for finding factors of A . 

• • ■ PROGRAM III. 



10 -READ A 
20‘ LET B=1 
■30 LETQ=INT(B/A) 
40- LET R=A-B*Q 
50 IF RoO' THEN 70 
60 PRINT B,Q,R 
70 LET B=B+1 
80 GO TO 30 
90 DATA 12 
100 END 



NOTE: 



to denote ’not equal to’ in computer 
language it is necessary to use < 
and > in that order. 



By use of examples on the board (Figure 3) the students may note that this 
program finds all the factors but continues on checking the next larger 
B. This is called an infinite loop. We need a decision box to shut B 
off when it gets to be the value of A . (Step 80 below) 

Program IV is adequate for finding factors of one A . For finding 
factors of several numbers, insert this step: 

95 GO TO 10 

and change step 80 to read: 

80- IF B>A THEN 95 
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PROGRAM IV: 



10 READ A 
20 LET B=1 


SAN4 


11:35 


PX TUE 08/22/67 


30 LET Q=INT(A/B) 


1 




12 


40 LET R=A-B*Q 


2 




6 


50 IF RoO THEN 70 


3 




4 


60 PRINT B,Q 


4 




3 


70 LET B=B+1 


6 




2 


80 IF B>A THEN 110 
90 G0 T0 30 
100 DATA 12 


12 




1 


110 END 


TIME: 


0 SECS. 





When the program is out of data it will tell you so . Note also in program 
4 that by printing B and Q we got a list of factors twice. The factors 
may be found by checking from B to SQR(A) and printing out B and the 
quotient Q. This will give you all the pairs of factors with no repetition 
as you found in Figure 4. The modified program will be: 

PROGRAM V: 



SANS 11:42 PX TUE 08/22/67 



1 

2 

3 



12 

6 

4 



0UT 0F DATA IN 10 



TIME: 0 SECS. 



10 READ A 
20 LET B=1 
30 LET Q=INT(A/B) 

40 LET R=A-B*Q 

50 IF RoO THEN 70 

60 PRINT B,Q 

70 LET B=B+1 

80 IF B>SQU(A) THEN 100 

90 G0 T0 30 

100 G0 T0 10 

110 DATA 12 

120 END 



The last change to be suggested for finding factors of a number is the FOR, 



NEXT statement, thus eliminating the decision box which is step 80 in 
program. 5. The change would be step 20 to read: 



20 FOR B=1 TO SQR(A) 



and step 70 would be changed to: 

70 NEXT B 

Step 80 and 90 would be left out. (Flow chart given at end. Figure 5) 
PROGRAM VI: 



10 READ A 

20 F0R B=1 T0 SQR(A) 
30 LET Q=INT(A/B) 

40 LET R=A-B*Q 
50 IF RoO THEN 70 
60 PRINT B,Q 
70 NEXT B 
80 G0 T0 10 
90 DATA 12 
100 END 



SAN6 11:47 PX TUE 08/22/67 



1 

2 

3 



12 

6 

4 



0UT 0F DATA IN 10 



TIME: 0 SECS. 



(The next section could be assigned as homework). 

You could now ask the students if, with prgram 6, they could find the 
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G.C.F. of two positive integers. They should see that they could enter two 
numbers in the data and look at the factors of each. Then find the common 
factors and simply take the largest. Then ask the students if they could 
find the G.C.F. a little more efficiently by using the division process 
explained in Euclid’s Algorithm. With a couple of examples on the board, 
the teacher could point out what changes are needed in program 6 to find 
the G.C.F. by this second approach. We need to read in both A- and B 
in this approach. Therefore B isn’t. going to be incrementing and we 
eliminate two steps. In the second division, the divisor B is changed to 
be the remainder (step 60) and the dividend A is changed to be the 
divisor B‘. (Step 50) This repeats until a zero remainder is found. 

When the remainder is zero, the last divisor B- is the G.C.F. (Step 80), 
(Figure 6). 



NOTE: that if A < B the computer will go through one more cycle to 

divide and end by placing A=B and,B=R thus switching the order.. 



11:53 PX TUE 08/22/67 



PROGRAM -VII: 

• 10- -READ A,B 
20- LET- Q=INT(A/B) 

- 30- LET- R=A-Q*B 
40‘ IF- R=0- THEN 80 
50- LET A=B 
^ 60- LET B=R 
'■ 70' G0- T0.2O 
80' PRINT B 
■ 90' G0 T0 10 

100- DATA- 12,16,5,8,48,30 

110- END TIME: 0 SECS. 

Ybu\may' want' to make the print step a little more detailed and insert a 
line: 



SAN7 

4 

1 

6 



0UT 0F DATA-IN 10 



75 PRINT' “THE GREATEST COMMON FAGT0R{= ” 

You- may also' want' a copy of your original A and B . To get this you 
would' have' to- insert a step 15 : 

15 PRINT- A,B 

The' flow- chart for this program is -given in Figure 7. 

I I i SUGGESTED HOMEWORK : 



1) Write the program for G.C.F. using Euclid’s Algorithm. 

2') After finishing the G.C.F. orally in class, the teacher could 
have the students write a program for L.C.M. of any two positive integers. 

3') Write a modified program to print all the factors of 48 and 72 
a:nd'also'to print ’’the greatest common factor of 48 and 72=”, and then 
print' the' G-,e.F. (this would combine program 6 and program 7 and could be 
assigned' after- the entire seven programs have been covered in class.) 
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III. SUGGESTEL QUIZ: 

1) Write a program using the integer function to find the greatest 
integer ^ the given numbers : 

a) 8.312 

b) -4.9 

c) 7.6001 

2) Write a program to find the squares of the first five primes. 

3) Write a program to divide 85 by all the numbers 5 through 12 and 
print out the quotient and the remainder in each case. (Use a FOR-NEXT 
statement in your program). 

4) Write a program to find the set of multiples of any number A . 

Use A=9, 13, 25. (NOTE: the set of multiples is infinite so find at 

least ten multiples and then print the three dots . ) 

IV. SUPPLEMENTARY PROGRAM FOR L.C.M. of two positive integers. This 
exercise could be used as supplementary material or could be used at a 
later time. (The student needs to understand the DIM statement and sub- 
scripted* variables) . 

PROGRAM- VIII: 

LISTNH 

10 D-IM M(200), N(200) 

20- READ N,M 

30- F0R- K=l- T0 100 



40- 


LET N(K+1)= 


=N*K 


50 


NEXT K 




SO- 


F0R-J=1 T0 


100 


TO- 


LET M<J+1): 


=M*J 


O 

GO 


NEXT J 




90- 


F0R- 1=1 T0 


100 



100 F0R- S=1 T0 100 
no IF N(I)=M(S) THEN 155 

- 120 NEXT S 
140 NEXT I 

150 PRINT- "N0- C0MM0N MULTIPLE YET*' 
155- LET- Y=M(S) 

-160 PRINT- N,M,’'LCM=",Y 

- 170' G0 T0 20 

180 DATA 12,16,130,25,52,104,3,72 
190 END 



SAN23 


10:39 


PX THU 08/24/67 




12 


16 


LCM= 


48 


130 


25 


LCM= 


650 


52 


104 


LCM= 


104 


3 


72 


LCM= 


72 



0UT- 0F DATA IN 20 ' 

v! SUPPLEMENTARY PROGRAM FOR FINDING THE PRIME FACTORIZATION OF ANY 
GIVEN POSITIVE INTEGER: It is assumed that the students have an under- 

standing of prime numbers and factorization. The teacher can then use 
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this program to introduce prime factorization. It is suggested that this 
take two class periods. 

The first period the teacher should introduce the concept of prime 
factorization. The following are questions suggested for a teacher to 
use to begin the class discussion. First ask, "What is meant if we are 
asked to prime factor a number like 27?" After the students are sure what 
is meant by prime factoring, then the.. next question the teacher may ask is 
"How would we prime factor 27?" When it is well established that to 
prime ■ factor a number we need to use division and the divisors must all be 
primes, then the class should be given some exercises to work where they 
have to prime factor given numbers by hand. This will help them realize 
■what processes the computer must use. 

*. : - The second class period the actual program may be written. Again it 

. ... is.. suggested .that the teacher ask questions to help the students get the 
.program started.. To get the student response that first it is necessary to 
READ the number that is to be prime factored into the machine, the teacher 
-..may simply ask, "Considering the exercises that you did yesterday and 
trying to relate them to what the computer will do, what must our first 
step be?" "Also after working some exercises the students will be aware 
ir that '.to prime factor a number it is necessary to use division. The 
teacher's next question could be: "Now that we have N in the machine 

what' .must - we .do to it?". After they have responded that division is the 
operation needed, ask the students, "What kind of numbers should we use 
as divisors?" Then ask, "How would we get the prime divisors in the 
- computer?"' To' actually develop: 

40.' FOR K= 2 TO N-1 NOTE: Refer to program on following page 



100- NEXT K 

To introduce the primes into the computer will probably require a suggestion 
from the teacher. However, the teacher is cautioned here to present this 
idea only as a suggestion when the students have met a dilemma and can’t 
produce an idea that will work. Then, let them discuss your suggestion to 
decide if they think it will work. 

From their knowledge of factorization the students will know that after 
we divide N we are interested in the quotient and that is the next number to 
be divided by a prime. Therefore the teacher could ask, "Now that we have 
divided N by one prime how should we proceed now to get the rest of the 
prime factors?" Statement 50 checks to see if the quotient found by dividing 
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N by a prime number K is an integer. 

50 IF NoINT(N/K)*K THEN 100 

Here again the teacher will probably have to aid the class discussion with 
a suggestion. Statement 50 is an IF... THEN command. IF the quotient is not 
an integer which makes the statement true then the computer goes to statement 
100 which increments K and the new K is tested to see if it can by a prime 
factor . 

If the quotient is an integer then K divided N evenly and K is a 
factor of N. Therefore we want to print K. Next replace N with the 
quotient N/K so we can then check to see what prime will divide it. 

It may happen that the N that is put into the computer is a prime 
number. To check for this, C=0 is put in as a counter. Each time a factor 
is printed, C is incremented by one. After each possible K is checked then 
C is compared to zero. If it is not equal to zero then N has some factors 
and is not prime. IF C=0 then N is prime so the computer prints "IS PRIME" 
next to N. 

Complete program and flow chart (Figure 8) are on the following pages. 

It is also suggested that the teacher take a number that the class wants 
prime factored and after the program is finished go' through the program 
using that number so the students may see exactly how the computer produces 
the factors. (An example follows using the Number 27). 

Also it may be of benefit to have the students look for and program 
another approach to the same problem. 

Here are some questions that may be used for a quiz after the class 
has completed the program. 

1) Why did we test primes to N-1 as factors of N and not to N or 
beyond? 

2) What would happen if we would put 1 in for N? 

3) When we divided N by a prime number. K we ,got a quotient. Why was 
it necessary to use the INT command to see if that quotient was an 
integer? 

4) After N is replaced by the quotient N/K, we check N/K without 
incrementing K. Why? Can you give' an example? 
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INT (.3333) = 0 

1 7 ^ 0.1 

THEREFORE TAKE THE NEXT K 
1 V 4 = .25 
INT (.25) = 0 
1 ? 0.1 

THEREFORE TAKE THE NEXT K 



THEREFORE TAKE NEXT K 

1 V 26 = .038 

INT(.038) = 0 
1 i 0.1 
C = 3 

THEREFORE END 
PRINTED 

27 3 3 3 



CU*DENVER**1 . 

SORRY- HAD T0 REL0AD WITH' SUNDAYS FILE 
‘ 0N^AT 12:58 PX TUE 08/22/67 TTY 3 

■ USER- NUMBER^-P61000 
• ' ■ ' SYSTEM--BASIC 

" NEW' 0R' 0LD-.-NEW 
- NEW-- PR0BLEM- NAME'=^^PRIME 
READY. 

• IP-^APE 
• • • ' ■ READY. 

• 10- READ N 
20 PRINT N; 

• 30 LET C = 0 

■ - - ■ 40- F0R X' =■ 2 T0 N-1 

50 IF N‘ <>• INT(N/K).*K- THEN 100 

■ ■ 60‘ PRINT K; 

70 LET^ e = C+1 

' • 80- LET- N = N/K 
' • 90 G0 T0 50 
100- NEXT K 
110- IF CoO THEN 130 
120 PRINT- "IS PRIME" 

■ ■ ’ 130 PRINT 

140 G0 T0 10 

150' DATA 45 , 89,23,48,120,10,12045 

160 END 

KEY 

READY. 

RUN 



PRIME 13:00 PX TUE 08/22/67 



45 

89 

23 

48 

120 

10 

12045 



3 3 

IS PRIME 
IS PRIME 
2 2 

2 2 

2 5 

3 



2 

2 



2 

3 



11 



3 

5 

73 



0UT 0F DATA IN 10 



TIME: 10 SECS. 



BYE 



:* 0FF AT 13:01 PX TUE 08/22/67, 



/ 



ONE: 

^7 V 2 = 13.5 
■INT (13.5) = 13 
27 + 13 * 2 

THEREFORE TAKE THE NEXT K 
27 T 3 = 9 
INT(9) = 9 
27 = 9*3 

THEREFORE PRINT 3 
LET N = 9 
9 V 3 = 3 
INT (3) = 3 
9 = 3*3 
THEREFORE PRINT 3 
LET N = 3 
3 V 3 = 1 
INT(l) = 1 
3 + 3X1 

THEREFORE PRINT 1 
LET N = 1 _ 

1 V 3 = .3333 



EXAMPLE 



O 
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DATA EVALUATION FOR THE "CORRECT ANSWER" 
FOR JUNIOR HIGH SCHOOL STUDENTS 
by Margaret Ann Grill 
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Purpose of Paper 

The purpose of this paper is to present a procedure for evaluation 
of student data obtained from experiments in ninth grade physical science. 
Since tedious numerical methods are often used in such an evaluation, 
much time can be saved by using the computer; and this will be illustrated. 
Introduction 

A junior high student cannot expect nor be expected to get accepted 
literature values (those found in standard references such as the "Hand- 
book of Chemistry and Physics") from his experimental endeavors or those 
values which are calculated from a given mathematical statement. If his 
results are compared to these values, the result is nothing more than 
the calculation of his percent error. However, he is eager to know if he 
has "the correct answer." An accepted "correct answer" for a junior 
high student should be based upon his: 

1. Limited number of trials. 

2. Limited technique. 

3. Limited experimental procedure c 

4. Limited equipment and materials. 

5. Limited mathematical ability. 

To determine the "correct answer" with the given limitations, the student 
will gain valuable experience if he participates in the evaluation of the 
data which will result in the "correct answer." If his results are outside 
of the accepted bounds, perhaps he will be motivated to analyze his work 
and find out how, where and why he "goofed." He should be encouraged to 
repeat parts of the experiment to verify his analysis. 

Procedure 

To illustrate evaluation of student data, it is assumed that all 
the students have completed the experiment and that the teacher has a 
record of all the data. Transparencies are prepai-ed prior to the class 
discussion. Each set of data to be evaluated will require one transparency. 
To facilitate the discussion, the data should be recorded on the upper 
portion of the transparency. This could be done by a student assistant. 

See Figure 1 for an example. Results ("answers") from an unsophisticated 
study of the inverse square law in ninth grade physical science will be 
used as data. To start the discussion, ask the students to search the data 
and find the smallest momber. Use this number as a label for the extreme 
left column of the graph on the transparency. Continue with the next 







TRANSPARENCY EXAMPLE 



Data: 



4.6, 


5.0, 


4c8, 


5c0, 


4.9, 


5.1, 


5 a J. , 


4.9, 


5.0, 


5.2 


5.0, 


5.0 , 


4.7, 


5.2, 


5cl, 


4.7, 


4.8, 


5.0, 


5.2, 


5.1 


4.9, 


6.0, 


5.5, 


4.5, 


5.0, 


4.9, 


5.6, 


4.8, 


5.5, 


4.8 



(unit in centimeters) 



M.5 4.6 4.7 4.8 4.9 5..0 5d 5.2 5.5 6..0 
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largest number, etc., until the columns of the graph are labeled in order 
from the smallest number at the left to the largest number on the right. 

See Figure 1 for an example. To enter the data on the graph, have a 
student read the first number. Mark an "x" to represen c this number in 
the bottom space of the appropriate column. If this number occurs again, 
it will occupy the next space up. After all numbers are represented on 
the graph, the occupied spaces could be filled in with a colored pen. 

The result would be a bar graph or a hologram. At this time, ask the 
students to recall the units for the numbers. For this example, the unit 
is the centimeter. Have them consider how much one number varies from the 
next, which columns have the most numbers, which have the least, etc. From 
this, they should arrive at some agreement on which numbers occur most 
frequently. They will probably select the numbers from 4.7 to 5.2. A 
discussion should follow deciding whether these "answers" are reasonable 
for the procedure of the experiment. Probably these are the "correct 
answers" for tnis one set of data. 

Ask the students how they would find the average (mean) of the num- 
bers. Introduce the concept of deviation from the mean. Guide the dis- 
cussion in such a way that they will understand the meaning of standard 
deviation. After they have some idea of this, present the mathematical 
expression for the calculation of the standard deviation. Ask how many 
of them would like to make these calculations for the nine sets of data 
from this experiment. Have them estimate how long these calculations 
would take. If a computer is available, introduce its use at this time 
for these repeated calculations. Put the tape (which has been prepared 
prior to the discussion) of the program into the computer, and run it. 

Note the calculation time for the computer, and compare this with the 
previous estimation by human calculations. See Figure 2 and following 
pages for this program, etc. 

Explain that the computer can do only that which it is instructed 
to do in its own language. Put a transparency of the flow chart for the 
program on the projector, and briefly explain how the program is logically 
executed. 

Use the standard deviation obtained from the computer calculation to 
determine the acceptable range for the "correct answers." For such small 
numbers, it would be best to allow only one standard deviation. The 
"correct answers" for this set of data would be from 4,7 to 5.3. Have the 
students compare these with the results obtained from the hologram. It 
would be well to have a short discussion of rounding off numbers and the 
possible errors involved in general. 
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COMPUTER PROGRAM 

FOR SORTING A LIST OF DATA, CALCULATION OF MEAN, 
DEVIATIONS AND STANDARD DEVIATION 
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GRILLl 15:03 PX THU 07/13/67 



10 DIM N(30), A(30), D(30) 

20 PRINT "N=", "D=" 

30 READ N 

40 F0R I = 1 T0 N 

50 READ A(I) 

60 NEXT I 

70 LET K = 1 

80 LET J = K + 1 

90 IF A(K) < = A(J) THEN 130 

100 LET T = A(K) 

110 LET A(K) = A(J) 

120 LET A(J) = T 
130 LET J = J + 1 
140 IF J <= N THEN 90 
150 LET K = K + 1 
160 IF K < N THEN 80 
200 LET S = 0 
210 F0R I = 1 T0 N 
220 LET S = S + A(I) 

230 NEXT I 

240 LET M = S/N 

250 LET T = 0 

260 F0R I = 1 T0 N 

270 LET D(I) = A(I) - M 

280 LET T = T + D(I)*D(I) 

290 NEXT I 

300 LET T1 = T/N 

310 LET T2 = SQR (Tl) 

320 F0R I = 1 T0 N 
330 PRINT A(I), D(I) 

340 NEXT I 

350 PRINT "M=" M, "S.D.=" T2 
355 PRINT 
360 G0 T0 20 
660 END 



In general, the listed statements 
are instructing the computer to do 
the following: 

10-60 Put data into memory. 

70-160 Sort the random numbers . 

200-230 Calculate the sum. 

240 Calculate the mean. 

250-270 Calculate the deviation 

of each number from the 

mean. 



280-310 Calculate the standard 
deviation . 

320-350 Print sorted numbers with 
the deviation of each, 
mean and standard deviation. 

355 Allow a blank space before 

printing again. 

360 Pick up next set of data. 

Data statements would be ntmbered 

from 361 to 659. 

660 Stop computing . 



The example below illustrates the data instruction using the first set 
of data. 



365 DATA 30 

370 DATA 4. 6, 5. 0,4. 8, 5. 0,4. 9, 5. 1,5. 1,4. 9, 5.0, 5.2, 5.0, 5.0, 
380 DATA 4.7, 5. 2, 5.1, 4.7, 4.8, 5,0, 5.2, 5.1, 4.9, 6.0, 5.5, 
390 DATA 4.5, 5.0, 4.9, 5.1, 4.8, 5.5, 4.8 




I 
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COMPUTER OUTPUT 
FOR THE NINE SETS OF DATA 



GRILLl 



11:15 



PX THU 07/13/67 



N= 

4.5 

4.6 

4.7 

4.7 

4.8 
4.8 
4.8 

4.8 
4c9 

4.9 
4.9 
4.9 
5 

5 

5 

5 

5 

5 

5 

5.1 

5.1 

5.1 

5.1 

5.1 

5.2 
5.2 
5.2 
5.5 
5.5 

6 

M= 5.01333 



D= 

-.513333 
-.413333 
-.313333 
-.313333 
-.213333 
-.213333 
-.213333 
-.213333 
-.113333 
-.113333 
-.113333 
-.113333 
-1.33333 E-2 
-1.33333 E-2 
-1.33333 E-2 
-1.33333 E-2 
-1.33333 E-2 
-1.33333 E-2 
-1.33333 E-2 

8.66667 E-2 

8.66667 E-2 

8.66667 E-2 

8.66667 E-2 

8.66667 E-2 
.186667 
.186667 
.186667 
.486667 
.486667 
.986667 
S.D.= .286046 



8 


-.84 




9 


.16 


8.2 


-.64 




9 


.16 


8.3 


-.54 




9 


.16 


8.4 


-.44 




9 


.16 


8.4 


-.44 




9 


.16 


8.4 


-.44 




9.2 


. .36 


8.5 


-.34 




9.2 


.36 


8.5 


-.34 




10 


1.16 


8.5 


-.34 




10 


1.16 


8.6 


-.24 . 




10 


1.16 


8.6 


-.24 




M= 8.84 


S.D.= 


8.7 


-.14 








8.7 


-.14 








8.8 


-3.99999 


E-2 






8.8 


-3.99999 


E-2 






8.8 


-3.99999 


E-2 






8.9 


CO 

o 

• 








8.9 


CO 

o 

• 








8.9 


CO 

o 

o 








8.9 


CO 

o 

o 









N= 


D= 




6.5 


-.673333 


$ 

1 


6.8 


-.373333 


1 


6.8 


-.373333' 


>( 


6.9 


-.273333 




6.9 


-.273333 


a 


6.9 


-.273333 


1 


6.9 


-.273333 




7 


-.173333 




7 


-.173333 


y 


7 


-.173333 


li 


7 


-.173333 


« 

25 

I 


7 


-.173333 


V 

$ 

a 


7.1 


-7.33333 E-2 




7.1 


-7.33333 E-2 




7.1 


-7.33333 E-2 




7.1 


-7.33333 E-2 


§ 

* 


7.1 


-7.33333 E-2 


i 


7.2 


2.66667 E-2 




7.2 


2.66667 E-2 




7.2 


2.66667 E-2 




7.2 


2.66667 E-2 


1 


7.4 


.226667 


■S 


7.4 


.226667 


1 


7.4 


.226667 


I 


7.5 


.326667 




7.5 


.326667 


s 


7.5 


.326667 




7.5 


.326667 


1 

£5 


8 


.826667 


i 


8 


.826667 


I 


1= 7.17333 


S.D.= .324482 


1 

1 

1 

X 



s 

I 



I 



I 
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N= 


D= 






N= 


D= 


9.5 


-.876667 






10.8 


-.82 


9.7 


-.676667 






11 


-.62 


9.8 


-.576667 






11 


-.62 


9.8 


-.576667 






11 


-.62 


9.9 


-.476667 






11.1 


-.52 


9.9 


-.476667 






11.2 


-.42 


10 


-.376667 






11.2 


-.42 


10 


-.376667 






11.2 


-.42 


10 


-.376667 






11.2 


-.42 


10 


-.376667 






11.2 


-.42 


10 


-.376667 






11.2 


-.42 


10 


-.376667 






11.2 


-.42 


lOcl 


-.276667 






11.3 


-.32 


lOcl 


-.276667 






11.3 


-.32 


10.1 


-.276667 






11.3 


-.32 


10.1 


-.276667 






11.3 


-.32 


10.2 


-.176667 






11.3 


-.32 


10.2 


-.176667 






11.4 


-.22 


CO 

4 

o 
1 — 1 


- 7.66666 E -2 






11.4 


-.22 


10.3 


- 7.66666 E -2 






11.5 


-.12 


10.3 


- 7.66666 E -2 






11.5 


-.12 


10.3 


- 7.66666 E -2 






11.7 


8.001 


10.4 


2.33334 E -2 






12 


.38 


10.5 


.123333 






12 


.38 


10.5 


.123333 






12.2 


.58 


11.2 


.823333 






12.5 


.88 


11.5 


1..12333 






12.7 


1.08 


12 


1.62333 






13 


1.38 


12 


1.62333 






13.2 


1.58 


12.6 


2.22333 






13.7 


2.08 


«= 10.3767 


S.D.= .725113 






M= 11.62 


S.D.= 


M= 


D= 










12 


-.813333 


14 


1.18667 






12 


-.813333 


14 


1.18667 






12 


-.813333 


14.8 


1.98667 






12.1 


-.713333 


14.9 


2.08667 






12.1 


-.713333 


15 


2.18667 






12.1 


-.713333 


M= 12.8133 


S.D.= .882698 






12.1 


-.713333 










12.2 


-.613333 










12.2 


-.613333 










12.2 


-.613333 










12.3 


-.513333 










12.3 


-.513333 










12.3 


-.513333 










12.3 


-.513333 










12.4 


-.413333 










12.6 


-.213333 










12.6 


-.213333 










12.6 


-.213333 










12.7 


-.113333 










12.7 


-.113333 










12.9 


8.66668 E -2 










13.1 


.286667 










13.1 


.286667 










13.2 


.386667 










13.6 


.786667 













N= 


D= 






N= 


12.9 


-1.08333 






13.8 


13 


-.983333 






14 


13.1 


-.883333 






14 


13.1 


-.883333 






14.1 


13.1 


-.883333 






14.1 


13.2 


-.783333 






14.1 


13.2 


-.783333 






14.1 


13.2 


-.783333 






14.1 


13.3 


-.683333 






14.2 


13.4 


-.583333 






14.3 


13.5 


-.483333 






14.3 


13.5 


-.483333 






14.3 


13.5 


-.483333 






14.3 


13.5 


-.483333 






14.5 


13.6 


-.383333 






14.5 


13.8 


-.183333 






14.5 


13.8 


-.183333 






14.5 


13.8 


-.183333 






14.6 


13.8 


-.183333 






14.6 


14 


1.66667 E-2 






14.7 


14 


1.66667 E-2 






14.8 


14.1 


.116667 






15 


14.8 


.816667 






15.5 


15 


1.01667 






15.8 


15.1 


1.11667 






16 


15.2 


1.21667 






16 


15.2 


1.21667 






16.1 


15.7 


1.71667 






16.7 


16 


2.01667 






16.7 


16.1 


2.11667 






16.8 


M= 13.9833 


S.B.= .926673 






M= 14 


N= 


D:= 








14.7 


-.966667 


16.9 


1.23333 




14.8 


-.866667 


17 


1.33333 




14.9 


-.766667 


17 


1.33333 




14.9 


-.766667 


17.4 


1.73333 




15 


-.666667 


17.5 


1.83333 




15 


-.666667 


17.5 


1.83333 




15 


-.666667 


M= 15.6667 


S.D.= .888194 


15 


-.666667 








15 


-.666667 








15.1 


-.566667 








15.1 


-.566667 








15.1 


-.566667 








15.1 


-.566667 








15.2 


-.466667 








15.2 


-.466667 








15.2 


-.466667 








15.3 


-.366667 








15.4 


-.266667 








15.5 


-.166667 








15.6 


-6.66666 E-2 








15.8 


.133333 








15.8 


.133333 








16.2 


.533333 








16.8 


1.13333 









1 



o 

ERIC 



D= 

-1.03333 
-.833333 
-.833333 
-.733333 
-.733333 
-.733333 
-.733333 
-.733333 
-.633333 
-.533333 
-.533333 
-.533333 
-.533333 
-.333333 
-.333333 
-.333333 
-.333333 
-.233333 
-.233333 
-.133333 
-3.33333 E-2 
.166667 
.666667 
.966667 
1.16667 
1.16667 
1.26667 
1.86667 
1.86667 
1.96667 
S.D.= .891191 
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Continue the evaluation of each set of data. Using the sorted 
list of numbers from the computer calculations will shorten the time for 
the preparation of more holograms. Have various students read these so 
they will see the "answers" from the computer. It would be interesting 
for the students to know that the computer sorted the random set of numbers 
into successive order starting with the smallest. Probably after the com- 
pletion of the second hologram the students will accept the results from 
the computer calculations for the remaining sets of data. 

If a computer is available in the school, students should be given 
a program for their calculations if the calculations are time-consuming. 

The use of the computer by the student would be a great motivator. 

In this experiment, two graphs are made and compared. The result 
shows that if the illumination is plotted versus the distance a definite 
curve is obtained, and that the illumination versus the inverse of the 
square of the distance results in a straight line. The more able student 
should be able to derive the equation for the straight line using; geometrical 
concepts and algebra. See Figure 3. Another computer demonstration could 
be presented to do the tedious arithmetic involved in deriving this equation 
by the method of least squares. For data from this experiment the linear 
equation is "y = 0.04x," since -0.00058 is essentially zero. See Figure 4 
and the following pages. The significance of the straight line obtained 
should be discussed. 

Summary 

A procedure involving the use of the computer to evaluate data re- 
sulting from students’ experimental work has been presented. 

The use of a computer to statistically evaluate data relieves the 
teacher of many tedious and time-consuming calculations. If this tech- 
nique were available for use on data collected from experimental work by 
students, the teacher should be more willing to do the large number of 
calculations required for this evaluation. Changes, then, in experimental 
procedure could be suggested which would result in more meaningful data 
for the student. This modified procedure would enable future students 
to better understand primary concepts resulting from experimental activities 
in introductory physical science. 
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Figure 3 
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FLOW CHART • * 

FOR METHOD OF LEAST SQUARES 



I 



i 



& 




Figure 4 
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COMPUTER PROGRAM 
FOR METHOD OF LEAST SQUARES 



Sl= 2.82 



S2= .108708 



A=-5. 87172 E-4 
.04 

1.92901 E-2 
1.29132 E-2 
9.24556 E-3 
7.43163 E-3 
6.10352 E-3 
5.10204 E-3 
4.56538 E-3 
4.05696 E-3 



S3= 1.5364 
B= .040423 



S4= 6.04501 E-2 



TIME: 



1 SECS, 



LISTNH 
10 READ N 
20 F0R I = 1 T0 N 
30 READ X(I),D(I) 

40 NEXT I 

50 F0R I = 1 T0 N 

60 LET Y(I) = 1/(D(I)*D(D) 

70 NEXT I 

80 LET SI = 0 

90 LET S2 = 0 

100 LET S3 = 0 

110 LET S4 = 0 

120 F0R I = 1 T0 N 



In general the listed statements are 
instructing the computer to do the 
following : 

10-40 Put data into memory. 

50-70 Calculate the square of the dis- 
tance and then the inverse. 

80-170 Calculate various sums, 

180 Print the sums. 

185 Allow a blank space (line) 

before printing again. 

190-210 Calculate values for the coeffi- 
cients of "x" and "y" , 

220 Print coefficients. 

224-226 Print y values. 

Data statements would be numbered from 
227-249 



250 



Stop computing. 



51 + X(I) 

52 + Y(I) 

53 + X(I)*X(I) 

54 + X(I)*Y(I) 



130 LET SI 
140 LET S2 
150 LET S3 
160 LET S4 
170 NEXT I 

180 PRINT "S1="S1,"S2="S2,"S3="S3,"S4="S4 
185 PRINT 

190 LET Z = (9*S3)-(S1*S1) 

200 LET A = ((S2*S3)-(S4*S1))/Z 
210 LET B = ((9*S4)-(S1*S2))/Z 
220 PRINT "A="A, "B="B 

224 F0R 1=1 T0 N 

225 PRINT Y(I) 

226 NEXT I 

230 DATA 9, 1.00, 5.0, .50, 7. 2, .33, 8. 8, .25, 10. 4, .20 
240 DATA 11. 6, .17, 12. 8, .14, 14.0, .12, 14. 8, .11, 15. 7. 

Since the program was listed after the calculations, the output 

is at the top. 



Note: 
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MATHEMATICAL EXPRESSIONS 






For the Mean 

• S/N • ■ • ■ 

where- M* = mean * . • • . 

..r S -= sum . . ... 

N = total number of values 

For the Deviations 

, D = A - M 

where D = deviation 
* ■ A’ = value from, data’ 

,M = mean 

For the- Standard Deviation '• ... 




SD = 



' where SD = standard deviation 
• ' ■ Eh = sum of all values from 5 = 1 to n 

t = 1 - 

n = tofal number of values 
D ■ = ' deviation for each value 



For ' the Method of Least ' Squares 

• • (i). -AN BEXj^ = Ey^^ 

( 2 ) 

/where. N* =, number of values used 

k = summations from k = 1 to N 
r , = values from’ data' ’ 

. . .-A,B =;Constants to -be .solved 

Let SI = sum of x values 

52 = sum of y values 

53 = sum of X- values 

• ■ 'S4 = sum of xy values 

Substituting in (1) and (2) above 

. - . — . . 9A: + SIB ,= -S2 . : ' 

SIA + S3B ■= S4 '■ 

Solution for A and B 

, (S2)(S3) - (S4)(S1) 

(9)(S3) - (S1)(S1) 

R - (9)(S4) - (S1)(S2) 

. ■^.(9)CS3) - (SlXSl). 

The resulting linear equation is "y = Bx + A” 
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INTRODUCTION TO FLOW CHARTINi^ 

By J. Groninger S T. Shoemaker 

A flow chart is a pictorial diagram that can help one see his way 
through any logical process, sometimes not associated with a computer problem 
or exercise. Flow charts are commonly used in many fields of science and 
business. Scientists and mathematicians use flow charts to assist them in 
preparing problems to be solved with the aid of digital computers. 

There are several translation steps necessary for preparing 
i.nformation for using computers. The task of using the computer to help 
solve problems can be separated into two steps: that of reducing the 

problem to a sequence of elementary steps (flow charting or preparing an 
algorithm); and that of converting to a procedural language. Each of these 
steps is a translation process. 




It should be remembered there is no one way of preparing a flow chart 
or program for solving a problem. Also, there is not a predetermined length 
or limit of complexity of a flow chart. The thoroughness of a translation is 
dependent on the desired depth of logic and difficulty of the program. 
Sample Flow Charts 

Several sample flow charts are shown to illustrate the intricacies of 
flow chart preparation and standardized symbols. Compare the general and 
expanded charts of one situation shown here. Both flow charts are correct. 
FLOW CHART SYMBOLS 



C 




Terminal : The beginning or end of a program. 



Calculation or Computing Symbol : Used to indicate a 

series of arithmetic operations or computer 
instructions . 




Used to indicate the input of information. 




Decision : A symbol used to indicate points where 

branching of alternate paths are possible based 
on variable conditions. 
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Output : Used to indicate an output of information. 




Connector t An entry or exit to another part of a 
program . 



Non-Math Situation 

Prepare a flow chart of starting an automobile or making a local 
telephone calll 
Mathematical Situations 

•Prepare a flow chart of an appropriate situation. 

General Mathematics — How to solve story problems or word problems.. 

Algebra I — Finding the solution set of an equation. 

Geometry — Determining the coordinates of the midpoint of a 
given segment . 

Algebra II — Finding prime numbers. 

Introduction to Flow Charting Teacher Information Sheet 

A preparation of a flow chart probably will not make the solution of a 
problem easier or faster, but instead will assist the student in understanding 
the logical steps necessary in finding the solution. In short, a flow chart 
will assist the person preparing it to ’’see" what is going on inside the 
actual solution. A flovr chart is the linking betvfeen a communicating 
everyday language and procedural language necessary for directing the 
necessary operations, often done .on a computer. Learning to do basic flow 
charting will prepare you for using procedural languages such as: BASIC, 

COBOL, EASYCODE, etc.. 

Flow chart construction demands patience and the desire to prepare a 
detailed plan for problem solution. Students may not immediately recognize 
the value of a flow chart until working on a computer or computer terminal. 
Therefore, it is advisable for students to have a "hands-on" experience 
trying the logical sequential pattern of his flow chart. 

The sample flow charts included are shown to assist the student in 
use of standardized symbols and the flow of a logically prepared diagram of 
solution planning. Therefore, we recommend you replace the examples 
provided or supplement them with situations and flow charts of your own 
choosing. 






Teaching Techniques and Recommendations 

A. Flow charts have been introduced by instructors in two ways: 

1. Through concrete problems, followed by a more abstract relation. 

2. When students are in search of a method of solution because 
other methods have failed or are not understood. 

B. The instructor must have a working knowledge, from experience, of the 
assigned situation or problem. It is necessary to understand all 
facets of a problem solution, since there is no single correct 

flow chart for a problem. 

C. On the situations to be completed herein the student is expected to 
progress at his own rate with minimal guidance from the instructor. 

D. A flannel board with flow chart symbol cut-outs is an excellent way 
of demonstrating the logic pattern of a problem solution. An over- 
head can be used in a similar manner with plastic cutouts drawn on 
the projector table and projected onto a chalk board. 

E. Flow chart templates can be purchased from many computer manufactures 
for a nominal cost . 
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f 



B C 
Id 



Start 



SAMPLES OF NON-MATH FLOW CHARTS 
Sharpening A Pencil 
(A General Flow Chart) 




Get 




Walk To 




Stick Pencil 


Permission 




Sharpener 




In Sharpener 






Turn 


► 


Take Pencil 
Out when 
Sharp 




Return to 


Crank 

Clockwise 




► 


Seat 












End 



3 




I 



(More Detailed Flow Chart) 









I 



^ star 





Turn Crank 
Clockwise 
Several 
■Times 




Stick 
Pencil ir 




Sharpen- 

er 


1 


f 




[No 



<s- 




Is 

^Pencil 
Sharp 



No 


Raise Hand 






for 






Permission 





1 


Wait 

A 




Bit 



Yes 
> 


Go To 


^ i 




Seat 


— ^ 







ERIC 
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SAMPLES OF MATH FLOW CHART! 
.AVERAGING NUMBERS 




DETERMINING WHETHER AN INTEGER 
IS EVEN OR ODD 




USING THE COMPUTER FOR 
PERMUTATIONS AND COMBINATIONS 
by Neil 0. Haflich and Ellis R. Mercer 
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To the Teacher 

The material presented in the following pages is intended to il- 
lustrate the use of the computer when working with permutations and com- 
binations. It is written with the assumption that the student has had 
no previous experience with the computer. This material should be useful 
even if no computer is actually available to the student. This material 
can be used in junior high or senior high school. 

The material is divided into three sections. The first section 
0 xpla.ins how a computer may be used in the calculation of factorials. 

It is assumed that the student is familiar with the formula ^P^ = nl. 

The student should have practice in calculating factorial. 

The second section deals with the calculation of the numbers of 
permutations of n things taken r at a time. The student should have 
had the experience of working with the formula ^P^ = 

The third section shows the calculation of the number of combina- 
tions of n things taken r at a time. The student should have had the 

n j 

experience of working with the formula = (n-r) Ir! * 

It is our recommendation that each of these sections involving 
the computer be taught following the corresponding sections in the 
chapter on permutations and combinations. For example j after teaching 
the permutations of n things taken n at a time* then teach how the com- 
puter could be used to calculate n factorial. 

This method gives the student reinforcement in both the computer 

and the mathematical computations. 

At the end of Section III it is suggested in the text that some 
students try to construct a flow chart and write a program for calculating 

C . 

n r 

One possible example is shown in Figures T-1 and T-2. It should be 

observed that in addition to solving for the program also gives the 

value for P which is an intermediate value in the calculation of C . 

nr “ ^ 

There are times when it is desirable to have data entered during 
the running of a program. This is particularly true when one person 
writes the program and enters it into the machine’s memory location 
while other persons are to supply the data. This may be done by an 
"INPUT" statement, which acts as a "READ" statement but does not draw 
numbers from a "DATA" statement. 

For example, in the program shown in Figure T-2, line 20, it reads: 
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Flow Chart for P and C 

nr nr 




Figure T-1 



102 



20 INPUT N,R 

This statement must come before N and R are used in the program. 

When the computer encounters this statement, it will t3^pe a question 
mark. The operator of the computer types two numbers, separated by a 
comma, presses the return key; and the computer goes on with the remainder 
of the program. 

Frequently an "INPUT" statement is combined with the "PRINT" 
statement to make sure the user knows what the question mark is asking 
for. The program would then read: 

10 PRINT "Y0UR VALUES OF N AND R ARE", 

20 INPUT N,R 

The machine will then type out 

Y0UR VALUES OF N AND R ARE? 

Program for P and C 
n r n r 

10 PRINT "Y0UR VALUES 0F N AND R ARE" : 

20 INPUT N,R 

3b LET Q=1 

40 LET K=1 

50 LET Q=Q*K 

60 IF N=R THEN 140 

70 IF K=(N-R) THEN 120 

■80 IF K=R THEN 150 

90 IF K=N THEN 170 

100 LET K=K + 1 

110 G0 T0 50 

120 LET T=Q 

130 G0 T0 80 

140 LET T=1 

150 LET V=Q 

160 G0 T0 90 

170 LET P=Q/T 

18Q LET C=P/V 

190 PRINT N,R,P,C 

200 G0 T0 20 

210 END 



Figure T-2 



Section ^ 
with the Compurer 

This section explains how the computer can be used to help you 
solve factorials. As you have experienced, the factorial of a small 
number results in a large number. It certainly would be nice if we 
had a machine to do the series of multiplications that is necessary to 
get a factorial even for factorials as small as 121 The machine can help 
us do this. 
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The first step in preparing ourselves to use the computer would be 
to draw a flow chart* The purpose of a flow chart is to enable you to 
see what mathematical operations are done by the computer. You can also 
see how each step is related to each other step. The flow chart could be 
thought of as a road map that you would consult in planning a vacation 
by automobile. A flow chart is not unique. As you could plan several 
roads to travel during your vacation, so you could plan several flow 
charts that would compute N factorial. One of these flow charts of N 
factorial is shown in Figure 1. 

The various symbols used in making flow charts each have a special 
meaning. The parts of the flow chart in Figure 1 on the preceding page 
are each labeled with a small letter. Refer to Figure 1 for the following 
explanation. 

The symbol at (a) has the sole purpose of showing where to begin 
using the flow chart. 

The symbol at (b) is for input. It tells you what data you need to 
supply to the computer. The machine takes the number you give it and 
puts it in a memory location which we will call N. Any time you call 
for N the machine then goes to this location to find what value you 
assigned to N. It has many memory locations for storing variables. 

You may think of a memory location as a sack with a label on it. - If' 
we are calculating 11 we would be telling the machine to go to the sack 
labeled N and put a 7 in it. 

The symbol at (c) indicates a computation box. A computation box 
indicates the instruction we must give the computer. Q = 1 inside the 
computation box means that we will tell the computer to label a sack Q 
and to put a 1 in the sack. The left hand side of the equal sign is the 
label for the memory location (sack) , and the right hand side is what is 
to be put in the location. We ask the computer to put 1 in the Q sack 
because the first factor in calculating any factorial is 1. 

The computation box at (d) means we instruct the computer to put 
2 in the memory location marked K. It will become clearer later that 
we want to put a 2 in this location because we want to calculate 21 first. 
2! will be used in the calculation of any higher factorial. We may omit 
the calculation of 01 and II as the result of both of these is 1. 

The Q = Q*K is the computation box at (e). It looks confusing until 
we remember that the Q on the left of the equal sign is simply the label 
of the sack (memory location), and the Q*K on the right of the equal sign 
is the instruction; so the statement will mean for the computer to go 
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Flow Chart N! 




Figure 1 
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to the memory location Q and to the memory location K, then to take the 
values it finds in the sacks and to multiply (*) the values together. 
Finally , it will store this new value of Q in the Q sack; and the value 
of K stored in the K sack remains unchanged at this time. 

The symbol at (f) means that we will instruct the computer to 
make a choice, so it is called a decision box. The K = N in the decision 
box means that the computer is to compare the values in the K and N 
sacks. In other words, the decision box asks the question, "Does K = N?" 
If the answer is no, the next step in the flow chart is to the computation 
box at (g)o The statement K = K + 1 in that box indicates that the 
computer will be instructed to go to the K sack, take the value of K 
that is stored there, increase it by one and put the new value of K back 
in the sack. 

The arrow from the box at (g) goes back to the box at (e). This 
means that the computer will take this new value for K and multiply it 
by the value of Q and put the new Q value in the Q sack. This cycle is 
repeated until the answer to the question "DOES K = N" is yes. 

When the answer is yes , the next step is indicated by the symbol 
at (h). This symbol is the one which is used to indicate that the computer 
will be instructed to print certain information. In the present example, 
the computer is directed to print the value of Q which turns out to be 
the value of N ! . 

Now take an example to see exactly how it would work. Let N == 4. 

This means we intend to have the computer calculate 41. Therefore, 4 
is given as the data for the computer to store in the N sack. Now- 
let Q = 1, and let K = 2. The instruction Q = Q*k will give us Q = 1*2 
which is 2. This is the value now stored in Q. The decision box now 
indicates that the computer will compare the value in the K and the N 
sacks. Since 2 does not equal 4, it will increase K by 1, which is 3, 

It now computes a new Q which will be 2*3 or 6, where 2 is the last 
preceding value for Q. Since K is now 3, which is still not equal Lo 4, 

K will be increased to 4. The new Q will be 6*4 or 24. When K is com- 
pared to N, they will be equal; so the computer will print the value of 
Q to be 24. This is 41. 
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One method of checking to see whether a flow chart will lead to 
correct results from the computer is to do what is called "playing 
computer." This involves boxes or squares to represent the memory 
location of the computer and then following the flow chart to see 
what happens to the values stored in each location. 



N 

For example, the first data location according to the flow 

4 

chart shown in Figure 1 is labeled N ; so if we supply a value 
of 4 for N we can represent it as shown. 

Q 

^ This box represents a value of 1 stored in the memory location 

labeled Q. 



K 

The next step according to the flow chart is to assign a value 
2 of 2 in sack K. This is represented as shown. 



Q 

The product of K*Q would now be 1*2 or 2, which is now the 

2 

value in the sack Q. 




Since K = 2 and N = 4 at this time, the decision box directs 
the computer to computation box. 



The computation box directs the computer to add 1 to K, which 
now has a 3 stored there. 



Q 

Then the spot labeled Q will have the value of K*Q or 2*3, 
which is 6. 







/ 

j 

I 

I 

t 



The flow chart is not in the form that gives instructions to the 

computer. Instructions to the computer must be written in a certain way 

in order for the computer to do the problem. When the instructions 

indicated in the flow chart are written in proper form, the instructions 

are called a program. The program for the flow chart shown in Figure 1 

might be written as follows: 

10 READ N 
20 LET Q=1 
30 LET K=2 
40 LET Q=Q*K 
50 IF K=N THEN 80 
60 LET K=K+1 
70 G0 T0 40 
80 PRINT Q 
90 DATA 4,8 
100 G0 T0 10 
no END 

The computer takes your program in numerical order. You should use 
statement numbers far enough apart so that you can add to your program with- 
out rewriting the entire program. For example, if you wish to add another 
"DATA” statement, you could number the statement 91; and the computer 
would fit it into your program in proper numerical order. 

The "READ" statement tells the machine what variable you want 
and to set aside a memory location for N. It then goes looking for a 
"DATA" statement. It now puts 4 in the N sack. 

"LET Q=l" tells the computer to set aside a Q sack and put 1 in 
that sack. 

"LET K=2" tells the computer to set aside a K sack and put 2 in 
that sack. 

"LET Q=Q*K" tells it to go to the Q and K sacks. Take the product 

of the two values, and put this new value in the Q location. Also put 

the K back in the K location. 

"IF K=N THEN 80" has the computer compare K and N, and if they are 

equal, then to go directly to statement 80; but if they are not equal, 

it will go to the next statement. 

"LET K=K+1" has the machine go to the K location, take K out of the 
K location, and increase it by one and then put it back in the K location. 

"G0 T0 40" tells it to go back to 40 and to go through again. 

"PRINT Q" has it print out the latest value of Q. 

"DATA" statement furnishes the values of the variables in the "READ" 

statement. 

"G0 T0 10" tells the computer to go back to statement 10. Here 
the computer takes the next value out of the "DATA" statement and assigns 
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that as the new value for N. N now has the value of 8. Now the entire 
process is repeated again. 

"END" tells the computer to quit reading the program. 

bection II 

Permutation of Things Taken ^ at £ Time 

In the preceding section we learned how a computer could be used 
in working with permutations. An example was given for the calculation 
of which is the number of permutations of 4 things taken 4 at a time. 
Since ^P^ is n! , the only thing required of a computer is to calculate nl. 
Similarly, if we now consider the problem of using a computer to solve 

P , again the primary operation is the calculation of factorials. 

n I 

Since ^P^ is t > program for the computer results in the 

calculation of n! and (n-r)! and then divides the results of these two 
quantities . 

The flow chart is given in Figure 2. This is followed by a tentative 
program for the computer. 



Program for P 
nr 



10 


READ N,R 


100 


LET T=T*K 


20 


LET 


Q=l 


110 


IF K=(N-R) THEN 140 


30 


LET 


K=2 


120 


LET K=K+1 


40 


LET 


Q=Q*K 


130 


G0 T0 100 


50 


IF K=N THEN 80 


140 


LET Q=Q/T 


60 


LET 


K=K+1 


150 


PRINT N,R,P 


70 


G0 T0 40 


160 


G0 T0 10 


80 


LET 


T-=l 


170 


DATA 5,3,10,2 


90 


LET 


K=1 


180 


END 



If we follow the flow chart and the program, it is evident that the 
first part is simply the calculation of N! , which is the same as the pre- 
ceding section. The last number to be stored in the memory location 
labeled Q will be N!. The next part of the flow chart is for the cal- 
culation of (N-R)!, and if this section is studied, it should be evident 
that it is very similar to the first section. It should be clear that 
the memory location labeled T is being used to store the value of CN-R)! 
as it is being generated. One difference in the two parts of the program 
is that in calculating N I , the first step gives 2 1 , while in the calcula- 
tion of (N-R) I the first step gives ll. This is necessary to take care 
of the cases where R is one less than N. If the program is used to 

"play computer" for the case of _P. the results could be written as 

b o 



follows : 
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Flow Chart P 
n r 




Figure 2 



Ill 



Line Values in Various Memory Locations 



Number 


N R 


Q 






K 


T 


p 


10 


5 3 


? 






? 


9 

• 


9 

# 


20 


5 3 


1 






? 


? 


9 


30 


5 3 


1 






2 


? 


9 


40 


5 3 


2 






2 


? 


? 


50 


(G0 T0 60 because K 


does 


not 


equal 


N) 






60 


5 3 


2 






3 


? 


? 


70 


(G0 T0 40) 














40 


5 3 


6 






3 


9 


? 


50 


(G0 T0 60 because K 


does 


not 


equal 


N) 






60 


5 3 


6 






4 


? 


? 


70 


(G0 T0 40) 














40 


5 3 


24 






4 


? 


? 


50 


(G0 T0 60 because K 


does 


not 


equal 


N) 






60 


5 3 


24 






5 


9 


? 


70 


(G0 T0 40) 














40 


5 3 


120 






5 


? 


? 


50 


(G0 T0 80 because K 


does 


equal N) 








80 


5 3 


120 






5 


1 


? 


90 


5 3 


120 






1 


1 


? 


100 


5 3 


120 






1 


1 


? 


110 


(G0 T0 120 because K 


does 


not 


equal N-R) 






120 


5 3 


120 






3 


1 


? 


130 


(G0 T0 100) 














100 


5 3 


120 






2 


2 


? 


no 


(G0 T0 140 because K 


does 


equal N-R) 






140 


5 3 


120 






2 


2 


60 



At this point the next instruction to the computer is to print 

the values of N, R and P. We see that P has the value of 60, which is 
n! _ 5! 

(n-R)! (5-3) !‘ 

This leads us to believe that the program is all right as it is 
written. If "playing computer" shows that incorrect results are obtained, 
then the flow chart and program must be corrected accordingly. If a 
computer were actually available, the next step would be to run the pro- 
gram on the computer to see if it would produce the correct results. 

If additional permutations were to be calculated, then new values 
would be assigned to N and R. 

The question might arise as to whether the program for the calculation 
N’ 

/ TT- p ^ t could be used to find the value of P , i.e., the case where R 
t JN-Kj i n n 

is equal to N. In this case, N-R would be zero. If we look at the 
steps obtained in "playing computer" for ^P^ and we come to line 110 
the first time and then compare K with N - R, we see that K is already 
greater than N-R. Consequently, they could never be equal. Thus, K 
would continue to become one greater; and no solution would be reached. 

The next question that might occur is whether or not the flow chart and 
the program might be modified to work for the case where "N = R." One 
possibility would be to change the question in the second decision box 
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to "K ^ (N-R)" instead of "K=(N-R)." Line 110 would then read "IF 
K ^ (N-R) THEN 140." Check to see if it would produce the desired result. 

Section III 

The program for might be written as follows: 

10 READ N,R 

20 LET Q=1 

30 LET K=2 

40 LET Q=Q*K 

50 IF K=N THEN 80 

60 LET K=K+1 

70 G0 T0 40 

80 LET T=1 

90 LET K=1 

100 LET T=T*K 

110 IF K=(N-R) THEN 140 

120 LET K=K+1 

130 G0 T0 100 

140 LET V=1 

150 LET K=1 

160 LET V=V*K 

170 IF K=R THEN 200 

180 LET K=K+1 

190 G0 T0 160 

200 LET C=Q/(T*V) 

210 PRINT N,R,C 
220 DATA 5,2 
230 END 

The flow chart and the program for are tentatively correct. 

This is the combination of n things taken r at a time. It can be seen 
that they are the same as the one for up to line 120. The remaining 
part is for the calculation of R; and of C = Q/(T*V). 

It is suggested for the student to "play computer" for the values 
of N = 5 and R = 2. 

It should be pointed out that flow charts and programs are not unique 

for any given problem. That is, there are probably a number of different 

programs which would produce the desired results. For example, in the 

N! 

program written for the procedure for calculating a fac- 

torial was used three times. Somewhere in the calculation of NI the 



computer would calculate an intermediate value which would be equal to 

(N-R)!. Also, the value of RI would be obtained during the calculation 

of NI. It might be possible to write a program which would select these 

N! 

intermediate values and store them for use in calculating 

It would be worthwhile for students who are interested in doing 
more with the computer to write a flow chart developing the idea in the 
preceding paragraph. 
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Flow Chart C 
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THE SOLUTION OF A TYPICAL PROBLEM CONCERNING MOTION 
OF OBJECTS IN A GRAVITATIONAL FIELD 
by Marion E. Hakewessel 

A ball is throvm into the air with an initial velocity V at an angle 
A with the horizontal at a height D above the ground. 

a) Plot the trajectory of the ball using k second intervals 
until it hits the ground. 

b) Determine the magnitude and direction of the velocity of 
the ball at the above k second intervals.* 

A problem such as the above can be easily programmed and solved by 
students. One must assume a preliminary knowledge of some computer 
language such as BASIC; only the simple READ, PRINT, LET, and IF ... THEN 
statements should be required. 

First the physics of the situation must be understood. Motions in a 

gravitational field are best analyzed if they are resolved into their 

horizontal and vertical com.ponents. The horizontal velocity remains 

constant, and acceleration is zero, their being no component of force in 

this direction. The vertical velocity increases under the constant 

2 

acceleration g = -9.8m/sec due to constant gravitational force. Thus we 
can set up two systems of equations: one for determining horizontal dis- 

placement and velocity; one for determining vertical displacement and 

I 

velocity. Then at selected time intervals we can ask the computer to 
print out the coordinates of the object's position with respect to a 
ground position and to add the velocities ^/■ectorially to get a resultant 
magnitude and direction. 

The suitable equations are as follows: (See Figure 1) Given 

initial velocity V in m/sec at angle A in degrees (positive or negative) 
with the horizontal at a height D in m from (above only) the ground 
surface. 

1) initial horizontal velocity: H = V*cos(A:7t/180) 

o 

2) initial vertical velocity: V = V*sin( A *71/180) 

o 

3) horizontal displacement after t sec: D^ = H *t 

2 

4) vertical displacement after t sec: D^ = D + " 9.8t /2 

5) horizontal velocity after t sec: = H^ 

6) vertical velocity after t sec: V„ = V - 9.8t 

2^2° 

7) speed at any instant: S = (V^ + V_ ) 

-1 

8) direct ion(from horizontal) at instant: N = tan ^^2^^1^ * 
Equations 3 to 6 are, of course, the familiar equations of motion with 
slightly different notation. Notice that the units must be consistent. 



*( Problem 13 and 16, page 341, 1960 edition PSSC text) 
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Notice, too, that the angles must be converted to radian measure before 
using the BASIC functions SIN and COS; also a conversion back to degree 
measure must be made if the BASIC arctangent function is used. 

The next step is to sketch a flow chart to help understand the logic 
of the program. Refer to. Figure 2. In the first box we read in the 
givens. Then we calculate the initial horizontal and vertical velocities. 

It -is at this point that we. want to "time" the flight; the computer will 
print all the variables at t=0, then recompute them a time increment k 
sec later. It will continue . this loop until the object "hits" ground’ 
(vertical displacement is less than zero). 

A program that follows the flow chart in BASIC language is in 
Figure 3. 

The results of a run using sample. data V = lOm/sec, A = 45 , D = 0 m, 
k = .1 sec, is shown above. Notice the following: When the .path begins 

to point down, the value of the angle from the horizontal begins to be 
negative. Also notice one negative value for vertical distance; the 
object apparently "hit" the ground between 1.4 and 1.5 sec. The conditional 
statement (line number 140) allows the computer one more calculation as 
long as the previous height above ground is not less than zero. Thus, if 
the values are plotted, the place where the object "hit" could be read 
from the graph. Students should also realize that the values for speed, 
are close to proportional to the lengths of the velocity vectors whose 
direction is given. The plotted points are closer together for slower 
speeds ... like the images on a "multiple flash" photograph. 

The program could be modified to include calculation of the exact' 
time of hitting the ground, the highest point on the trajectory, exact 
"range", etc. If it were desired to generalize the program to include 
cases where the initial position of the object was "below" the ground, an 
additional IF ...THEN loop would have to be included to determine whether- 
any positive values of D2 have yet been found. (There is always the 
possibility that the object never gets "up to" the ground surface.) Of 
course this program is. valid only "near" the earth’s surface where g is 
a constant. If the value of g as a function of height above the earth 
were substituted for -9.8, then the program could be further generalized, 

For example, define the following; 15 DEF FNG(D2) = .667E-11 * M/(D2*D2) 
where M is the mass of the earth and is the distance from tho center 
of the earth. Use FNG(D2) in place of a constant value of g. 
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Fig. 2 
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TRAJEC 16 : 13 



PX M0N 08/14/67 



m TIME 






H0RIZ D 


VERT D 


SPEED 


DIRECTION 


Bn ^ 






0 


0 


10. 


45. 


K| 1 .1 - 






.707105 


.658108 


9.3328 


40.7419 


m 






1.41421 


1.21822 


8.72485 


35.86 


• 3 






2.12132 


1.68032 


8,18936 


30.2944 


ifjj .4 






2. -82842 


2.04443 


7.74139 


24.0192 


Bill .5 






3.53553 


2.31054 


7.39685 


17.0684 


B 






4.24263 


2.47865 


7.17067 


9.56139 


Bn 






4.94974 


2.54876 


7.0742 


1.70985 


■01 * ® 






5.65684 


2.52086 


7.11274 


-6.20605 


B . 9 






6.36395 


2.39497 


7.28413 


-13.8924 


Bn 






7.07105 


2.17108 


7.57937 


-21.103 


BmI ^ ^ 






7.77816 


1.84919 


7.98473 


-27.6779 


1.2 






8.48527 


1.4293 


8.48444 


-33.5489 


R 1.3 






9.19237 


.911405 


9.06292 


-28.7193 


ill 






9.89948 


.295513 


9.70608 


-43.2376 


101 1.5 






10.6066 


-.418379 


10.4019 


-47.1732 


1 


3 












Bi 


10 


READ V,A,D,K 








1 


20 


PRINT "TIME","H0RIZ D","VERTD" 


, "SPEED", "DIRECTI0N" 






30 


LET 


H0= V*C0S(A*3. 1416/180) 






■ 


40 


LET 


VO = V*SIN(A*3. 


1416/180) 








50 


LET 


T = 0 








P 


60 


LET 


D1 = H0*T 








L 


70 


LET 


D2 = D + V0*T - 


9.8*T*T/2 






p 


80 


LET 


VI = HO 










90 


LET 


V2 = VO -< 9.8*T 








In 


100 


LET 


S = SQR(V1*V1 + 


< 

ro 

< 

ro 






1 


110 


LET 


N = 180*ATN(V2/V1)/3.1416 






ll 


120 


PRINT T,D1,D2,S,N 








p 

1 


130 


LET 


T = T + K 









140 IF D2> = 0 THEN 60 
150 PRINT 
160 G0T0 10 

170 DATA 10,45,0, .1,12.6,37,6.9, ,2 
180 END 
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THE USE OF THE COMPUTER IN 
FINDING THE APPROXIMATE VALUE OF PI 
by Harry A. Hurley 
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Time after time while teaching math, the writer had been asked by 
students, "How are we able to say that Pi has an approximate value of 
3. 14158470+?" or, "How do we know that the ratio of the circumference 
of a circle to its diameter is approximately 3.14158470+?" It was 
therefore decided to use the computer in an attempt to answer these 
questions. This was not a new computation, but by the use of the com- 
puter it was thought it could be demonstrated; and the task of computing 
each step could be shortened. 

By geometric proof, it could be shown that, given a regular in- 
scribed polygon of n sides of s length, the length of a side of a 

regular inscribed polygon of 2n sides is equal to / ~ ~t 

/ 2r^ - r/ 



4r^ - (s 

n 



Graphically, it was easy to show that as the number of the sides 
of the regular inscribed polygon was doubled, the closer the perimeter 
of the polygon came to the circumference of the circle. The above for- 
mula was used in computing the length of a side of each new polygon. A 
hexagon was used as the polygon to start with, and since 1 was used as 
the radius, 1 was the length of the side of the first polygon. 

To do this manipulating and computing by hand would have been a 
gargantuan job; the computer did it quickly. 

As stated above, to compute the value of Pi, a regular inscribed poly- 
gon of six sides was used as the first for our computations. The formula 
stated above was used to determine the length of the side of a regular poly- 
gon of 2n or 12 sides. The length of a side was multiplied by the number of 
sides to obtain the perimeter of the polygon. From this the ratio of the 
perimeter to diameter was computed. The following results were obtained: 



No. of 


Length of 


Perimeter 


Approximate 


Sides 


Side 




Value of Pi 


6 


1 


6 


3 


12 


0.517638 


6.21166 


3.10583 


24 


0.261052 


6.26526 


3.13263 


48 


0.130806 


6.2787 


3.13935 


96 


6.54382 E-2 


6.28207 


3.14103 


192 


3.27235 E-2 


6.28291 


3.14146 


384 


1.63624 E-2 


6.28315 


3.14158 


768 


8.18144 E-3 


6.28335 


3.14167 


The 


above table showed that when 


the number of 


sides increased beyond 


n X 2^ the 


approximation went over the 


value of Pi. 


The reason for this 



was that the computer rounded off to the nearest 8th decimal place. This 
caused an increase in the error which the machine could not correct. By 
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chopping manually, it was possible to carry the computations two steps fur- 
ther to an increase in the number of sides to 768, 1536 and, by a drastic 
chopping, 3072 and 6144. 



For a check to determine if there was an error created by the choice 
of the original number of sides, the same program was run where n = 4. The 
following results were obtained. The manual chopping was not perfon ed on 
this running of the program. 



No. of 


Length of 


Perimeter 


Approximate 


Sides 


Side 




Value of Pi 


8 


.765358 


6.12286 


3.06143 


16 


.390176 


6.24282 


3.12141 


32 


.196032 


6.27302 


3.13651 


64 


9.81342 E-2 


6.28059 


3.14029 


128 


4.90819 E-2 


6.28248 


3.14124 


256 


2.45428 E-2 


6.28296 


3.14148 


512 


1.22717 E-2 


6.28312 


3.14156 



This showed that there was no error in the number of sides chosen as a start, 
for the computations became too great at the same step. 

Given a side of a regular inscribed polygon of n sides, to find the 
side of a regular inscribed polygon of 2n sides, see Figure 1. 

Given : 0 0 with radius r, AB(s) a side of a regular inscribed polygon of 

n sides, and AC(x) a side of a regular inscribed polygon of 2n sides. 

To Find : the value of x. 

Solution: 



1 . 

2 . 

3. 

4. 



5. 



6 . 



7. 



8 . 



Statements 



Reasons 



Draw CO intersecting AB at D 1. 

^ 2 . 

AO = BO 3. 

CO is J_ bisector of AB 4. 




Two points determine a line. 

Def. of regular polygon 
Radii are = 

If two points are equidistant from the 
end points of a line segment in their 
plane and are on opposite sides of the 
line segment , then they determine the 
perpendicular bisector of the line segment. 
CO is X bisector of AD. 

Pythagorean theorem 



Subtraction 







I 



i 

] 

f. 

j 

! 

\ 

\ 
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The proof of the above formula was left to the students. To show 
how great a task the complete computation would be, the students were asked 
to carry out the computation of the approximation of Pi for a step or tvro. 

The flow chart for this computation is shown in Figure 2. 

Program for Computing Approximate Value of ^ 

10 READ S,N 

20 PRINT "N","S","C',"P" 

30 FOR K = 1 TO 10 

40 LET S = SQR(2-SQR(4-S*S)) 

50 LET N = 2*N 
60 LET C = S*N 
70 LET P = C/2 
80 PRINT N,S,C,P 

90 NEXT K 
100 DATA 1,6 
110 END 

So that the lopping could be done by hand, the following changes 
were made in the above program: 

10 INPUT S,N 
30 (TAKE 0UT) 

90 (TAKE 0UT) 

100 (TAKE 0UT) 
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Figure 1 
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CONSTRUCTING FLOW CHARTS AND ^miTING PROGRAMS FOR SIMPLE IRCBLEM? 

BY Dan H. Jarrell and Royce D. Jackson 

CONTENTS 

I . INTRODUCTION 

A . Computer language 

B. Time sharing systems 

C . Flow charts 

II. Construction of flow charts 

A . Symbols 

B. Loops and indexes 

C . Examples 

D. Exercises 

III. Writing computer programs 

A. Loops and print statements 

B. Writing BASIC programs from flow charts 

1 . Examples 

2. Discussion of approaches to improve programs 

C. Exercises 

IV. Expanding computer applications 

A. Formulas or mathematical models 

B. Examples 

1. Application in geometry relating circles 

2 . Simple series 

3. Game of twenty-three matches 

We are presently living in the early stages of the age of xha cicmputer. 
The public has a wide range of ideas concerning the computer. Man has 
devised languages in order that he may communicate with The ccmpi_tep<. 

Computer procedures for the solution of mathematical problems are usually 
written in an algorithmic language, a language that computer sciencistj 
have especially designed for expressing problem '-solving prccedur-s^ The 
most commonly used algorithmic languages are ALGOL (Algorithmic Language) 
and FORTRAN (Formula Translator). 

In this discussion we will describe a language called BASIC, which is 
a modified form of ALGOL. BASIC is designed especially for instruccicral 
use, it is simpler and easier to learn than ALGOL or FORTRAN. 

For educational purposes, a computer for each school is financially 
impossible at the present time. Thus, the Time Sharing System is presently 
being used for this purpose. This system consists of a large centrax 
computer that has many input-output stations, all of which can be in use 
at the same time. This means that a number of individuals can simiuiraneously 
use the computer on a "shared-time" basis , 

Each school is equipped with a teletype. The teletype is connected by 
long distance telephone lines to the computer. The teletype itself is 
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similar to a regular typewriter. 

Basically we use the time sharing system in the following manner: 

1. Write a flow chart and then write the problem solving procedure 
in BASIC language. 

2. T 3 rpe the program on the teletype without being in contact with 
the computer. The teletype will punch the program on tape. 

3. Then call the computer. The teletype is equipped with a system 
similar to a telephone. After getting in contact v;ith the computer, 
run the prepared tape through the teletype tape reader, storing the 
program in the computer. 

4. Direct the computer to run the procedure. The computer will run 
the procedure and type the results on the teletype. 

5. Debug the procedure, if necessary, and try again. 

Therefore, we must first of all learn to write problem procedures in 
the BASIC language and learn to use the teletype to communicate with the 
computer . 

In our discussion that follows the reader will gain some knowledge in 
using "Basic Language". However, it will be necessary to consult reference 
manuals for a more advanced study. 

Introduction' to Flow Charts 

Since it is evident that man- to -computer communication must be concise, 
consistent, and' complete to yield a satisfactory relationship, it also 
follows that the analysis and interpretation of a problem must have the 
same characteristics to facilitate its "translation into a workable program 
for the computer. A logical sequence of steps that can be understood and 
acted upon by a computer must be followed to gain a satisfactory solution. 
This sequence is obtained by a detailed analysis of the relationships 
found in a particular problem. Such an analysis may well be done by using 
a "flow chart" to indicate the required steps. A flow chart is a pictorial 
display indicating the chain of steps involved in the solution of a problem. 
Such a chart presents an overall view of a procedure that is frequently 
much easier to group than a description in words or in ordinary algebraic 
expressions. In addition to the sequence of steps involved, the chart 
includes conditional and repeating loops which give directive instructions 
at certain points in the sequence of steps which are vital to the 
efficiency of the total solution procedure. Usually the more complex 
problems require a more complex flow chart. 

This study in computer programming is primarily concerned with BASIC 
language, which adapts quite readily to program translation from flow 
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charts. However, flow charts can be readily adapted and utilized in 
conjunction with other languages and different models of computers. 

Symbols denoting specific operations or considerations in procedure have 
been standardized to a great degree to promote and facilitate universal 
acceptance and utilization of flow charts. This symbolic language is also 
understood by mathematicians v;ho may not be familiar with any computer 
language but who utilize their skills to analyze complex problems and 
forward the analyses by flow charts to a programmer . 

Once a flow chart is constructed, the major part of the total 
procedure in solving a problem by a computer is completed. A good flow 
chart gives sufficient details in information to write a workable computer 
program. No step that is necessary to a solution may be skipped or placed 
out of its proper sequence in a chart. Logic, proceeding step by step and 
in order, is vital to -successful problem solving. Mathematics is a 
science based' on organization and procedure. Construccion of flow charts 
and- writing computer programs can be most useful tools in teaching 
organization and procedure in problem solving and in making more concrete 
a student's understanding of the fundamental arithmetic operations, 
effective manipulations of terms and expressions, and recognition of 
relationships between quantities. 

Flow Chart' Symbols 

Symbols used in construction of flow charts are illustrated and 
explained below. A few S3^ibols are adequate to portray most problems that 
yield' themselves to computer application. 

1. The beginning of the program and points at which the computer is 
started are indicated by an oval. 

2. Points at which data values are either read into the computer or 
printed' out as the computer asks for information are indicated in and 
enclosed' by a pentagon resembling a rectangle with the top left-hand 
corner- clipped off. (It is given the familiar shape of the IBM card.) 

3. Statements instructing the computer to calculate a certain 
quantity or to set it equal to another quantity are enclosed in a 
rectangle . 

4. A point at which a decision is made and the program must choose 
between two or three branches to other parts of the program is indicated 
by conditional statements enclosed in a diamond. The appropriate branch 
is followed depending on the answer to the question or conditional 
restrictions in the decision box. 




5. If a flow chart is too large to fit on a single page or column, or 
if a branch line is difficult to draw between two points, a letter enclosed 
in a small circle serves as a connection link between sections. 

Loops and Indexes 

Loops in a flow chart are indicated by lines connecting parts out of 
their natural sequence, using arrows to denote the direction and to what 
point the computer should go to resume operation on the program. Repeating 
loops instruct the computer to repeat a sequence of steps that may include 
only a few steps or the entire program. An example of the latter is this: 
When one set of data is utilized in the program and an answer printed, it 
may be desired that a new set of data be entered and processed through the 
entire program. To accomplish this a loop is set up to direct the computer 
to the input statement after completion of the steps using the initial set 
of data. Conditional loops are used in conjunction with decision (diamond) 
boxes and may also be used with repeating loops. 

Indexes are used as counters or as subscripts to variables to enable 
counting- the number of passes through a particular loop and/or to 
differentiate between variables as they are read in or calculated . Loops 
and' indexes are illustrated in the examples listed following this section. 

In these examples of problems shown, certain descriptive or indicative 
words such as "let, calculate, add, input, and output" included in the boxes 
are normally left out in constructing flow charts. These are used to 
familiarize the student with the function of the symbols used. Succeeding 
examples and exercises will delete these words since the symbols used 
indicate the type of operation or statement. The examples illustrated do 
not necessarily indicate the most efficient procedures ; they are intended 
illustrate ways in which the symbols may be used in a program. 

Construct flow charts to represent the procedures in solving 
the following for Y, V, or A. 



Exercises: Constru 

the fol 


Al. 


y = 2a-b 


A2. 


A = ttR^ 


A3. 


V = 1/3-irR^h 


A4. 


3ax = 2d 


A5. 


l/3y - b = 5c 


A6. 


2A = (b^ + b^) 


A7. 


y. = 2B,^ 



2A = (bj^ + b^) h (set up a loop to repeat the program for more than 

one set of values for b^^, h^, & h) 

(set up loops to read 5 values of B and subscript 
both B and Y) 



A8. 



BY = A - C 
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(set up a conditional loop to avoid dividing by zero 
and another loop to use more than one set of data 
for A5B5 6 C) 

2 

A9. y = AX + BX + C (find points (x,y) to enable the drawing of the graph 

of the given function) 

2 

AlO. y. = 3C (write a statement to let C = C + 1 below the given 

^ statement, print out y., and set up a loop repeating 

the calculation of y. for 5 different values of C) 

■^1 

c = ab and values for a, b, 6c solve for d 
loops are included in this example) 

The oval indicates the starting point of the 
program . 



To solve for d, values for a,b, 6 c must be given. 
These are either read in from a "data" line or 
given as input information as requested by the 
computer . 



3. The desired quantity d is equated to the proper 
expression requiring rewriting the given equation 
in terms of d, and calculation is made. 




Example Al: Given: d - 

(no 



START 






INPUT 

A,B,C 



1 . 




4. The computer prints out the calculated value for d. 




5 . The sentence is completed so the program is 
terminated . 
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Example A2 : Given y = 1/X , solve for y 

A conditional loop and a repeating loop are used here. 




7. 



The program is started. 



Values for X are read from the data line or 
values are input as requested. 



The computer is instructed to compare the input 
value of X with zero . If it is equal to zero 
(the next step shows that this would mean 
dividing by zero, which is not permitted) a 
loop instructs the computer to return for 
another input value of X. If the comparison 
shows that X is not equal to zero then the 
decision box directs to the next step in the 
sequence . 

y is set equal to the indicated quantity and a 
calculation is made. 



The calculated value for y is printed. 



A loop is set up to direct the computer to 
the beginning of the program where a new 
value of X is taken and the entire program 
sequence is followed as before. 



When all values of X are read from the data 
line, or input values are terminated, the 
computer is instructed to terminate operations. 



o 
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Example A3: Given X 




solve for X for 



N values of A and of B 




Start the program. 



Again values are input as requested or are read in 
from data lines. N is used to limit the counter i. 



”i” is used to count the number of passes through the 
sequence and also as a subscript to XjA, S B 



”X” is set equal to the indicated quantity and 
calculated. It is also subscripted to differentiate 
between the various values calculated corresponding 
to the different input values of A and B. 



Here the values of A and B corresponding to the 
calculated value of X are printed as well as the 
X value. This enables correlating the calculated 
X value with the proper A and B values by inspection 
of the output since they are now grouped. 



This is not an ordinary statement of equality but a 
statement meaning that rhe value of i is to be changed 
by adding 1 to whatever value i held before reaching 
this point. If this is the first pass through the 
loop, then i held the value of 1 and so 1 is now added 
to that, giving i a new value of 2, Each pass through 
the loop at this point increases the value of i by 1. 



In the decision box i is compared to N , If i N , 
then a repeating loop is set up between the decision 



box and 



^A^ , and this loop is repeated until the 

value of i exceeds N. The program then proceeds to 
the next and last step. 



The program is now terminated. 



Loops and Print Statements in BASIC Programs 

In writing BASIC programs, if the flow chart contains loops, the 
program will also contain loops . There are several ways to set up the 
types of loops used. A conditional loop may use an "if-then" statement. 
"If" certain conditions are met, "then" the computer is directed to a 
specific point to resume operation on the program or to terminate it. 

A type of repeating loop is the "for" and "next" statements. The "for" 



statement introduces the loop and a "next" statement at the end of the 
loop directs the computer to the beginning of the loop for another value 
of the variable used as a counter and/or a subscript. The "for" statement 
sets the bounds on the variable. Passes are made through the loop until 
all values of the variable, up to and including the limiting value, have 
been used. Any number of these loops may be used in the same program. A 
loop may be entirely inside another such loop, but they must not intersect. 
The inner loop must be completed first before reaching the "next" 
statement that' is a part of the outside loop. Repeating loops, and also 
bypass loops, may be set up using a "go to" statement. In a repeating 
loop the "go to" statement directs the computer to a point in the program 
where a sequence of steps is repeated. In a by-pass loop the computer is 
instructed to by-pass certain steps and "go to" a certain point. 

The print statement in a program directs the computer to print out 
a message or quantity. If asked to print X, the value of X is printed 
and- not the letter "X". To print a message, symbols, or a letter as such, 
enclose them- in quotation marks. Following are examples of loops and 
print- statements in excerpts from hypothetical programs. 

Example Bl; 



5C FOR I = 1 TO N 
60 READA(I), B(I) 
70 NEXT I 



(A(k) is read: 
A sub I) 



"I" is used as a subscript, N as the upper limit of the 
subscript. A(I) and B(I) in the first pass through the 
loop will read two values from a data line for A(I) and 
B(I), store it in the computer memory, and then proceed 
to line 70. It is then directed to the next "I" in line 
50, which is now 2. Two values are again read for A(2) 
and B(2) and stored. This procedure continues until N 
values of A (to A(N)) and of B have been read. This 
terminates this loop and the computer proceeds to the 
next step in the program. 



Example B2 ; 

20 FOR I = 1 TO N 
30 FOR K = 1 TO M 
40 LET X(I,K) = 
A(I) * B(K) 

50 NEXT K 
60 NEXT I 



Here a loop is inside another loop. As the program 
progresses through these loops, in the first pass the 
first I, (I 1), is taken which is a subscript to both 
A and X. In the next step the first K,(K 1), is taken 
and these values are used to calculate the first X(I,K), 
which is X(l,l) A(1)*B(1). The next step directs to the 
next K, (K 2), in line 30. This time X(l,2) A(1)*B(2) is 
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calculated. To complete this phase M values of X are 
(XI, K) is read: calculated using only the A(l) value. The computer then 

X sub I,K) proceeds to the next step which directs to the next I, 

(I 2), in line 20. With this new value for I, the 
computer proceeds to the K Loop again and calculates M 
values of X again, this time using the A(2) value of each 
of the values for B from B(l) to B(M). This procedure 
continues until all values of I through N have been used. 



Example B3 : 

40 LET I = I T 1 Line 50 instructs the computer tc go to line 20 and 
50 GO TO 20 continue the program at that pcint. 



Example B4 : 

70 GO TO 90 

80 PRINT "NO SOLUTION" 
90 PRINT "FOR X=" X(K) 



Here line 70 instructs the computer to bypass line 
80 and "go to" line 90, where it prints: FOR X=4, 
(where we assume the value for X(K) to be 4), If 
executed, line 80 would print: NO SOLUTION 



Example B5 : 

70 IF K<N THEN 20 Line 70 instructs the computer tc compare the K 

80 PRINT’ "X("K"’)=" X(K) value at this point to the N value and make a 

decision. If K ^ N, then ix goes to line 20 and 
continues the program at that point. If K > N, 
then the computer proceeds to the next step 
(suppose K = 3 and X(3) = 12.6) and prints: 

X(3) = 12,6. 



0N AT 15:43 PX M0N 08/21/67 TTY 7 

USER NUMBER — P61000 

SYSTEM— BAS 

NEW 0R 0LD— NEW 

NEW PROBLEM NAME— P0LY1 

READY. 

TAPE 

READY. 

10 INPUT N 

20 LET S=(N-2)*180 

30 PRINT S 

40 END 

RUN 

WAIT. 

P0LY 15:44 PX M0N 08/21/67 

? 8 
1080 

TIME: 1 SECS. 

RUN 

P0LY 15J45 PX M0N. 08/21/67 

? 12 
1800 

TIME: 1 SECS. 

RUN 

P0LY 15:45 PX M0N 08/21/67 

? 24 
3960 

TIME': rSECS. 

BYE 

*** 0FF AT 15:46 PX M0N 08/21/67. 



EXAMPLE C 




EXAMPLE D 



0N AT 15:47 PX M0N 08/21/67 TTY 



USER NUMBER— P61000 

SYSTEM— BAS 

NEW 0R 0LD— NEW 

NEW PR0BLEM NAME— P0LY2 

READY. 

TAPE 

READY. 

10 -INPUT N 

20 LET S=(N-2)*180 

30 PRINT S 

40 G0 T0 10 

50 END 

KEY 

READY. 

RUN 



P0LY 15j48 PX M0N 08/21/67 

? 8 
1080 
? 12 
18Q0 
? 24 
3960 
? ST0P 



ran 3 SEC. 




ST0P. 
READY ^ 



BYE . 

*** 0EF AT 15:49 PX M0N 08/21/67. 





EXAMPLE 



0N AT 15:5C PX M0N 08/21/67 

USER NUMBER— P6 1000 
SYSTEM— BAS 
NEW 0R 0LD— NEW 
NEW PR0BLEM NAME— P0LY3 




ST0P. 

READY. 

BYE 

*** 0FF AT 15:53 PX M0N 08/21/67. 



O 
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ON 15:54 PX M0N 08/21/67 

USER NUMBER— P61000 

SYSTEM— BAS 

NEW 0R 0LD— NEW 

NEW PR0BLEM NAME— P0LY4 

READY. 

TAPE 

READY. 

10 READ N,Q 

20 IF N>=3 THEN 60 

30 PRINT "N0 S0LUTI0N" 

40 LET N=N+1 

50 G0 T0 20 

60 LET S=(N-2)*180 

70 PRINT S 

80 LET N=N+1 

90 IF N<=Q THEN 20 

100 DATA 0,-25 

110 END 

KEY 

READY. 

P0LY 15:56 PX M0N 08/21/67 
N0 S0LUTI0N 
N0 S0LUTI0N 
N0 S0LUTI0N 
180 
360 
540 
720 
900 
1080 
1260 
1440 
1620 
1800 
1980 
2160 
2340 
2520 
2700 
2880 
3060 
3240 
3420 
3600 
3780 
3960 
4140 

TIME: 0 SECS. 

BYE 



TTY 7 



EXAMPLE F 




0FF AT 15:56 PX M0N 08/21/67. 
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Exercises 

i'Jrite Basic Programs Using the 
Following Flow Charts 







sura 



Discussion of Approaches to Iraprove a Program 

All four examples were programs for the same problem, finding the 
of the interior angles of a polygon. In example C we* had to type "RUN" 
each time before we could input a new value of N. In this example, we used 
three polygons. We were on the computer three minutes and used three 
seconds of actual calculation time. : 

In example D our only addition to program A is statement 40, which 
sends the computer back to statement 10 for more input. We note that this 
time we do not have to type "RUN" each time we wish to input a new value of 
N. We used the computer two minutes and again used these records of 
actual calculation time. 

In example E we have a more efficient program simply because statement 
20 will not allow the operator to use a value of N less than three.. In 
examples D and E notice that the computer continued waiting for data until 
we typed "STOP", 

in example F notice that we have replaced statement ten with the term 
"READ". What is the purpose of this program? Here we desire to find the 
sum of the interior angles of each of several polygons beginning with a 
triangle and increasing the number of sides by one each time. N is the 

number of sides and Q is the maximum number of sides we wish to 
use. Again' note that if N is less than three we immediately increase N by 
one in statement' 40 and statement 50 sen'^s the computer back to statement 
• 20 with' a' new' value of N. When N is greater than or equal to three we 
then go' to* statement' 60 and compute'- the sum of the interior angles for 
the particular' value of N'. Notice in statement 70 that we increase N by 
one'. ' If- N' is' less than or equal to Q we return again to stateinent 20 and 
go' through the program again' finding the sum of the interior angles for 
our new "value of N. This loop will continue until N becomes greater than 
Q, then the program will be terminated. Notice that we used the computer 
two' minutes and had less than one second of acj-tual calculation time. Yet 
we- found- the' sum' of' the' interior angles for twenty-six different values of 

in the' foregoing examples and exercises in flow char-ts and in BASIC 
programs, ■ observe that a mathematical model or formula is .set up initially 
to' represent' a' problem or- a. desired quantity to be calculated. Using this 
formula,' or' formulas, a procedure (flow chart) is devised to obtain the 
solution.' ' -In' simpler problems these formulas are quite obvious. In more 
comp'lex* problems it is sometimes difficult to set up a formula , that will 
generate the desired quantity in .a manner that lends. itself to computer 





application. Often a series of terms are involved in a generation. To 
reserve a separate part of the program for each term desired would be 
prohibitive in time and space required for writing such a program. Use of 
loops to generate such quantities is much more efficient but requires a 
certain amount of insight and skill in mathematical procedures. To develop 
skill in this area requires much practice and application. Expanding 
computer application to include more complex geometric problems and those 
related to sequences are illustrated here in a limited sense to give an 
insight into the possibilities that exist, or may exist through future 
research, in the field of computer applications. In-depth study may be 
pursued as the teacher or student chooses. On completion of this 



elementary study presented here, the student is encouraged to continue in 
individual investigations and studies. Following are examples illustrating 



applications in geometry, sequences, and even in recreation. 



CIR 


16:05 PX TUE 


08/22/67 


TANGENT 


EXTERNALLY 


3 


9 


4 


5 


2 


2 


5 


10 


7 


2.23607 


10 


8 


2.23607 


CIRCLES 


ARE C0NCENTRIC 




CIRCLES 


ARE THE SAME 




0UT 0F DATA IN 10 




TIME: 


1 SECS. 





2 5 

TANGENT INTERNALLY 
EXTERNAL N0N INTERSECTING 
INTERNAL N0NINTERSECTING 
INTERSECT 



LISTNH 

10 READ R1,R2,X1,Y1,X2,Y2 

20 IF RKR2 THEN 160 

30 LET- D=SQR((X1-X2)i2+(Y1-Y2)i2) 

40 IF D=0 THEN 110 

50 IF R1+R2=D THEN200 

60 IF R1-R2=D THEN 220 

70 IF R1+R2<D THEN 240 

80 IF R1-R2>D THEN 260 

90 IF R1-R2<D THEN 280 

100 G0 T0 10 

110 IF R1=R2 THEN 140 

120 PRINT "CIRCLES ARE C0NCENTRIC" 

130 G0 T0 10 

140 PRINT "CIRCLES ARE THE SAME" 

150 G0 T0 10 
160 LET T=R1 
170 LET R1=R2 
180 LET R2=T 
190 G0 T0 50 

200 PRINT "TANGENT EXTERNALLY", Rl, R2 ,D 
210 G0 T0 10 

220 PRINT R1,R2,D, "TANGENT INTERNALLY" 

230 G0 T0 10 

240 PRINT R1,R2,D, "EXTERNAL N0NINTERSECTING" 
250 G0 T0 10 

260 PRINT R1,R2,D," INTERNAL N0NINTERSECTING" 
265 G0 T0 10 
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i 



% 



iD 



.1 






I O 

ERIC 



280 PRINT R1,R2,D, "INTERSECT*’ 

290 G0 T0 10 

300 DATA 3, 2, 2, 3, -1,-1, 9, 4, 2, 3, -1,-1, 2, 2, 2, 2 
305 DATA -1,-2, 10, 7, 1,1, -1,0, 10, 8, 1,1 
310 DATA -1,0, 6, 5, 5, 4, 5, 4, 6, 6, 5, 4, 5, 4 
670 END 



3 



'i 



I 



NEW 

NEW PR0BLEM NAME— SUMIOO 
READY. 

TAPE 

READY. 

10 PRINT"SUM 0F THE FIRST 100 INTEGERS" 

20 READ N 

30 LET K=1 

40 LET S=0 

50 LETS=S+K 

60 LET K=K+1 

70 IF K<=N THEN 50 

80 PRINT"S="S 

90 DATA 100 

100 END 

KEY 

READY. 

RUN 

SUM 100 8:42 PX THU 08/24/67 

SUM 0F THE FIRST 100 INTEGERS 
S= 5050 

TIME: 0 SECS. 

BYE 

*** 0FF at 8:42 PX THU 08/24/67. 



SQRRTS 8:39 PX THU 08/24/67 
SEQUENCE 0F SQUARE R00TS 



y 



X 



SQR[X] 

1 



I 



2 


1.41421 


3 


1.73205 


4 


2 


5 


2.23607 


6 


2.44949 


7 


2.64575 


8 


2.82843 


9 


3 


10 


3.16228 


11 


3.31662 


12 


3.4641 


13 


3 . 60555 


14 


3.74166 


15 


3.87298 


16 


4 


17 


4.12311 


18 


4.24264 


19 


4.3589 


20 


4.47214 


21 


4.58258 


22 


4.69042 


23 


4.79583 


24 


4.89898 


25 


5 


TIME: 0 


SECS. 


LISTNH 




5 PRINT" 


SEQUENCE 0F SQL 


10 READ 


N 


20 PRINT 


"X","SQR(X)" 



i 



3 



I 



3 



•J 



s 






S 



9 



30 F0R X=1 T0 N - 
40 PRINT X,SQR(X) 

50 NEXT X 
60 DATA 25 
70 End 

NOTE: LISTNH and the program printed 

below the .solutions serve as 
a means of making corrections 
in the original program and 
having the computer tjrpe out 
a corrected program. LISTNH 
is typed to instruct the 
computer to accomplish this. 







Installments where Monthly Pa3nnents Pay Interest Due for the Previous 
Month and the Balance of the Payment is Applied to the Principal: 
$ 100006 % 

PAYM 10:08 PX THU 08/24/67 

5 PRINT"N0.", "PAYMENT"," INTEREST", "0N PRINCIPAL" , "BALANCE 0N P" 

10 READ R,D,A 

20 LET K=1 

30 LET J=1 

40 LET B=R*D/12 

50 LET C=A-B 

60 LET D=D-C 

70 PRINT K,A,B,C,D 

80 LET K-K+1 

90 IF J>=2 THEN 160 

100 IF D>=A+B THEN 40 

no LET C=D 

120 LET A=C+B 

130 LET J=J+1 

135 LET D=0 

140 G0 T0 70 

150 DATA .06, 1000,60 

160 END 

RUN 



PAYM 


10:08 PX 


THU 08/24/67 






N0. 


PAYMENT 


INTEREST 


0N PRINCIPAL 


BALANCE 1 


]. 


60 


5. 


55 


945 


2 


60 


4,725 


55.275 


889.725 


3 


60 


4.44863 


55.5514 


834.174 


4 


60 


4,17087 


55.8291 


778.344 


5 


60 


3.89172 


56.1083 


722.236 


6 


60 


3.61118 


56.3888 


665.847 


7 


60 


3.32924 


56.6708 


609.177 


8 


60 


3.04588 


56.9541 


552.223 


9 


60 


2.76111 


57.2389 


494.984 


10 


60 


2.47492 


57.5251 


437.459 


11. 


60 


2.18729 


57.8127 


379.646 


12 


60 


1.89823 


58.1018 


321.544 


13 


60 


1.60772 


58.3923 


263.152 


14 


60 


1.31576 


58.6842 


204.468 


15 


60 


1.02234 


58.9777 


145.49 


16 


60 


.727449 


59.2726 


86.2173 


17 


60 


.431087 


59.5689 


26.6484 


18 


27.0795 


.431087 


26.6484 


0 



TIME: 1 SECS. 

Game of 23 Matches 

In the 23 match game there are two players. Each in turn must 
pick up 1, 2, or 3 matches. The player who is forced to pick up the 
last match loses. 
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in 



gi\ 



k 





23GAME 12:56 PX WED 08/23/67 
10 PRINT"23 MATCH GAME:K0M PEWTER VS.Y0U" 

20 LET N=23 
30 LET N=N-2 

40 PRINT"I,K0M PEWTER, TAKE 2 MATCHES" 

50 PRINT"THERE ARE"N"LEFT.H0W MANY D0 Y0U TAKE?" 

60 LET 1=1 
70 INPUT K 
80 IF 4-K>3 THEN 160 
90 IF 4-K<l THEN 160 
100 LET N=N-K 
110 IF N=0 THEN 240 
120 LET N=N-(4-30 

130 PRINT" I , K0MPEWTER ,TAKE"4-K"MATCHES . " 

140 PRINT"THERE ARE"N"MATCHES LEFT.H0W MANY D0 Y0U TAKE?" 

150 G0 T0 70 
160 LET 1=1+1 
170 IF I<3 THEN 212 
180 IF I >3 THEN 220 

190 PRINT"THAT'S TWICE Y0U’VE CHEATED-0NCE M0RE AND I QUIT’] INPUT K 
200 PRINT"N0W, CH00SE EITHER 1,2,0R 3 MATCHES" 

210 G0 T0 70 

212 PRINT" Y0U CHEATED. CH00SE 0NLY 1,2,0R 3 MATCHES" 

215 G0 T0 70 
220 PRINT"! QUIT" 

230 G0 T0 250 

240 PRINT"H0W AB0UT THAT, I W0N AGAIN. WH0*S NEXT?" 

250 END 
RUN 

23GAME 12:59 PX WED 08/23/67 
23 MATCH GAME:K0M PEWTER VS.Y0U 
I,K0M PEWTER, TAKE 2 MATCHES 
THERE ARE 21 LEFT.H0W MANY D0 Y0U TAKE? 

? 1 

I,K0MPEWTER,TAKE 3 MATCHES. 

THERE ARE 17 MATCHES LEFT.H0W MANY D0 Y0U TAKE? 

? 3 

I,K0MPEWTER,TAKE 1 MATCHES. 

THERE ARE 13 MATCHES LEFT.H0W MANY D0 Y0U TAKE? 

? 5 

Y0U CHEATED. CH00SE 0NLY 1,2,0R 3 MATHCES 
? 2 

I,K0MPEWTER,TAKE 2 MATCHES. 

THERE ARE 9 MATCHES LEFT.H0W MANY D0 Y0U TAKE? 

? .5 

THAT'S TWICE Y0U'VE CHEATED-0NCE M0RE AND I QUIT 
N0W, CH00SE EITHER 1,2,0R 3 MATCHES 
? 2 

I,K0MPEWTER,TAKE 2 MATCHES. 

THERE ARE 5 MATCHES LEFT.H0W MANY D0 Y0U TAKE? 

? 3 

I,K0MPEWTER,TAKE 1 MATCHES. 

THERE ARE 1 MATCHES LEFT.H0W MANY D0 Y0U TAKE? 

? 1 

H0W AB0UT THAT, I W0N AGAIN. WH0'S NEXT? 

TIME: 13 SECS. 




LET N=N-(4-K) 
STOP 
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THE INTERSECTION OF CIRCLES 
by Nick R. Kalianov & Pearl Price 

Our purpose is "to investigate, using analytic methods, the problem of 
the intersection of two circles lying in the same plane. 

Consider two circles wirh centers and B(x 2 ,y 2 ) radii r^^ 

and r^ respectively. Let d = m(AB). Assume r^^ 

In Figure 1, a coordinate system is set up so that A is the origin. 



AB is the x-axis; and the x-coordinate of B is nonnegative. The equations of 

2 2 2 2 2 2 
the circles are: (l)x + y ~ ^2. (x-d) + y “ ^ 2 * 

For a Doint P to belong to both circles, the coordinates of P must 

tr 2 2 

satisfy both (1) and (2). If we subtract (2) from (1), we get (3) (x + y ) - 
[(x - d)^ + y^] “ " ^2* ^ point P belongs to both circles if and 

only if its coordinates satisfy both (1) and (3). Equation (3) can be 
simplified algebraically to (3’) 2dx = r^ - r^ + d . Show this simplification 
We can conclude from (3’) that if d = 0 there are points common to both 
circles only if r^^ = r^, i.e., the circles are the same. Also, if d = 0 and 
r>^ f r^j then the circles do not intersect. 

Let's assume that d ^ o. Then P belongs to boxh circles if and only if 



(i) x^ + y^ = r^ and (2’) x = 



2 2 ^ ,2 
+ d 

2d 



Note that (2’) is the equation 



2 2 ^ ,2 
+ d 



2d 



0 . 



of a li.n6 p6]?p6ndicul3.r* *to *tli6 x axis • Sinca r*^ ^ 

Hence we can reduce the problem to the intersection of a circle (with the 
radius r^^) and a line (where the distance between the line and the center of 



the circle is 



2 2 ^ ,2 

■^1 - ^^2 ^ 

2d 



). 



We may now conclude that, if d 0, the two circles intersect in two 



points if and only if 



2 2 ^ ,2 

- r.2 + d 



< r. 



Simplifying the inequality, we 



2d ' ‘1 

find that the two circles intersect in two points if and only if d< r^^ + r^ 
and d ■> r^^ - r^. 

The two circles will intersect in one point if and only if 

2 2 ^ ,2 
r» — r + Q 

1 2d ” ^1‘ Simplifying the equality, we find that the two circles 

will intersect in one point if and only if d = r^^ - v ^ ov d - r^^ + • 

The intersection of the two circles will be the empty set if and only 

2 2 ^ ,2 

r — r + d 
1 2 

if ^ >^1 



Simplifying the inequality, we find that the two 



i 



Ceric 









'1 
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circles do not intersect if and only if d < r^^ - r^ or_ d > r^^ + r^ . 

We have now completed our investigation. See Figure 2 for a summary 
of our analysis . 

Figure 3 is a flow chart that graphically depicts an overall view of 
the above procedure. The flow chart must be logically complete and must 
lead step by step to the desired answer. Do you see that it does? 

As an exercise for the student, select four pairs of circles, noting 
the coordinates of the center and the radius of each circle. For each pair 
of circles, follow the flow chart through until you reach a print 
statement. Check the result of the flow chart graphically. 

We have written out the program in BASIC directly from the flow 
chart. Study this program, using the flow chart as a guide, to see that it 
is correct. The program includes data for five pairs of circles that are 
externally tangent, are the same, intersect in two points, do not intersect, 
and do not intersect, respectively. 

As an exercise for -the student, run this program with data statements 
of your own. 

Exercises 

1. In our discussion of the intersections of circles, we made use of 
previous knowledge concerning the intersection of a circle and a 
line. Construct a flow chart, and write a program that will 
determine whether or not: 

a. the circle and line intersect in two points, 

b. the circle and line intersect in one point 

c. the circle and line do not intersect. 

2. Using the program above, describe the relative position of two 
circles of radii 5 and 9 if the center line has the length 

(a) 14 (b) 4 (c) 12 (d) 16. 

PROGRAM 

10 Read x^^, y^^, x^, r^^, r^ 

20 If r^^> = r^ then 50 

30 Print " Pick r^^> = r^” 

40 Go to 10 

50 Let D = SQR( (x^ - + (y^ - y^^)^ ) 

60 If = x^ then 120 
70 If D = r^^ + r^ then 190 



80 If D = - r >2 then 210 

90 If D > + T 2 then 230 

.100 If D > r, ** ^2 

110 G?5 tj5 150 

120 If then 140 

130 tj5 70 

140 if then 170 

150 Print "Circles are Concentric" 

160 G^ t^ iO 

170 Print "The circles are the same" 

180’ G(z5 10 

190 i^inj; "The circles are externally tangent" 
200 G^5 t;$ 10 

210 . I^int "The circles are internally tangent" 

220 X.i> 10 

230 R?iht "One circle lies outside the other" 

240 G^5 tjJ 10 

250 Print. "The. dircles intersect in two points*-' 
.260 G?5. tjz5 10 

.270 Data- .0 ,.0^j0 j 2 ,1 jl; ,0,0 ,0,0, 4,4 ' 

.280 Data: 2i0^0-,2,3,2, 0,0,0, 1^4,1 





300 .End 
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‘ Fig. 3 - Continued 
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FEEDBACK AROUND AN INTEGRATOR 
by Dean C. Larsen 
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This unit presents an additional section for use with the "Feedback" 
chapter (C-1) in the Engineering Concepts Curriculum Project text. The 
unit contains a presenration of a computer model of feedback at a level 
which would be understandable to the average, college bound high school 
student. There is a short discussion of why the material is included in 
the text and how it is presented. Following this discussion, there is a 
detailed description of how one might build a model of a dynamic feedback 
system in the computer. It is hoped that this model will help students 
understand the concept of feedback. It is assconed that each student has 
some exposure to BASIC computer language and that section B of the text 
has been covered. Section B deals wirh the model of dynamic systems, 
and it is assumed that students recognize the significance of the area 
under a rate of change curve. 

During the last ten years , engineers have begun to recognize feed- 
back as the underlying principle behind any system in which there is goal 
seeking behavior. It is this universal applicability of the feedback 
concept to both physical and nonphysical systems which requires its 
presentation in this text. The essence of feedback is that the output 
of ^ system is used to modify the input into the system so that the output 
is changed to conform to some predetermined goal . Many qualitative exam- 
ples of feedback are introduced at the first of the chapter to help 
students develop their understanding of this simple but elusive concept. 

In thinking about feedback, consider the system to be a student 
and all elements relating to his progress in a difficult course. The 
output of the system is the set of grades which the student earns. The 
input is all the study and work which the student does. How does the 
concept of feedback relate to this system? 

Assume that the goal for the system is a cumulative grade average 
of A in the course for this student. During the first week of school 
the student goes to class , takes a few notes and in general tries to 
understand the material. At the end of the week the teacher gives a 
quiz. The student gets some grade. This is the first output. Now 
feedback begins to take over and should direct the system toward the goal. 
The student makes a comparison between the grade (output) and the kind 
of study he is doing (input). This comparison causes the student to 
change his input in a direction which will force the output toward the 
desired goal. If the quiz grade is low, the feedback causes the student 






to adjust his rate of study (input) so as to bring the level of output up. 
If the grade (output) is higher, less correction in the rate of study is 
needed. Note that in this qualitative example of feedback, output changes 
input to produce changed output. This is feedback. 

Before leaving this example, note how the quality of the feedback 
control would be improved if the student were aware of the "SUM" of all 
of his grades in the course as opposed to just being aware of the "last" 
grade earned. Consider the student's situation if all grades were F 
except the last one, which was A. If the feedback does not keep track 
of the "SUM" of all grades earned, the student might reduce his rate of 
study on the basis of the last A and yet be on the verge of failing the 
class and missing the goal. Remember this idea of "SUM." The student 
would be better off keeping track of the sum of his grades rather than 
just his most recent grade. 

Many examples of feedback producing goal seeking behavior are used 
in the text. Feedback in physical systems is discussed, and students 
should have a feel at this point for the presence of feedback in a machine 
directing the system toward a goal. Following this treatment, several 
quantitative examples of feedback around a scalor or multiplier are 
presented. The following example shows the nature of this quantitative 
description. 

Consider the system WITHOUT feedback (Figure 1). The desired out- 
put is 5/2 times the. input at the SAME INSTANT. The output is nothing 
but a scaled up picture of the input at the same instant. It is obvious 
that in this system, output is very sensitive to changes in system 
characteristics. For example, if the value of the scalor A is reduced 
by 50% to 5/4, the output also drops by 50%. 

Consider the system WITH feedback (Figure 2). The output of this 
system will be the same as the output of the previous system. The feed- 
back loop will cause this system to be less sensitive to changes in system 
characteristics. It is important to realize that in this system, the 
input at any instant is a function of the output at that same instant. 

That is, output modifies input at the same instant. Hence, we can sayt 

y = (x - 0.4y)(A) 

This means that the output y is a function of the input x which is 
modified by a factor of the output y at the same instant. Both of 
the y's in our equation represent output at the same instant, so we 
can solve the equation for y. 
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y = (x - 0.4y)(A) 
y + 0.4yA = xA 
y(l + 0.4A) = xA 





Multiply numerator and denominator of the fraction on the right by 1/A. 



Since the output of this system is to be the same as the output of the 
first system, A is made very large. Then 



The text spends considerable time discussing examples such as this. 
The fact is stressed that the output is relatively independent of large 
changes in the scalor A. 



A 50% drop in the value of the multiplier A causes very little change in 
the output y. 

The preceding discussion is to inform the reader in a general way 
about the method of presentation in the chapter. There is a problem with 
this presentation in that most feedback controlled systems are not simple 
scalors. The output is not a simple multiple of the input at the same 
instant. The fact is that most controlled systems act more like inte- 
grators than scalors. The following material is written to give stu- 
dents insight into the nature of feedback around an integrator. It also 
provides an interesting example of computer modeling. The material is 
written as it might appear in the text. 

The careful reader will observe a difference between the qualitative 
examples of feedback (the car on page C-1.2, the ship on page C-1.9) and 
the latter quantitative examples starting on page C-1.14. The systems 
which are treated with algebra involve feedback around a scalor. Before 
reading further, see if you can determine how these systems differ from 
the qualitative examples at the first of the chapter. 




y = (5/2)x 



because is very small. 



Let A = 1000 



Let A = 2000 
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The truth of the matter is that many of the descriptive examples 
are not true scalors. These systems produce an output which is not 
a simple scaled up replica of the input. What one must be aware of at 
this point is that in many real systems , the output is determined by all 
previous values of the input and not by the most recent value. Consider 
the following as an example of this idea. 

Suppose elevator number 13 is dropping at full velocity past the 
25th floor of the "Mile Hi" building at this very instant. The operator 
is inputing a velocity signal with the control handle, and the elevator 
is producing a change in location. The operator senses how far the 
elevator has moved by viewing the floor numbers as they pass by a little 
window. 

25-24-23-22-21-20-19-18-17-16- 

The operator remembers the system goal. 

"Stop this car at the ground floor after traveling down 25 stories." 
Quickly the operator makes a comparison between output (9 stories gone) 
and the input signal. This input signal is really the position of the 
control handle. This "feedback comparison" causes the operator to change 
the control handle position (input). The elevator slows but does not stop. 
It continues to drop at a reduced velocity. 

-r--15 14 13 12 11 10 9 8 7 6 
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Again the operator senses the output (19 stories gone). A quick compari- 
son of this output with the input (control handle position) again causes 
him to change the input. The elevator continues to drop but at an even 
lower velocity. 

Eventually this feedback procedure will produce the sought after goal of 
25 stories displacement in a downward direction. Notice that the output 
of the .elevator (distance traveled) is a result of all previous input 
velocity signals prior to the time the elevator stops. This elevator 
is really an integrator. 

The output of an integrator at any instant is determined by all 
previous values of the input up to that instant. Since most of the sysr 
terns we control in the real world are dynamic, that is they change with 

f 

time, we can see that they. are really integrators and not simple scalors. 

Let us now develop a systematic example of feedback around an integrator. 

V = velocity in ft. /sec. — this is the input signal. 

D = displacement in ft. — this is the output signal. 

F = feedback signal 
E = corrected velocity signal 
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Figure 3 represents feedback around an integrator. The variables 
are defined above. Let the sysrem represent an automobile. The inte- 
grator will total or "SUM" up rhe distance traveled by the automobile 
after some arbitrary smarting Time. The feedback loop senses the total 
distance traveled by the car and corrects or adjusts the velocity accor- 
dingly, This feedback signal indicares how much the original velocity 
should have decreased at the end of any time interval. As rhe total 
distance D increases, the feedback loop reduces the velocity to a lower 
and lower value. The sealer constant -4/10 was picked because of the 
behavioral goal which we have chosen for the system. The method used to 
pick this constant is beyond the scope of this book. However, methods 
are available for picking constants in light of chosen goals. Given an 
initial velocity of 88 ft. /sec., the goal for this system is to slow the 
car so that the velocity is zero when the car has traveled 220 feet. 
Assume the following values for the variables at time T = 0 sec. 

V = 88 ft. /sec, — Input 

T = 0 sec. — Time 

D = 0 ft. The car has not yet gone any distance. 

Now, from Figure 3, we can see: 

F = -0.4(D) 

F = 0, no feedback signal 

E = V + F 

E = 88 ft. /sec. This is the corrected velocity signal. 

Now, assume that E holds constant for one second (see Figure 4). At 
the end of one second, what is the value of D? D is represented by the 
area under the corrected velocity curve (see Figure 5). 

D = 1 sec. X 88 ft. /sec. 

D = 88 ft. 

During the first second the car traveled 88 ft. The new values of the 
variables after one second are : 

T = 1 sec. 

D = 88 ft. 

F = -0.4 X 88 

F = -35.2 ft. /sec. 

Now, let us see what happens if E x-emains constant at this new value 
for one second (see Figure 6). You notice that the total area under 
the curve has increased during the time interval T = 1 sec. to T = 2 
sec. However, the increase in area during this time interval is less 
than during the Time interval T = 0 sec. to T = 1 sec. This is because 
the car travels at a lower corrected velocity of 52.8 ft. /sec. during 
the interval T = 1 sec. to T = 2 see. 

Increase in area under the curve = 1 sec, x 52.8 ft, /sec, 

= 52.8 ft. 






X 


SCALOR 


y 


input 


(10/4) 




output 



y = (lQ/4)x 



Figure 1« Without Feedback 




y = (x - 0.4y)A 



Figure 2. With Feedback 
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Figure 3. Integrator with Feedback 
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Increase in distance is then = 52.8 ft. 

Total Distance = 88 ft. + 52.8 ft. 

= 140.8 ft. 

The new values of the variables after 2 seconds are: 

T = 2 sec. 

D = 140.8 ft. 

F = -0.4 X 140.8 

F = -56.32 ft. /sec. 

I • . I 

This feedback signal is added to the original velocity signal to get a 
new value of E, the corrected velocity signal. 

E = V + F 

E = 88 - 56.32 

E = 31.68 ft. /sec. 

Observe what happens if E remains constant at this new value for 1 sec. 

(see Figure 7). 

Increase in Distance = 1 sec. x 31.68 ft. /sec. 

= 31.68 ft. 

Total Distance = 88 ft. + 52.8 ft. + 31.68 ft. 

= 172.48 

Figure 7 shows that the corrected velocity signal E gets smaller. This 
means that the total distance traveled will approach some limit as the cor- 
rected velocity approaches zero. This limit on total distance should be 220 
ft. in accord with the system goal. How can this be determined? One way 
would be to make calculations by hand like those above for many 1 sec. inter- 
vals and check to see if D approaches 200 ft. Since these calculations are 
very tedious , it would be helpful to write a computer program in BASIC and 
let the computer become a model of the d3mamic feedback system. By studying 
this model, it can be determined . if the system is seeking the goal. 

Study the flow chart (Figure 8) to make sure you understand the repete- 
tive nature, of the calculations vjhich the computer will make. If you do not 
understand the flow chart, go back and study the. hand calculations shown above. 

A computer progr^ and subsequent run for the system are shown in 
Figure 9. Notice how the system "seeks" the goal. Velocity is corrected 
until it approaches zero. Distance traveled increases and approaches 220 
ft. ihis is feedback in a dynamic system. Output is measured and used to 
alter input which produces new output. The system displays goal seeking 
behavior. 

This model could now be compared to the real world and improved so 
that it more accurately represented the physical system. For example, we 
could use a time interval T2. of 0.00001 sec.. This would have the effect of 
smoothing out the curves in Figures 4-7. This changed model would more 
accurately represent the real world, as automobiles cannot instantaneously 
change from a velocity of 88 ft. /sec. to 52.8 ft. /sec. Other changes in 
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the model could be introduced to improve the representation of the 
real world. In any case, we can new see the importance of feedback 
in controlling systems. 
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FDBMDL 15:27 PX M0N 07/17/67 

100 REM D = INITIAL VALUE 0F DISPLACEMENT WHEN SYSTEM IS STARTED 

110 REM B = FEEDBACK L00P SCAL0R 0R MULTIPLIER VALUE 

120 REM V = INITIAL VALUE 0F VEL0CITY AT TIME T 

130 REM E = C0RRECTED VEL0CITY SIGNAL 

140 REM F = FEEDBACK SIGNAL 

150 REM D1 = CHANGE IN DISPLACEMENT 

160 REM T = INITIAL VALUE 0F TIME WHEN SYSTEM IS SET IN M0TI0N 

170 REM T1 = T0TAL TIME SYSTEM IS M0DELED 

180 REM T2 = TIME INTERVAL SIZE DURING WHICH E IS C0NSTANT 
190 RlAD T,T2,D,B,V.T1 



200 DATA 0,1,0,-0.4,88,20 
210 PRINT "t","V",'’E","D","F" 

220 PRINT "SEC", "FT/SEC", "FT/SEC", "FT", 

230 LET F=B*D 

240 LET E=V+F 

250 PRINT T,V,E,D,F, 

260 LET D1-E*T2 
270 LET T=T+T2 
280 LET D=D+D1 
290 IF T<=T1 THEN 230 
300 END 

RUN 



FDBMDL 


15:28 


PX M0N 07/17/67 


T 


V 


E 


SEC 


FT/SEC 


FT/SEC 


0 


88 


88 


1 


88 


52.8 


2 


88 


31.68 


3 


88 


19.008 


4 


88 


11.4048 


5 


88 


6.84288 


6 


88 


4.10573 


7 


88 


2.46344 


8 


88 


1.47806 


9 


88 


.886837 


10 


88 


.532102 


11 


88 


.319261 


12 


88 


.191557 


13 


88 


.114934 


14 


88 


6.89605 E-2 


15 


88 


4.13764 E-2 


16 


88 


2.48258 E-2 


17 


88 


1.48954 E-2 


18 


88 


8.93736 E-3 


19 


88 


5.36239 E-3 


20 


88 


3.21746 E-3 


TIME: 


1 SECS. 





FT/SEC" 



D 


F 


FT 


FT/ SEC 


0 


0 


88 


-35.2 


140.8 


-56.32 


172.48 


-68.992 


191.488 


-76.5952 


202.893 


-81.1571 


209.736 


-83.8943 


213.841 


-85.5366 


216.305 


-86.5219 


217.783 


-87.1132 


218.67 


-87.4679 


219.202 


-87.6807 


219.521 


-87.B084 


219.713 


-87.8851 


219.828 


-87.931 


219.897 


-87.9586 


219.938 


-87.9752 


219.963 


-87.9851 


219.978 


-87.9911 


219.987 


-87.. 9946 


219.992 


-87.9968 



Figure 9 
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A SHORT DESCRIPTION OF A NARRATED SERIES OF SLIDES 
WHICH PRESENT TWO INTRODUCTORY COMPUTER LESSONS 
FOR JUNIOR HIGH SCHOOL 

Narration and Slides prepared by: 

Donald Marcotte 

K. Wayne Mahan 

Dan Colvin 

Photographic assistance and art work provided by: 

Jefferson County Instructional Aids Department 

The Jefferson County Schools have two small instructional 
computers (Bitran Six) which are used in eighteen junior high schools 
to illustrate what a computer is and how it functions. The computers 
are moved from school to school and also provide motivation for 
study of numeration systems. 

As their project for the NSF Institute Mr. Marcotte, Mr. Mahan, 
and Mr. Colvin prepared two complete lessons with tape recorded nar- 
ration and color slides to present to junior high classes. 

The first lesson provides a brief introduction to computers, 
a specific introduction to a ’’teaching computer” (The Bitran Six), a 
presentation of ' the’ numeration system used' by computers, and simple 
programming' procedures for’ the Bitran Six. 

I 

The second lesson* presents programs for addition, subtraction, 
multiplication, and division, and suggests other programs for student 
work . I 

Copies of tapes and slides or printed copies of the narration 
with accompanying sketches of the slides may be obtained at cost from 
the Jefferson County Public Schools Instructional Aids Center, 

7655 West 10th Avenue, Lakewood, Colo. 80215. 



A COMPUTER ORIENTED APPROACH TO 
LINEAR EQUATIONS IN ONE UNKNOWN 
by Robert A. Norton and Roe E. Willis 

Today’s junior high school student is destined to have some con- 
tact with computers and computer programming before high school graduation. 
This is the impetus for writing this text. 

The basic tool in this text is the flow chart which is used as a 
means of teaching solutions of single variable linear equations. The 
text is written for the student to read^ reason through ^ and work out 
the exercises. Eighth or ninth grade students are best suited for the 
content, although an aggresive seventh grader could be expected to succeed. 

The flow chart is the basic step for computer programming. The 
self-teaching properties of the flow chart coupled with its motivational 
value will give the student a good basic foundation. 

This is not intended as a package text. It is a preliminary step 
in using a computer oreinted concept to teach a section of elementary 
algebra. Other topics readily lend themselves to this technique. In- 
equalities, systems of linear equations, and quadratic equations are 
such examples. 

The Flow Chart 

The equation x + 3 = 13 is asking you to find a number- to replace 
X which will make the sentence (equation) true. The obvious answer is 
10. The aim here is learning to solve much more difficult problems. 

This requires an organized approach to solving a simple problem and will 
lead to more complicated problems using this basic approach. 

Figure 1 is the flow chart solution of this equation. 

There are five mathematical sentences in the rectangular boxes. 

If X is replaced by 10 in each of them, the equations are satisfied or 
true. When two or more sentences have exactly the same truth set, the 
sentences are equivalent . We have a way of making certain every step we 
make gives an equivalent equation. The properties which are given in 
Appendix 1 are the necessary tools for this. A step is not allowed 
unless there is a property to support it. When we have to justify or 
give a reason for each step this is mathematics.- 

From here on you, are required to justify each step you make by 
stating the property used. Remember equivalent equations mean they 
have the same truth set. 

Figure 2 shows the first solution and gives the properties for 
each step. 
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A lot of work for a simple problem? Yes, it is; but now a 
pattern is set. Every problem like this is solved by the same steps. 
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The foregoing problem involved only properties of addition. 
Here is a similar problem involving only multiplication. Solve 
-3x =24, and give the properties for each step. Refer to Figure 3. 




C 



END 



3 



Figure 3 
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Exercises ^ 

Flow chart the solutions, and give the properties for each step 
in the solutions of the following problems, 
a) 3 = X - 18 b) 4x = -22 c) c - 4x = 9 

If your flow chart solution does not agree with the answer given, 

carefully check and compare your steps. You are developing two tech- 
niques: 1) steps in solving equations and 2) an understanding of why 

we can make each step. 

We slipped in a property in Problem la which had not been mentioned. 
The properties only provide for addition and multiplication. Subtraction 
is defined with addition. For example, 3 - 8 is the same as 3 + (-8). 

Division is defined in a similar way — 9 v (-3) is the same as 9C~^). 

This provides us with a well defined set of properties. These are suf- 
ficient for solving any linear equation using rational numbers. 

The next equation to be considered is one like Ic involving both 
addition and multiplication properties. The solution of 13 = 4 - 5x 
is shown in Figure 4. 

In checking the answer, the order of operations is: 1) Do all 

operations within the parentheses first, 2) all multiplications, and 
3) all additions. 

Exercises II 

Flow chart the solutions, and give the ' properties for each step, 
a) 3x + 2 = 7 b) 9=5-^ c) 2^x - 17 = 3 

The next problem to be solved is 2x + 3x = 20. The left side of 
the equation 2x + 3x can be rewritten as (2 + 3)x. The distributive 
property justifies this step. It simplifies further to 5x using the 
closure property of addition. This problem could occur in the following 
forms and have the same basic solution: 1) 2x + 3x = 20, 2) 20 - 3x = 2x, 
and 3) 3x = 20 - 2x. The flow chart solution for 2) is shown in Figure 5. 

Now that the distributive property is explained, START means to 
look for all terms containing the variable, x has been used as our vari- 
able throughout, but y, z, etc., could have been used as well. In plan- 
ning the solution ahead, all terms containing the variable must end up 
on the same side. Use of the distributive property then enables these 
terms to be collected into a single term. Figure 6 shows the solution of 
3x - 13 = 7x + 4x. 





Figure 4 
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An ecjuation liks 3(x + 2) = 21 involves the distributive property 
in reverse order.. Refer to Figure 7. 
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It should be mentioned at this point, for 3 + 7x + 8 = 21, 

3 and 8. should be associated and added into a single term. The 
result would be 7x + 11 - 21 or 11 + 7x = 21. 3, 8, 21 are numeri- 

cal terms and 7x a variable term. Always plan your solution 
ahead with two objectives in mind: 1) Collect all numerical terms 

into a single term, and collect all variable terms into a single 
term; and 2) pick the shortest route with the least number of 
steps. 

EXERCISES III 

Flow chart the solutions, and give the properties for each 
step in the following exercises: 

a) 3(x - 4) = 21 b) 4x - 8x = -72 c) 4 - 3x = 2x + 9 
d) ^ - 2 = 4(2 - e) 2(x - 3) = 7 - 2x 

3 D 

These few pages have exposed you to a computer oriented 
concept — the flow chart. The flow chart is a way or means of 
organizing a solution. More important is the actual mathematics 
involved. The solution appears to be our main objective, but 
like any problem if the steps have no support or reason the solution 
means nothing. Problem solution is formally named applied mathe- 
matics. One branch can be emphasized more than the other, but 
it is impossible to study one without in some way considering the 
other. 



Merging of applied and pure mathematics has given us today’s 
high speed computers. The applied mathematician provides the 
steps it must go through, while the pure mathematician provides 
the logic or reasoning required to achieve these steps. 

Once you have successfully understood the properties and 
their use, it is to your benefit to drop this formality and develop 
efficient equation solution methods for yourself. 



Consider Thirds: 



~ - 2 = 4(2 - -) 
3 6 
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This is a quick and efficient method. Why all the careful 
step-by-step reasoning? If every problem you encounter is no 
more difficult or different than these, there would be no purpose. 
However, a new and completely different pixDblem can be successfully 
solved if a step-by-step solution is organized with a property 
for each step. 

The computer has a language it understands. If the engineer 
can break his problem down to the ccmiputer’s language, the solution 
is performed in seconds compared to days. The flow chart is the 
first step in organizing a problem before putting it into computer 
language . 

Already computers have reached the speed of light, which 
isn^t fast enough for complex weather problem solutions to be 
ready for the evening news. Finding different methods for solving 
the same problem and continued redesigning of the computer network 
will keep the two branches of mathematicians busy for some time. 

Don’t forget — they had to begin where you are now. 
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Properties and Definitions Used to Establish 
Simpler Equivalent Equations in the Rational Number Field 

Properties in the Rational Number Field 

1) Closure property of addition: If a and b are any rational 

numbers, then a + b is exactly one rational number. 

2 2 

Examples: 2 + ~ = 2— ; 4 + 10 = 14 

2) Closure property of multiplication: If a and b are any rational 

numbers, then (a)(b) is exactly one rational number. 

Examples: (2)(|-) = |- ; ^ I (3)(2) = 6 

3) Commutative property of addition: If a and b are any rational 

numbers, then a + b = b + a. 

2 112 

Examples : ’ 3 ‘‘'' 5 =' 5 +J J 5 + 7 = 7+ 5 

4) Commutative property of multiplication: If a and b are any 

rational numbers, then (a)(b) = (b)(a). 

Examples: (2)(3) = (3)(2) ; (f)(4) = (4)(|) 

5) Associative property of addition: If a, b, and c are any 

rational numbers, then (a + b) + c = a + (b + c). 

Example: (2+5) +4=2+ (5+4) 

Verify: 7 + 4 = 2 + 9 

11 = 11 True 

6) Associative property of multiplication: If a, b, and c are 

any rational numbers, then ((a)(b))c = a((b)(c)). 

Example: ((2)(3))(4) = (2)((3)(4)) 

Verify: (6)(4) = (2)(12) 

24 = 24 True 

7) Identity property of addition: If a is any rational number, 

then a+0=0+a=a. 

o 9 

Examples: 5+0=5 ; 0+~=~ 

3 3 

i.e., zero is the identity element for addition. 

8) Identity property of multiplication: If a is any rational 

number, then (a)(1) = (l)(a) = a. 

22 22 

Examples: (•—)(!) = ; (1)(2) = 2 

i.e., one is the identity element for multiplication. 
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9) Inverse property of addition: If a is any rational number 

and b is its opposite (additive inverse), then a + b = 0. 

Examples: 2 + (-2) = 0 ; = o 

o o 

10) Inverse property of multiplication: If a is any rational 

number except zero and b is its multiplicative inverse, 
then (a)(b) = 1. 

Examples: (5)(~) = 1 ; (2~)(|-) = 1 

11) Distributive property of multiplication over addition: If 

a, b, and c are any rational numbers, then a(b + c) = (a)(b) + (a)(c). 
Example: 3(4 + 5) = (3)(4) + (3)(5) 

Verify: 3(9) =12+15 

27 = 27 True 

12) Multiplication property of zero: If a is any rational niomber, 

then (a)(0) =.0. 

Examples: (0)(2) = 0 ; = 0 

13) Addition property of equality: If a^ b, and c are any rational 

numbers and a = b, then a + c = b + c and is equivalent to a = b. 

4 4 

Examples : ^ and 2 + 5 = — + 5 are equivalent . 

X = 5 and x + 3 = 5 + 3 are equivalent. 

14) Multiplication property of equality: If a, b, and c are any 

rational . numbers , c 0, and a = b, then (a)(c) = (b)(c) and 

is equivalent to a = b. 

15 ' 15 

Examples: 5 = ~ and (5) (7) = (-~)(7) are equivalent. 

I 

X = 5 and- 3x = 15 are equivalent. 

Definitions in the Ratio nal Number Field 

t 

1) Solution set: replacement values for the variable which make 

the statement true. 

2) Division by zero: undefined in the field of rational numbers. 

3) Equivalent statements or equations: statements which have 

the same solution set. 

4) Subtraction defined in terms of addition: a - b = a + (-b) 

Examples: 2-4=2+ (-4) ; 1 - | = i + (_ i) 

o o o 3 

5) Division defined in terms of multiplication: ^ = (^)(a) 

Examples: | = (|)(2) ; | = (i)(x) 
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PRIME NUMBERS AND "BASIC" COMPUTER LANGUAGE 
by Ray Platt and Ken Marine 

*1" ■»-»i -4- “ ■— I — ^ 

JLLl UXOil 

This is a project for one or two class periods in a Junior or Senior 
High matheiDatics course. The purpose is to show that a computer can help in 
determining prime numbers. The individual students should knew the terms 
square root, natural numbers, factor, divisor, multiple, even and odd. 

The computer approach will be incorporated with the prime numbers topic 
used in most mathematics text books. Two flew charts and two programs are 
given for each section of this topic. The students might generate others 
that fit their individual program or procedures. 

This project consists of three parts. Part I is a brute force attack 
on a number. Part II is a more sophisticated method cf determining if a 
number is prime. Part III is a section to list prime numbers with the use 
of the Sieve of Eratosthenes. 

Definition 

A natural number with only two different factors is called a prime 
number . 

PART I 
Brute Force 

First ask the students bow they can tell if a number is prime. Give 
them some numbers such as 1, 2, 4, 7, 9, and 15. Then give, them a large 
number such as 529 which seems to be prime but is dsvisible by 23. Then 
ask the students how they worked on the numbers given. This should lead 
to a flow chart such as Flow Chart 1 or- 2. In this attack, a number is 
divided by every number starting with 2 up to the given number. Program 
1 prints out which numbers are prime and which ones are not and lists two 
possible factors. Program. 2 prints out factors of the given number other 
than one and the number itself if such factors exist. If there are no such 
factors, it prints out one and the number itself. Inspection of the results 
reveals which numbers are prime. Students of your’ class might be asked to 
compare the run time of this type of program with a more refined program of 
Part II. Note: Try a large prime such as 2351. 

PART II 

Division by Numbers Less than or Equal to xhe 
Square Roox of a Given Number 

After the students start dividing by all the numbers less than the 
number in question, you can ask them to find a shorter method. Take a 
number such as 53. Play computer using one of the programs given in Part I. 



o 
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This might lead to the idea that a limit of the square root of the given 
number can be set on the divisor. (See Flow Charts 3 and 4). The programs 
given here show two ways in which the square root step is used. One 
program starts with the square root and works down to one. The other works 
from 2 as the first divisor. (See Programs 3 and 4). Note: The run time 
is the same for Part I and Part II programs. Try other numbers to make 
a good comparison. 

PART III 

The Sieve of Eratosthenes 

The Sieve of Eratosthenes is commonly introduced by having the students 
make a chart (Figure A) which lists the numbers from 1 to 100. First the 
number one is crossed out, because it does not fit the definition. Then 2 
is circled as a prime. Then each multiple of 2 is crossed out. Then 3 is 
circled, and all of the multiples of 3 are crossed out. This process is used 
until the numbers that are circled are considered to be prime. This idea 
leads directly to the flow charts for this section (5 and 6). In each case 
the only divisors used are prime numbers. In Program 5 the square root limit 
is used. Note also the different methods used for storing the primes and 
printing the primes. (See Programs 5 and 6). It must be mentioned that the 
dimension of the total of primes (see step one of the programs) is a point 
for the students to determine. You also may wish to print out the total at 
the end of your list of primes. 

PROGRAM 1 



10 -READ M 
20 LET N=2 
30 LET K=INT(M/N) 

40 LET L=M-K*N 
50 IF L<>0 THEN 80 

60 PRINT M,"IS N0T PRIME;FACT0RS ARE",K,N 

70 G0 T0 10 

80 LET N=N+1 

90 IF N<M THEN 30 

100 PRINT M,"IS PRIME" 

110 G0 T0 10 

120 DATA 31,121,123,113 

130 END 

RUN 

WAIT. 



Ml 15:01^ PX THU 08/17/67 



31 


IS PRIME 






121 


IS N0T PRIME ;FACT0RS ARE 


11 


11 


123 


IS N0T PRIME ;FACT0RS ARE 


41 


3 


113 


IS PRIME 
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0UT 0F DATA IN 10 
TIME: 0 SECS. 



PROGRAM 2 




10 READ N 
20 LET D=2 
30 LET Q=N/D 

40 IF INT(Q)*D=N THEN 70 
50 LET D=D+1 
60 G0 T0 30 

70 PRINT N;"=";Q;"*";D 
80 G0 T0 10 

90 DATA 31,121,123,113 

100 END 

KEY 

READY. 



RUN 

WAIT. 



PI 14:49 PX THU 08/17/67 




31 


= 1 




31 


121 


= 11 




11 


123 


= 41 


rv 


3 


113 


= 1 




113 


0UT 0F DATA 


IN 


10 



TIME: 0 SECS. 



PROGRAM 3 



10 READ M 
20 IF M<2 THEN 130 
30 LET N=2 
40 LET K=INT(M/N) 

50 LET L=M-K*N 
60 IF LoO THEN 90 

70 PRINT M, "ENTER ONLY INTEGERS >=2" 
80 G0 T0 10 
90 LET N=N+1 

100 IF N<=SQR(M) THEN 40 
110 PRINT M,"IS PRIME" 

120 G0 T0 10 

130 PRINT M,"IS WR0NG" 

140 G0 T0 10 

150 DATA 31,121,123,113 

160 END 

RUN 

WAIT. 

L 15:03 PX THU 08/17/67 






31 


IS PRIME 


121 


IS N0T PRIME 


123 


IS N0T PRIME 


113 


IS PRIME 


0UT 0F 


DATA IN 10 


TIME: 


0 SECS. 



PROGRAM 4 



10 READ N 

20 LET D=INT(SQR(N)) 

30 IF INT(N/D)*D=N THEN 90 

40 LET D=D-1 

50 IF D=1 THEN 70 

60 G0 T0 30 

70 PRINT N;"IS PRIME" 

80 G0 T0 10 

90 PRINT N;"IS N0T PRIME" 

100 PRINT D;"*";INT(N/D);"=";N 

110 G0 T0-1O 

120 DATA 31,121,123,113 

130 END 

RUN 

P3 14:59 PX THU 08/17/67 



31 


IS 


PRIME 


121 


IS 


N0T PRIME 


11 


* 11 = 121 


123 


IS 


N0T PRIME 


3 


?’« 1 


n = 123 


113 


IS 


PRIME 



0UT 0F DATA IN 10 
TIME: 0 SECS. 



PROGRAM 5 



10 DIM A(50) 

20 LET A(l)=l 
25LET A(2)=2 
30 F0R M=3 T0 99 
40 F0R N=2 T0 SQR(M) 

50 IF A(N)=0 THEN 100 
60 LET K=INT(M/A(N)) 

70 LET L=M-K*A(N) 

80 IF L-0 THEN 150 

90 NEXT N 

100 F0R 1=1 T0 50 

no IF A(I)oO THEN 140 

120 LET A(I)=M 

130 G0 T0 150 

140 NEXT I 

150 NEXT M 
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160 F0R 1=2 T0 50 
170 IF A(I)=0 THEN 200 
180 PRINT A(I); 

190 NEXT I 
200 END 



RUN 



M 




15:10 


PX THU 


08/17/67 












2 


3 


5 


7 


11 


13 


17 


19 


23 


29 


31 


37 


41 


43 


47 


53 


59 


61 


67 


71 


73 


79 


83 


89 


97 


















TIME: 


2 


SECS. 



















PROGRAM 6 



05 DIM X(50) 
lOLET Z=1 
20 READ N 
30 F0R A=3 T0 N 
40 F0R M=1 T0 Z 
50 LET X(l)=2 
60 LET Q=INT(A/X(M)) 

70 IF A-(X(M)*Q)=0 THEN 110 

80 NEXT M 

90 LET Z=Z+1 

100 LET X(Z)=A 

110 NEXT A 

120 F0R M=1 T0 Z 

130 PRINT X(M); 

140 NEXT M 
150 DATA 100 
160 END 



RUN 

P2 14:56 PX THU 08/17/67 

2 3 5 7 11 13 17 19 23 29 31 

37 41 43 47 53 59 61 67 71 73 79 

83 89 97 

TIME: 2 SECS. 



Summary 

Those programs and charts are guides in the construction of material 
for determining prime numbers. In each case refinements can be made on 
these programs. During the writing of this project, tests were made on 
larger numbers. It was found that the run time can be cut by several 
seconds if only odd numbers are considered for primes. This step can 
easily be inserted with a Step Statement of 2. 
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Each section of this project can be used separately if so desired. 

All three sections will work in connection with most mathematics texts 
which cover prime numJbers. 




FIGURE A 
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FLOW CHART 1 



FLOW CHART 2 
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FLOW CHART 3 
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AN APPROACH TO THE INTRODUCTION OF 
COMPUTER PROGRAMMING IN 
JUNIOR HIGH SCHOOL MATHEMATICS 
by Helen E. Schweizer 

The purpose of this project is to develop materials for a junior 
high school student which would reinforce and extend his concepts of num- 
ber and number operations ari help him to develop a familiarity with com- 
puter programming. 

It is assumed that the student has identified and operated with the 
sets of natural, whole, and non-negative rational numbers (both common and 
decimal fractions) and has an awareness of negative numbers; that he has 
had experience v^orking with several number bases — place value, s 3 rmbols, 
conversion from one base to another; and that he has made limited use of 
an exponent. Primarily, the number operations v^hich he would have used 
are addition, subtraction, multiplication, and division. 

The study guides which follow begin with a mathematical definition 
appropriate for this level, or with a suitable concept, and attempt to 
relate this to computer techniques. They are intended as curriculum guide- 
lines and, in general, might be used as a basis for class discussion or 
for individual study by a student. In some instances, an attempt has been 
made to provide student material and problem exercises, but the knowledge- 
able teacher will realize that these need to be augmented with appropriate 

class work. 

Study Guide I 
Number 

1. What is the first counting (natural) number? Write the second 

counting number. How could it be obtained from the first number? 

Write the third counting number. How could it be obtained from 

the second? 

How can each counting number be obtained from the previous counting 
number? 

What is the largest counting number you know? 

(Teacher: Do not accept "infinity." By showing that infinity does 

not obey our rules for number, demonstrate to the student that infinity 

is not a number.) 

2. The set of whole numbers contains all the counting numbers and one 
other number. What is this number? 



193 



Computer 

With certain ccxnmands the computer can produce those elements of 
the set of counting numbers or the set of whole numbers which are speci- 
fied. This is the program for the first 100 v/hole numbers. 

1 LET X = 0 

2 PRINT X 

3 LET X = X + 1 

4 IF X <= 99 THEN 2 

5 END 

The following comments explain the above program in detail. 

1. Each line begins with a number. 

2. The sequence of the line numbers indicates to the computer the 
order of steps to follow. 

3. X is a variable . A variable in BASIC may be indicated by any 
letter. A variable stands for a number. 

4. "LET" is one of the statements used to assign a value to a variable; 
it is noL used to state an equality. The value of the variable does 
not change until the next LET statement is encountered with a value 
for that variable. 

PRINi X tells the computer to print each value for X. 

6. The "IF" "THEN" statement places a limit on the size of the set. 

7. "END" tells the computer to stop reading the program. Without this 
command the computer looks for more information. 

8. The computer would perform with the above program as follows: 

1) Let X = 0 because of line 1. 

2) Print 0 because of line 2. 

3) Change the value of X to 0 + 1 (or 1) because of line 3. 

4) Determine whether the latest value for X would be less than or 
equal to 99 because of line 4. 

5) Because the value of X would be less than 99, return to line 2 
and print 1. 

6) The present value of X would be increased by 1 because of line 3. 
The new value of X would be 2. 

7) Continue to run through the loop of lines 2, 3, and 4 until X 
would not be less than or equal to 99. 

8) At this time continue to the next line, line 5, and stop. 
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Computer Activity I 

A teacher may be able to create interest and promote student involve- 
ment by using students to simulate a computer program. Some students, 
particularly those who clearly understand the computer program for whole 
numbers, could function as programmers. The other students in the class 
could respond to the commands of the program. 

A group of six students with the assignment of a computer program 
for the first 15 counting numbers might perform as follov;s : 

1. One student would write a complete program similar to that 
for the set of whole numbers . 

2. Each of the other 5 students would be assigned a different line 
of the program; his identification would be the number of his 
line. He vfould determine his function in the assigned program. 

3. The programmer would check individually with each of the 5 
students to be certain each had the correct command and under- 
stood his responsibility with it. 

4. Lines #1 and #3 would be given blank cards. Each card would 
be used to complete an operation v^ithin the computer. 

5. Line ^2 would be at the blackboard to print the results requested 
by the program. 

This is the program: 

1 LET X = 1 

2 PRINT X 

3 LET X = X + 1 

4 IF X < = 15 THEN 2 

5 END 

The program might run as follows: 

1. The programmer would instruct the program to "start." 

2. ifl would write "1" on a card and would give the card to §2, 

3. #2 would print "1" on the blackboard and would give the card to #3. 

4. #3 would keep that card. On a new card he would write "2", the 

result of adding 1 to the previous value of X. He would give the card 
to #4. 

5. #4 would determine that 2 is less than or equal to 15. He would give 
the card to #2. 

6. The loop would continue to repeat until X > 15. Then #4 would give 
the card to #5. 

7. The computer wouj.d stop. 
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Upon completion of the program, each student would make a diagram 
(flow chart) to show the purpose of each command and its location in the 
total program. A discussion of these flow charts should lead to the 
recognition of the need for uniformity of symbolism and to an agreement on 
the symbols to be used in the flow charts for start , input , instructions , 
decisions, output, and stop. 

Many variations of this activity are possible. The responsibilities 
within a program may be rotated. A related mathematical concept or a 
completely different concept may serve as the basis for developing a new 
program. As new computer commands are learned, attention may be given to 
the various program possibilities which are available for a given program. 

In any case, the complexity of the activity should be dictated by the ability 
of the group involved. 

The above activity might also be utilized to develop a greater 
student appreciation of the time-saving aspect of the computer. A 
comparison can be made between the running times of the class and the computer 
for a particular program. Actual computing time can also be compared. 

Computer Activity II 

1. Write the first 15 even counting numbers. 

2. How is the second even counting number obtained from the first? 

3. How is each successive even counting number obtained from the previous 
even counting number? 

4. How could the program for the whole numbers be changed so that even 
counting numbers could be obtained? 

5. Write a program for the even counting numbers between 117 and 202. 

6. Should any changes be made in the program in to obtain the odd 
counting numbers between 117 and 202? Explain your answer. 

Computer Activity III 

1. Write five different multiples of three. What characteristic is common 
to these five numbers? 

2. Write the first five multiples of three which are larger than 39. 

3. Begin with 40, and count by 3's to 60. Are these numbers multiples 
of 3? Explain your answer. 

4. Write a program for consecutive multiples of 3. 

5. Write a program which will obtain the numbers for §3. 
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Study Guide II 
Number 

In base 10 the counting and whole numbers may be written by using 
10 different S 3 nnbols which are called digits. These digits are 0, 1, 2, 

3, 4, 5, 6, 7, 8, 9. The position of these digits in the numeral deter- 
mines their value. The place values of the first six places to the left 
of the decimal point in base 10 may be expressed as follows: 

100000 10000 1000 100 10 1 

10 ^ 10 ^ 10 ^ 10 ^ 10 ^ 10 ° 

The value of 29056 could be expressed as: 

2(10000) + 9(1000) + 0(100) + 5(10) + 6(1) 
or as: 

2(10*^) + 9(10°) + O(IO^) + 5(10^) + 6(10°) 

In base 8 the counting and whole numbers may be written by using 
eight different S 3 mibols (0,^ 1, 2, 3, 4, 5, 6, 7) and the place values 
may be expressed as powers of eight. The place values of the first six 
places in base eight are as follows: 



base 


8: 


100000 


10000 


1000 


100 


10 


1 


base 


8: 


10° 


lo'^ 


10° 


iq2 


10^ 


10° 


base 


10: 


32768 


4096 


512 


64 


8 


1 


base 


10: 


8^ 


O 


8° 


8^ 


8^ 


8° 




The 


value of 2314g 


may be 


expr-essed 


in base 


8 as 


follows 



2(10°) + 3(10^) + I(IO^) + 4(10°) 



or in base 10 as : 

2(8°) + 3(8^) + 1(8^) + 4(8°) 

2(512) + 3(64) + 1(8) + 4(1) 

1024 +192+8+4 
1228 

In base 2 the counting and whole numbers may be written by using 2 
different S 3 mibols (0, 1) and the place values may be expressed as powers 
of two. The place values of the first six places in base two may be ex- 



pressed as 
base 2: 


follows : 
100000 


10000 


1000 


100 


10 


1 


base 2: 


,^101 


n^lOO 


n^ll 




, 1 


0 


10 


10 


10 


10 


10 


10 


base 10: 


32 


16 


8 


4 


2 


1 


base 10: 


2^ 


2 “ 


CO 

CN 


CN 

CN 


2 ^ 


2° 
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The value of lllOOl^ may be expressed in base 2 as follows : 

1(10^°^) + 1(10^°°) + 1(10^^) + 0(10^°) + (10^) + 1(10°) 
or in base 10 as: 

1 ( 2 ^) + 1 ( 2 ^) + 1 ( 2 °) + 0 ( 2 ^) + 0 ( 2 ^) + 1 ( 2 °) 

1(32) + 1(16) + 1(8) + 0(4) + 0(2) + 1(1) 

32+16+8+0+0+1 

57 



Computer 

The computer places the following restrictions on the numerals 
which it will accept; 

1. All numerals should be base 10 numerals. If they are not, they are 
read as if they were in base 10. Knowledge of base 2 and base 8 

is necessary for understanding the internal operations of the com- 
puter. 

2. Numerals which contain commas are not accepted. 

3. The maximum number of digits which will be accepted varies with the 
computer. Assume that the maximum number of digits accepted by this 
computer is six. If a numeral contains more than 6 digits, it may 
be rewritten with 6 or less digits and a power of 10. 



4650000 



465000(10^) 

46500(10^) 

4650(10°) 

465(10^) 

46.5(10°) 

4.65(10°) 

.465(10^) 

.0465(10°) 



The computer will accept the powers of 10 only if they are indicated 
with the letter E followed by the exponent. The above would be rewritten 
as follows for computer use; 



465000(10^) 


465000E1 


46500(10^) 


46500E2 


4650(10°) 


4650E3 


465(10^) 


465E4 


46.5(10°) 


46.5E5 


4.65(10°) 


4.65E6 


.465(10^) 


.465E7 


.0465(10°) 


.0465E8 



The letter E must alwa 3 ?’s be preceded by at least one digit. 



o 
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Computer Activity I 
Play computer! 

V/hich of the following will you accept? 



1. 


9405 


2. 


32410- 

o 


3. 


208,507 


4. 


1059000 


5. 


lOllOO^ 


6. 


7690000000 


7. 


23173 


8. 


60781 


9. 


lOOllOl^ 


10. 


32000000 


• 

« — 1 
t — 1 


5623 


12. 


10000000000 



Rewrite the numerals which were rejected in a form which will he accepted. 

Study Guide III 
Number Operations 

1. List the operations which may be performed with counting numbers. 

2. Given the counting number A and the counting number B, v/ill the 

answers to each of the following problems always be a counting number? 
A + B 
A - B 
A X B 
A V B 
A® 

3. If the answer to any of the above is no, write the conditions which 
would be necessary to obtain an answer which is a counting number. 

Computer Activity ^ 

Make 5 copies of the following program : 

1 READ A, B 

2 PRINT A, B, A B 

3 DATA 12, 2, 15, 3, 10, 5 

4 DATA 8, 4, 2, 5, 9, 2 

5 GO TO 1 

6 END 

1. Replace the blank separating A and B in line 2 in each of the programs 
with a different one of the following number operation symbols : 



ERIC 




2 . 



Examine the results which are printed by the computer to determine 
the operation which is indicated by each of the symbols. Write 
each symbol with the name of the operation that it indicates . 

3. Are the results for each operation counting numbers? Do these 
results verify your earlier conclusions concerning this? 

Interpretation of the above program: 

1. "READ A, B" identifies the variables with which the computer will 
operate. The variables are separated with commas. 

2. The READ statement must always be followed later in the program 
vrith a DATA statement . The DATA statement gives the values of the 
variables in the order in which they are listed in the READ state- 
ment and separates them with commas. 

3. In this program the operation to be performed on the variables A 

and B is indicated in the PRINT statement. The computer will print 
the values of A, B, and A B. 

4. After the computer has used the 12, 2 values for A, B, it will go 

to line 5 which tells it to GO TO 1. The READ statement tells it 

to look for new data so 15, 3 are used as the new values for A, B. 

A, B, A B are printed. The computer will continue to work through 

the loop of lines 1, 2, 3, 4, 5, When there is no more DATA in 
line 3, it will read the DATA in line 4. When there is no more 

DATA in line 4, it will print "OUT OF DATA IN 4" and stop. 

Computer Activity II 

1. Solve the following problems by computer. 



1) 


5*2+4 


6) 


40/10*5 


11) 


23-7*3 


16) 


16/8+2 


2) 


5+2*4 


7) 


11-2+3 


12) 


23*7-3 


17) 


3+2*5 


3) 


8/2+2 


8) 


11+2-3 


13) 


4+2+3 


18) 


3*2+5 


4) 


8+2/2 


9) 


90-15/5 


14) 


4+2+3 


19) 


6+4-3 


5) 


40*10/5 


10) 


90/15-5 


15) 


16+8/2 


20) 


6-4+3 


By 


examination 


of the 


answer of each 


of the 


: above problems , ( 


determine 



which operation was performed first in each problem. Circle the 
symbol for that operation. 

3. What order do you think the computer would follow in solving a problem 
in which all five operations were involved? 

4. Write a problem which involves all five operations. 

Solve this problem by computer to test your answer to #3. 
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COMPUTER-ORIENTED SOLUTION OF A BALLISTIC PROBLEM 
by Bill Tipton and Marvin Miller- 

The solution of problems involving two dimensional motion (non- 
straight line motion), v^herein numerous computations must be made under a 
variety of given conditions, is readiuy adaprable xo a computer program. 
Application of the computer gives a flexibility for dealing with a great 
number of changing situations in a given type cf problem. The following 
example illustrates a basic situaxicn v/herein a bcdy is prcjected with a 
given initial velocity in some direction as determined by the angle 
measured from the horizontal. As in most problems in the elementary 
stage, only those factors which are readily determined by secondary school 
laboratory techniques are considered. Thus air-friction, coriolis 
effect, etc, are not considered in this example. If these are to be 
introduced at a later stage of investigation, the basic program would 
remain the same, making rt necessary only to add the instructions in the 
program to account for any changes cf motion due to these factors. 

This example will require the student xc develop an expression such 
that the maximum range is a function of the initial velocity (V) and the 
angle of projection (A) if he wishes to work with only one function. In 
this example we wish to examine other variables, i.e., time of flight (T) 
and maximum altitude (Y); therefore, other functions are necessary. The 
conversion of radians to degrees (and having this result in tabular form) 
is also an extraneous portion of the pregram so far as determining the 
range; however, it may be of interest to illustrate the value of the 
angle in each of the two units.. 

As may be noted, the use of constants having seven place accuracy 
results in an error of four parts in ten thousand after ninety solution 
cycles. However, this result may be beneficial; as it illustrates the 
error introduced into a computation by use of rounded numbers. 




involving other types of motion at the earth’s surface. Also, the reader 
can either extend or limit the scope of the program as given here by 
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making desired changes in the initial velocity or by reading in only 
those angle values for which he desires a solution. For example, if the 
initial velocity is 3000 cm/sec at an angle of 30 degrees, this would 
necessitate a change in the program as follows: Assign the new value 

(3000) in line number 10. Also assign the new value (30) in line number 
50 and remove lines numbered 130 and 140. The resulting program would 
then determine the range, altitude, and the time of flight of the 
projectile for that particular angle of elevation and initial velocity. 

The student should familiarize himself with the meaning of these 
formulas: A = D(. 01745329), where A is the measure of an angle in radians, 

and D is the measure of that angle in degrees. (The con- 
version factor used here is only approximate , not exact . ) 

2 

R - 2V cos A sin A/G, where V is the initial velocity 

(expressed in cm/sec), A is the measure of the launch angle 

(in radians), G is the gravitational force toward the earth 
2 

(980 cm/sec ), and R is the horizontal range of the projectle 
(in cm); 

T = R/V cos A, where T is the total time in flight (in seconds); 

1 2 

Y = V sin A (t) - j Gt , where t is one-half of T (this would be 
the time elapsed from launch until the projectile reaches its 
apogee); and Y is the height of the apogee (the maximum 
height attained, in cm). 

In the following discussion of the program, it has been assumed that 
the student is familiar with the BASIC language which is a very widely 

used medium for teletype communication with a time-sharing computer such 
as the G-E 235. 



Having understood these formulas (and perhaps performed a few trial 
computations by conventional methods), the student must concern himself not 
only with the problem of instructing the machine to substitute his data 
(known values for variables V, D, and G) into the above formulas, and sub- 
sequently calculate corresponding values for A, R, T, and Y, but also he 
( 2 ) 

will want to program the necessary steps which will cause the machine 



to arrange the output information according to an orderly, easy-to-read 

format. And, as is mentioned elsewhere, if the student desires to test 

varying conditions in the launch angle and initial velocity, he should 
(3) 

create a loop in the program which will automatically cause the machine 
to repeat its cycle of computations until all the given data is exhausted. 

The part of the problem identified as (1) above (i.e., substituting 
known values for some variables and computing corresponding values for the 
other variables) is achieved by lines 10 and 50-110 in the program. 







rgppg 






202 

The instructions which cause the output to be printed in tabular 
form (objective 2 above) are lines 20-45 and 120. (Line 45 simply 
causes the teletype to skip one vertical line space). 

Note that the commas in line 120 instruct the machine to space the 
columns of numerical output data in 15-space horizontal intervals. That 
is, the left sides of the 5 columns will begin at the first, six tenth, 
thirty-first, forty-sixth, and sixty-first spaces on the 75 character 
width of the printing space. Therefore, care must be taken when typing 
lines 30 and 40 to cause the column headings to be properly placed. 

Objective 3 is fulfilled by lines 130 and 140. This loop is 
conditional - the computer continues to increase the measure of the 
launch angle by 1 degree until the angle becomes 91 degrees. The data 
resulting from such a case would duplicate that of an 89 degree angle. 
Therefore, we have chosen the condition that the angle exceed 90 degrees 
as that which directs the computer to "jump" out of the cycle of: 

(1) substituting 6 evaluating, (2) printing information, and (3) repeating 
with a larger angle measure, by line 150, the end statement. 



10 LET V= 1000 

20 PRINT"IF THE INITIAL VEL0CITY IS A C0NSTANT" V"CENTIMETERS/SEC" 

30 PRINT" AND IF THE ANGLE IS THEN TIME IN FLIGHT RANGE ALT" 

40 PRING"DEGREES RADIANS SEC0NDS CM CM" 

45 PRINT 

50 LET D=0 

60 LET G=980 

70 LET A = D*. 01745329 

80 LET R= 2*Vi2*C0S(A)*SIN(A)/G 

90 LET T = R/(V*C0S(A)) 

100 LET T=T/2 

110 LET Y=V*SIN(A)*T-1/2*G*T*T 

120 PRINT D,A,T,R,Y 

130 LET D=D+1 

140 IF D<=90 THEN 70 

150 END 

RUN 

BALIST 10:25 PX M0N 08/21/67 



IF THE INITIAL VEL0CITY 
AND IF THE ANGLE IS 
DEGREES RADIANS 

0 
1 
2 

3 

4 

5 

6 



IS A C0NSTANT 1000 
THEN TIME IN FLIGHT 
SEC0NDS 



CENTIMETERS/SEC 



RANGE 




CM 


CM 


0 


0 


35.6117 


.155401 


71.1801 


.621416 


106.662 


1.39748 


142.013 


2.48263 


177.192 


3.87557 



0 

1.74533 E-2 
3.49066. E-2 
5.23599 E-2 
6.98132 E-2 
8.72665 E-2 
(etc. ) 



0 

1.78086 E-2 
3.56117 E-2 
.053404 
7.11801 E-2 
8.89344 E-2 
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"IF THE INITIAL VELOCITY IS A CONSTANT" V "CENTIMETERS/ SEC 
AND IF THE ANGLE IS THEN TIME IN FLIGHT RANGE ALT 
DEGREES RADIANS SECONDS 



CM 



CM 



D = 0 






G = 


980 



A = D(. 01745329 



R = 2V^ Cos A Sin A/G 


- 




T = R/(V Cos A) 1 


( T = 


T/2 i 






Y = V Sin A • T - 


2 

G T 






A COMPUTERIZED APPROACH TO GRAPHING PARABOLAS 
by Lavere C. Wilson 
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2 

In this unit, functions of the form [(x,y): ax + bx + c = y] 

where and c_ are real numbers and a 0 will be discussed. The 

graph of each such function is called a parabola. Although there are 
parabolas which are not graphs of functions , this discussion will involve 
only those which are. 

The parabolas under discussion here have either a maximum or a 

minimum value, a vertex point, and an axis of symmetry. The parabola 

in Figure 1 has a minimum value of N, the coordinates of the vertex point 

are (M,N), and the axis of symmetry is defined by the equation x = M, 

where M is the x coordinate of the vertex point. In Figure 2, the vertex 

point and the axis of symmetry are the same as in Figure 1. The parabola 

in Figure 2 has a maximum value represented by N. 

2 

In the equation, ax + bx + c = y, where a ^ 0, the coefficients 

a, and £ determine certain things about the parabola. If b = 0 and 

2 

c = 0, the resulting equation, ax = y, is the equation of a parabola 

which has its vertex point at the origin, (0,0). The axis of symmetry 

is the y-axis , as in Figure 3 , If b = 0 , resulting in the equation 
2 

ax + c = y, the value £ is positive, the curve is moved up £ units; 
and if £ is negative, the curve is moved down |c| units. See Figures 
4 and 5. 

2 

If b 5 *^ 0 in the equation ax + bx + c = y, then this causes a hori- 
zontal displacement of the parabola from a position symmetric with the 
y-axis . 



In summary, the following may be said about the effect of the values 
of a, and £ upon the graph of the function C(x,y): ax +bx+c=y,a=0] 

The value of £ determines whether the parabola will have a maximum or mini- 
mum value. If a < 0, the parabola will have a maximum value, and if a > 0 
the parabola will have a minimum value. The value of b determines whether, 
the axis of symmetry is to the left or right of the y-axis, and the value of 

£ affects the vertical position of the parabola, moving its graph up or down. 

2 

By completing the square with the equation ax +bx+c=y, (aT^O), 
information about the vertex point, the maximum or minimum value, and the 



axis of symmetry may be obtained. 

2 

ax + bx + c - y where a / 0 

a(x + — x) + c = y 
d. 

/2^b ,b.2. ^ /t).2 



( 1 ) 

( 2 ) 



(3) 






^ ^^rrs^-rfS^^^ ’*^. -p- 
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In equation (4), if x = -b/2a, then either a maximum or a minimum value 

will be obtained which determines the vertex point c Since the axis of 

symmetry passes through the vertex, then x = -b/2a determines the value 

of the X coordinate of the vertex point c Hence, x = -b/2a is the equation 

for the axis of S 3 nnmetry. The value of the y coordinate of the vertex 

2 2 

point is c - (b /4a ) . Whether this vertex point represents a maximum 

or a minimum point of the parabola may be determined by observing the 

2 

value of ^ in ax + bx + c = y. The parabola in question will always be 
a curve that is S3nnmetric with respect to the line called the axis of 
S3nnmetry . 

The preceding discussion is representative of The type of discussion 

found in many Algebra II texts. Students are usually asked to analyze the 
2 

equation ax + bx + c = y and drav7 graphs for several specific examples. 
Although this is a relatively simple task and certainly can be accomplished 
quite well without the use of a computer, this topic does provide an op- 
portunity to illustrate how posing problems in a general form and de- 
veloping a computer program can save considerable calculation time. To 
calculate points on the graph of y = 1.495x - 0.832x + 15.632 for values 

of X incremented by 0.25 would require more calculation than most students 
would be willing to do. However, equations which arise from "real" 
problems are likely to be even more complex than this. In addition to 
saving calculation time, developing computer programs encourages students 
to analyze carefully the mathematical relationships involved. The re- 
mainder of this unit will illustrate how flow charts and computer programs 
could be used in presenting this topic. 

Presenting problems to a computer requires defining certain symbols 
and language that the computer is able to understand or interpret. As a 
preliminary step to preparing a computer program, construction of a 
flow chart is often helpful. Some common flow chart s 3 nnbols are presented 
in Figure 6. 

A flow chart may serve as an outline for the writing of a program. 

A program is a set of directions that is used to tell a computer how to 
provide an answer to some problem. The program must be presented in a 
language that is understood by the computer, and the program must be com- 
pletely and precisely stated. Some of the programming S 3 nnbols are listed 

in Figure 7. A program designed to yield information for plotting the 

2 

graph of the equation ax = y could be written as indicated by Program I. 
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FLOW CHART SYMBOLS 




START OR STOP 



TO INTRODUCE MATERIAL TO COMPUTER 



TO GIVE COMMANDS 



TO MAKE DECISIONS 




TO OBTAIN OUTPUT RESULTS 



FLOW CHART CONNECTOR 



Figure 6 
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PROGRAMMING SYMBOLS 



= 


IS 


EQUAL TO 


< 


IS 


LESS THAN 


<= 


IS 


LESS THAN OR EQUAL TO 


> 


IS 


GREATER THAN 


>= 


IS 


GREATER THAN OR EQUAL TO 


<> 


IS 


NOT EQUAL TO 



+ ADDITION 

SUBTRACTION 
* MULTIPLICATION 

/ DIVISION 

f RAISE TO THE POWER 

( ) FORM OF INCLUSION 

Figure 7 



PROGRAM I 

10 READ A 
20 LET X = -5 
30 LET Y = A*X*X 
40 PRINT X,Y 
50 LET X = X + 1 
60 IF X <= 5 THEN 30 
70 G0 T0 10 

80 DATA l,-l,2,-2,3,-3,.5,-.5,o25,-.25 
90 END 



Compare Program I with Flow Chart I, and notice the similarities. 

The program is written in a language which resembles ordinary mathematical 
notation. This language, called BASIC, is designed for complete and 
precise specification of a program for execution by a computer. 

Observe the following about the program: 

1. The program uses only capital letters, since the teletype- 
writer has only capital letters . 

2. Each line of the program begins with a number. These numbers 
are called line numbers and serve to identijfy the lines , each 
of which is called a statement . Line numbers serve to speci:fy 
the order in which the statements are to be performed by the 
computer. 

3. Each statement starts, after its line number, with an English 
word. This word denotes the type of statement. 



209 




( START j 


\ 




/ 




INPUT 


A 



THIS STEP INTRODUCES THE VALUE OF A. 



LET X = -5 



THIS STEP ASSIGNS THE ORIGINAL VALUE OF X. 



LET Y = A*X*X 



-THIS STEP CALCULATES Y = AX' 



PRINT (X,Y) THIS STEP CALLS FOR THE PRINTING OF (X,Y). 




NO 



STOP 



) 



i 




LET X = X + 1 






J 


f 



THIS STEP INCREASES THE ORIGINAL VALUE OF X 
BY ONE UNIT. (THIS IS ^ AN EQUATION. THE* 
EQUALS SIGN IS USED HERE TO MEAN THE FOLLOWING: 
REPLACE THE PREVIOUS VALUE ASSIGNED TO THE 
VARIABLE X BY THAT VALUE PLUS ONE.) 

THIS STEP CHECKS TO SEE IF X £ 5 AFTER X HAS 
BEEN INCREASED. 

a. IF YES, IT REFERS BACK TO THE CALCULATION 
PORTION OF THE CHART TO HAVE THE NEW VALUE 
OF Y CALCULATED. THIS PART OF THE CHART 
IS CALLED A LOOP. THE FIRST TIME THROUGH 
THE LOOP THE VALUE OF X IS -5, THE SECOND 
TIME THROUGH THE VALUE OF X IS -4, THIRD 
TIME, -3, FOURTH -2, AND SO ON UNTIL X = 5. 
THE VALUE OF Y FOR X = 5 IS CALCULATED AND 
THEN THE STEP “LET X = X + 1" SETS X = 6. 
SINCE 6 IS NOT LESS THAN OR EQUAL TO 5 THIS 
ENDS THE LOOP AND LEADS TO THE STOP STATEMENT, 

b. IF NO, IT GIVES DIRECTIONS TO STOP. 



IrERiC 
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A line-by-line analysis of the program will now be presented. The 
first statement, 10, is a READ statement. It must be accompanied by one 
or more DATA statements. When the computer encounters a READ statement 
while executing a program, it will cause the variables listed after the 
READ to be given values according to the next available numbers in the 
DATA statements. In this case, A is read in statement 10; and the first 
value in the line 80 DATA statement is assigned to A. That is, A is 
assigned the value 1. 

The computer then proceeds to statement 20 and assigns the value -5 

to the variable x. The computer then proceeds to. statement 30 where it 

2 

is directed to compute a value of y according to the formula ax = y. 

(The asterisk is used to denote multiplication.) In general, a LET 
statement directs the computer to set a variable equal to the numerical 
value of the formula on the right side of the equal sign. The computer 
then proceeds to line 40 where it is directed to print x and y, which 
in this case are -5 and 25, respectively. The computer then proceeds to 
line 50 where it is directed to replace the value of x by x + 1. In this 
case, it replaces -5 by -4. Now x = -4. The computer then proceeds to 
line 60 where it encounters the statement, IF x <= 5 THEN 30. If the 
answer to the question "Is x <= 5?" is "No," the computer goes on to the 
next statement. If the answer is "Yes," as it is in this case, the 
computer goes to the line number stated after "THEN," in this case line 30. 

This begins a "loop" which will result in the computer calculating values 

2 

of ax = y and printing them out for values of x ranging from -5 to 5. 

Now, assume the computer is at line 40 and prints out 5,25. It now pro- 
ceeds to line 50 and assigns x the value 6. It then proceeds to line 60. 
Since 6 is not less than or equal to 5, it proceeds to the next line in 
the program. Line 70 directs the computer to "G0 T0 10." It returns to 
line 10, where it is directed to "READ A." It then proceeds to the first 
DATA statement and reads the next value listed in the DATA statement, 
which is in this case -1. Now the computer proceeds through the program 
again printing out values of -lx = y for x ranging from -5 to 5. It 
then returns to line 10 where it is directed to "READ A." It assigns A 
the next unused value in the DATA statement, in this case 2»and proceeds 
through the program again. After assigning all the values in the DATA 
statement to A, the computer will print "OUT OF DATA IN 10" and stop. 

Every program must terminate with an END statement. 

This program could be used by students to produce sets of ordered 
pairs for graphing. Students could draw the graphs and discuss the effect 
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of changes in A on the graph of ax = y. Time saved in computation could 
be used for discussion of flow charting, programming, and analysis of graphs. 

PROGRAM II 

10 READ A,C 

20 IF A = 0 THEN 90 

30 LET X = -5 

40 LET Y = A * Xi2 + C 

50 PRINT X,Y 

60 LET X - X + 1 

70 IF X <= 5 THEN 40 

80 G0 T0 10 

90 PRINT "NOT A PARABOLA" 

100 G0 T0 10 

no DATA 1,5, 1,-5, 2, 7, 2, -7, .5, 7. 25, -.5, -7. 25 
120 DATA -2, -10, -4. 5,20 
130 END 

Flow Chart II and Program II are for calculating points on the 
2 

graph of ax + c = y. The only real difference between Flow Chart II 
and Flow Chart I is the inclusion of a check on the value of A. This 
flow chart could be presented in blank form to students on a ditto or 
overhead projector and completed as a class project. 

In Program II, when the computer starts on line 10, "READ A,C," it 
searches for the first DATA statement in the program and assigns the first 
listed value to A and the next listed value to C. In this case, it would 
go to line 110 and pick up the value 1 for A and 5 for C. The next time 
it gets to line 10 it would assign the next two listed values to A and 
C respectively. In this case A = 1, and C = -5. This process continues 
until all DATA in line 110 is used. Then the computer continues to the 
next DATA statement and assigns values to A and C respectively from that 
statement until it runs out of data. 

Line 40 introduces the arrow "t" which, in this example, means 
raise x to the second power. Line 90 introduces the quotation marks 
following a PRINT statement. When the computer reaches this line, it 
prints out exactly what is enclosed between the quotation marks and 
proceeds to the next line in the program. DATA statements may be placed 
anywhere in the program following the READ statement, but are usually 
inserted just prior to the END statement. 

Students could be asked to construct Program II using Flow Chart II 
as a guide. Graphs obtained from the selected data could be analyzed to 










FLOW CHART II 



CALCULATE 
Y = AX^ + C 



PRINT 
"N0T A PARAB0LA" 



^ START ^ 






{ 

1 INPUT 


A,C 



PRINT 


YES ^ 


"N0T A PARAB0LA" 






LET X = -5 



PRINT X,Y 



CALCULATE 

9 







LET X 


= X + 1 








FLOW CHART II 
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determine the effect of changes in the values of A and C. 

PROGRAM III 

% 

Dotted lines indicate missing parts of this program. Complete 
the program using Flow Chart III as a guide. 

10 READ A,B, 

20 IF THEN 90 



30 


LET X = -B/(2 * 


, ) 


- 5 




40 


LET Y = A * X * X + 




+ C 




50 


PRINT 








60 


LET 








70 


IF X <= -B/(2 * A) + 


5 ■ 


IHEN 




80 


G0 T0 10 








90 


PRINT ” 




IT 




100 


G0 T0 10 








no 


DATA 1,-4, 4, 1,4, 4,1, 


-8 


,16,1,8,16 




120 


DATA 2, -7, 0,2, 7, 0,-2 




7,7,-2,7,15 




130 


DATA 3. 9, -4. 2, -27. 6 








140 


- - - - — — 














PROGRAM IV 




10 


READ A,B,C 








15 


PRINT A,B,C 








20 


IF A = 0 THEN 130 








30 


LET X = -B/(2 * A) 








40 


LET Y = A * X * X + B 




X + C 




50 


IF Y > C - B * B/(4 * 


A 


* A) + 25 THEN 


10 


60 


IF Y < C - B * B/(4 * 


A 


* A) - 25 THEN 


10 


70 


LET R = X + B/(2 * A) 








80 


PRINT "("x","Y”)" 








90 


PRINT "("X - R","Y")" 








100 


LET X = X + 1 








no 


IF X <= -B/(2*A) + . 


10 


THEN 40 




120 


G0 T0 10 








130 


PRINT ”N0T A PARAB0LA" 






140 


G0 T0 10 








150 


DATA 1,4, 4, 0,2, 3, 1,0 


,0 


,-3,2, -5, -3. 2, -a 


i,14 


160 


END 









The axis of S 3 nnmetry of a parabola whose equation is ax + bx + c = y 
is determined by the equation x = -b/2a. Flow Chart III and Computer Program 
III provide for computation of ordered pairs with integral values of x which 
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lie in the region five units on either side of the axis of S 3 miinetiy of the 

2 

parabola whose equation is ax + bx + c = y. Programs I and II calculated 
values of y for the interval -5 <= x <= 5. Program III calculates values 
of y for the interval (~b/2a) - 5 = x = (-b/2a) + 5. Flow Chart IV and 
Program IV specify limits on the value of y, as well as x. Any limitation 
could be chosen. In Program IV y is calculated for integral values of 
X <= -b/2a + 10 so long as the value of y is less than or equal to the 
parabola's minimum value plus 25 or greater than or equal to the parabola's 
maximum value minus 25. Symmetry of the curve is used to produce pairs 
of points rather than one point at a time. 

In order to be very proficient in programming and working witii 
. computers, it is essential to have a knowledge of mathematics. A computer 
can perforin many operations very rapidly, but it can only do what it is 
told to do. Analyzing a problem carefully in terms of its mathematical 
structure and then writing a program to perform the necessary calculations 
can save a great deal of tedious ccsnputation and also leads to a more thor- 
ough understanding of the mathematics involved. 



o 
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FLOW CHART IV 
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APPLYING THE caiPUTER TO A UNIT ON CIRCLES IN GEOMETRY 

by Karen L. Winquist 

This project is an application of the COTiputer to a unit on circles 
in geometry. The problems and programs presented here are meant to be 
supplementary materials which will enrich and expand a unit on circles. 
While also stimulating the student *s interest in the subject matter, it 
is hoped the material presented will motivate the students to make up 
programs of their own and to investigate many other problems connected 
with circles and also spheres. 

Prior to presenting the problems here, the students should have 
studied the definitions and proofs of the theorems and corollaries in- 
volved in a unit on circles. It is furthermore assxmied that they have 
a knowledge of BASIC language and are fajniliar with programming elementary 
problems and working on the time sharing system of the computer. 
Definitions that Should Be Known for a Unit on Circles 



1. 


Circle 


13. 


Tangent circles externally 


2. 


Sphere 




and internally 


3. 


Center of circle 


14. 


Congruent circles 


4. 


Radius of circle 


15. 


Central angle of a circle 


5. 


Diameter of circle 


16. 


Minor and major arcs 


6. 


Concentric circles 


17. 


Semicircle . 


7. 


Chord 


18. 


Angle inscribed in an arc 


8. 


-Secant 


19. 


Angle intercepting, an. arc - 


9. 


Interior of circle 


20. 


Quadrilateral inscribed in a 


10. 


Exterior of circle 




circle 


11., 


Tangent to . a . circle 


21. 


Quadrilateral circumscribed : 


12. 


Tangent plane to a sphere 




about a circle 


Special Theorems that Should Be 


Known for 


Working with the Problems 



of this Project 

1. The Power Theorem — Giyen a circle C, and a point Q of its exterior — 
let Lji^ be a secant line through Q, intersecting C in" points R. and S. 

Let be another secant line through Q , intersecting C in points 
U and T. Then (QR)(QS) ;= (QU)(QT). ’ , ' 

2,. Given a tangent segment to a circle, and a secant line through 

2 

Q, intersecting the circle in points R and S — then (AR)(QS) = QT . 

3. Let RS and TU be chords of the same circle, intersecting, at Q = — 
then (AR)(QS) = (QU)(QT). 

2 2 2 

4c The graph of the equation (x - a) + (y - b) = r is the. circle 
with center (a, b) and radius r. 

5. Every circle is the graph of an equation of the form 

2 2 

X ' + y + Ax = By + C = 0 . 

2 2 

6. The graph of the equation x +y +Ax+By+C=Ois(l)a circle, 

(2) a point, or (3) the empty set. 
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Program 3^ — Finding the Coordinates of Points of a Circle with its 
Center at the Origin 



2 2 2 

The graph of the equation X + Y - R is a circle with center at 



the origin and radius R. The following program will output coordinates 
of points belonging to the circle for X = 0 up to the length of the radius 
in steps of N. Lines 80 and 90 of the program can be changed for various 
radii and steps. The output is in the first quadrant and will be printed 
in two vertical columns labeled X and Y. It is not necessary to have the 
computer print the values for the other three quadrants as they can be 
figured as follows: 

-X,Y for the second quadrant 
-X,-Y for the third quadrant 
X,-Y for the fourth quadrant 

One example of a result is as follows for a circle of radius 3 
figured where step N = .2. 



X 


Y 


0 


3 


.2 


2.99333 


.4 


2.97321 


.6 


2.93939 


.8 


2.89137 


1. 


2.82843 


1.2 


2.74955 


1.4 


2.6533 


1.6 


2.53772 


1.8 


2.4 


2. 


2.23607 


2.2 


2.03961 


2.4 


1.8 


2.6 


1.49666 


2.8 


1.07703 


3 


2.72958E-4 



Notice that in line 40 when X = 3 with R = 3, Y should = 0; however, 

the computer outputs the square root of 0 - 2.72958E-4, a very small 

number close to zero. This can be eliminated by inserting the following: 

45 LET Z = INT(1000*Y)/1000 
50 PRINT X,Z 

Now all numbers will have three decimal places, giving Y = 0 for X = 3. 

10 READ R,N 



20 PRINT "X","Y" 



30 F0R X = 0 T0 R STEP N 
40 LET Y = SQR(Ri 2-(X*X)) 

50 PRINT X,Y 
60 NEXT X 
70 G0 T0 10 
80 DATA 
90 DATA 
100 END 

Optional for 3 decimal place numbers: 45 LET Z = INT(1000‘<Y)/1000 

50 PRINT X,Z 
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R , N 



X 

11 

o 

c 




1 ^ . 


Y = SQR((R)(R)) - ((X)(X)) 





X , Y 




Program 2 ^ — The Graph of + AX + BY + C = 0 

As pr>oved by a theorem, the equation of all circles is of the 



2 2 9 9 
foim X +Y +AX+BY+C=0; however, conversely, X + Y + AX + BY + C = 

is *not always the equation of a circle. The following program shows that 

the converse is not always true. A, B, and C are the coefficients of X, Y, 

and the constant respectively. Inputing A, B, and C from any equation of • 

the above standard form this program will print what kind of graph would be 

obtained from that equation? As one of the theorems stated, the graph 

will be one of three possibilities : 
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(1) a circle 

(2) a point 

(3) the empty set 

if it is a circle, the computer will print the coordinates of its center 
and the length of its radius. If it is a point, the computer will print 
the coordinates of the point. 

It is assumed here that in studying the proof of the above theorem, 

the student learned the various conditions yielding a circle, point, or 

- 4C 

the empty set frcxn the expression . However, if necessary 

as a review, the expression is derived as follows: 

X2 + + AX + BY + C = 0 

X2 + + AX + BY = -C 

X2 + AX + y2 + BY = -C 

Completing the square: [X^ + AX + (^)^] + [Y^ + BY + (^)^l = -C + 

(X +f)2 + (Y + |)2 - C 

(X + f)2 . (Y . |)2 = + 



As proved by a theorem, the above expression is the equation of a circle 
with its center at the point (-A/2,-B/2) with radius equal to the square 
of 



10 BEAD A,B,C 

20 LET X= -A/2 

30 LET Y= -B/2 

40 LET D=(A*A+B*B-4*C)/4 

50 IF D>0 THEN 110 

60 IF D=0 THEN 90 

70 PRINT "GRAPH IS THE EMPTY SET" 

80 G0 T0 10 

90 PRINT "GRAPH IS THE P0INT" X,Y 
100 T0 T0 10 
no LET R=SQR(D) 

120 PRINT "GRAPH IS A CIRCLE WITH CENTER AT" X,Y 

130 PRINT "THE RADIUS 0F THE CIRCLE IS" R 

140 G0 T0 10 

150 DATA 

160 DATA 

170 END 



Suggested data might be -6, -8, 21 

8 ,- 2,-8 
6, -2, 10 
10,-4,33 



Will be circle with center at (3,4) R=2 
Will be circle with center at (-4,1) R = 5 
Point (-3,1) 

Empty Set 
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A,B, 
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X = -A/2 

\ 

Y = -B/2 

■■■ I ~ 



D = (A^ + -40/4 




Program ^ — Finding the Length of a Tangent Segment Using an Application 
of the Power Theorem 

This program will give the length of a tangent segment AB to a 
circle knowing the coordinates of points D, C, and B of sec^t segment 
DB. DB can be any secant segment to the circle from point B, inter- 
secting the circle at points D and C as shown in Figure A. In the pro- 
gram (XljYl) must be the coordinates of point D, (X2,Y2) the coordinates 
of point C, and (X3,Y3) the coordinates of point B. 

D2 = BC (distance from point B to point C) 

D3 - BD (distance from point B to point D) 

D1 = AB (distance from point A to point B or the length of the tangent) 

If D2 and D3 distances are already known, then they can be substituted 

in lines 20 and 30. 



10 

20 

30 

40 

50 

60 

70 

80 

90 



READ X1,Y1,X2,Y2,X3-,Y3 

LET D2 - 3QR((X2-X3)+2+(Y2-Y3)+2) 

LET D3 = SQR({Xl-X3)+2+(Yl-Y3)+2) 

LET D1 = SQR(D2*D3) 

PRINT "LENGTH 0F THE TANGENT SEGMENT = "D1 
G0 T0 10 
DATA 
DATA 
END 
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PTOgram ^ Showing the Tangent Length to Remain Constant for a Given 
Circle and Given External Point 

In this problem, use the same program as that used in Program 3. 

The computer will show that the length of the tangent segment remains 

the same for various secants drawn to the circle from the given external 
point . 

As shown in the diagram. Figure B, keep point B constant; and find 
the coordinates of D and C for various rotations of the secant. Program 
1 may be helpful in getting the coordinates more exact. Using Program 3, 
place these coordinates in the data statements 70-80. The output will 
show the tangent length remaining fairly constant. It is very difficult 
to find the exact coordinates for all rotations of the secant; and, 
therefore, the tangent lengths will not be equal but will be fairly close, 
showing the student that secant rotations from point B do not affect the 
tangent length. One example that could be used as shown on the diagram 
is as follows (the constant point B being (5,0)): 

70 DATA -3, 0,3, 0,5,0 

71 DATA -.5,2.95,2.75,1.25,5,0 

72 DATA -2.75,1.25,2.95, .33,5.0 

73 DATA -.5,-2.95,2.75,-1.25,5,0 

Lengths of the tangent segments respectively, then, are: 

4 

4.00802 

4.03733 

4.00802 

The 4 is the exact value. 







224 



Program _5 — Discovering the Relationship between Two Circles in the 
Coordinate Plane 

In this problem one is given two circles in a coordinate plane 
and must know the coordinates of the center and radius for each, 

(^IjYl) must be the coordinates of the center of the first circle. 
R1 must be the radius of the first circle » 

(X2,Y2) must be the coordinates of the center of the second circle. 
R2 must be the radius of the second circle. This second circle 
must be the larger circle if one is larger. 

If the two circles have equal radii, then it does not matter which 
circle is first or second. 

In the program E is the distance between the centers of the two 
circles. If the sum of the lengths of the radii is less than E, then the 
circles will not intersect. If, on the other hand, the sum of the lengths 
of the radii is greater than E, then four cases are possible as follows: 

1) Circles are externally tangent if E is equal to the sum of 
the radii. 

2) One circle lies within the other circle with no intersection 
if the larger radius minus E is greater than the other radius. 

3) Circles are internally tangent if the larger radius minus E 
is equal to the other radius. 

4) Circles intersect in two points if the larger radius minus E 
is less than the other radius. 

10 READ X1,Y1,R1,X2,Y2,R2 
20 LET D=(Xl-X2)t2+(Yl~Y2)t2 
30 LET E=SQR(D) 

40 LET F = R1 + R2 

50 IF E = F THEN 90 

60 IF E<F THEN 110 

70 PRINT "CIRCLES D0 N0T INTERSECT" 

80 G0 T0 10 

90 PRINT "CIRCLES ARE EXTERNALLY TANGENT" 

100 GO TO 10 
no LET G = R2-E 
120 IF G - R1 THEN 160 
130 IF G<R1 THEN 180 

140 PRINT "ONE CIRCLE LIES WITHIN OTHER CIRCLE WITH NO INTERSECTION" 

150 G0 T0 10 

160 PRINT "CIRCLES ARE INTERIvTALLY TANGENT" 

170 G0 T0 10 

180 PRINT "CIRCLES INTERSECT IN TWO POINTS" 

190 G0 T0 10 
200 DATA 
210 DATA 
220 END 
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PrHjgrain ^ Plotting a Semicircle and Other Stored Programs 

This problem is a stored program that will plot any semi-circle. 
When the computer asks "new or old," type old. The "old problem name" 

IS XYPLOT***. Upon calling up the old program XYPLOT*** the computer 
will type the following: 

This program sill plot single-valued functions of X, with X 
oil the vertical axis . To use , type : 

10 LET Y = (THE FUNCTION TO BE PLOTTED) 

During running, the program will ask for YMIN and YMAX (the limits 
on the horizontal Y-axis), for XMIN and XMAX (the limits on the 
vertical X-axis), and for DELX, the increment to be used along 
the X-axis. Note: Lines 11-99 of the program may be used as 

desired to express complicated functions. 

To plot a semicircle for line 10 type: 10 LET Y - SQR( (R*R)-(x*X)) 

R is the radius. In a case where, for example, the radius is three: 
YMIN=0 

YMAX=3 

YMIN=-3 

XMAX=3 

DELX IS arbitrary; however, .1 generally works to give a graph closest 
to a semicircle. 

The following are some other stored programs that might be appli- 
cable to a unit on circles in geometry. 

CIRCLE*** Divides circles into n equal parts 
PLOTTO*** Simultaneously plots one to six functions 
TWOPLOT*** Simultaneous plot of two functions 
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Program 7 — Finding the Area of a Circle Incrementing with Narrow 
' Vertical Rectangles 

This program finds the area of a circular region incrementing with 
narrow vertical rectangles. In the program R is the radius of the 
circle, and S is the change in X or the width of each rectangle. For 
each vertical rectangle, then, the area will equal 

length*width 
= Y * X 
= SQR(R^-X^)*(S) 

For more accurate answers the program finds the value of Y halfway 
between any X^ and The program figures the area of each rectangle 

and adds it to the previous area (line 60). This area is 1/4 of the 
area of the circular region, and therefore in line 80 this area is multi- 
plied by 4 to obtain the total area of the circular region. In other 
words, for a circle with center at the origin, it figures the area in 

the first quadrant and then multiplies that area by 4. 

Two programs are listed here. The first will print the area after 
each rectangular area is added. The student can then see the area growing. 
For small increments the print part of this program gets to be quite 
lengthy. Therefore, the second program prints only the final area. In 
this program it is interesting to use S = ,1, .01., .001 respectively to 
see the total area get closer ^d closer to the answer one obtains by merely 
miiltiplying n*R*R, the usual formula for the area of a circular' region. 

10 READ R 
20 READ S 
25 PRINT "X", "A" 

30 LET A = 0 

40 FOR X=S/2 T0 R STEP S 
50 LET Y=SQR((R*R)-(X*X)) 

60 LET A=A+S*Y 
65 PRINT X,A 
70 NEXT X 
80 LET F-4*A 

100 PRINT "TOTAL AREA OF CIRCLE="F 
110 DATA 
120 DATA 
. 150 END 

Program for printing only the total area of the circle. 

Omit lines 25 and 65. 
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